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ABSTRACT 



This thesis examines the pool boiling heat transfer 
characteristics of HCFC-124 and HCFC-124/oil mixtures with up 
to 10% (by weight) miscible alkylbenzene oil. One smooth and 
4 enhanced tubes were tested: a 19 and 26 low integral -fin 
tube (GEWA-K) ; a modified finned tube (TURBO-B) ; and a porous 
coated tube (HIGH FLUX) . The tests were carried out using the 
procedure used for CFC-114 at the same saturation temperature 
of 2.2 ”C. This allowed for direct comparison of the pool 
boiling heat transfer characteristics between the two 
refrigerants . 

The smooth and GEWA-K 19 fin per inch tube performance in 
pure HCFC-124 and HCFC- 124/oil mixtures ranged between 10 to 
50% better than in pure CFC-114 and CFC- 114/oil mixtures for 
all heat fluxes. The HIGH FLUX and TURBO-B tubes were similar 
in performance. With pure HCFC-124, the finned tubes typically 
provided enhancements in the heat transfer coefficient between 
2 and 3 times that of a smooth tube. The HIGH FLUX and TURBO-B 
surfaces typically provided additional enhancements 2 times 
that of the finned tubes. With the addition of oil, the heat 
transfer increased from the smooth and finned tubes but 
decreased from the HIGH FLUX and TURBO-B tubes. The HIGH FLUX 
and TURBO-B tubes therefore exhibited enhancements less than 
the finned tubes at high oil concentrations and high heat 
fluxes . 
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NOMENCLATURE 



A area 

A b Cube outside surface area of active boiling section 

Ac cross-sectional area of tube 

c p specific heat 

D diameter 

D.j_ tube inside diameter 

D 0 outside diameter of unenhanced ends 

D x diameter at the position of the thermocouple 

D 2 diameter of test section to the base of fins 

g gravitational acceleration 

h heat - transfer coefficient 

h fg latent heat of vaporization 

I current 

I s current reading by AC Current Sensor 
k thermal conductivity of liquid refrigerant 

k c thermal conductivity of copper 

L active boiling tube length 

Lu non-boiling length of the test tube 

Nu Nusselt number 

p tube outside wall perimeter 

Pr Prandtl Number 

Q heat - transfer rate from boiling surface 

Q f heat - transfer rate through one non-boiling end 

Q h heat - transfer rate from cartridge heater 

q heat flux 

Ra Rayleigh number 

T temperature 

T avg average wall temperature at the thermocouple location 
T c critical temperature 



x 



T f film temperature 

Tn temperature of the respective (n) thermocouple 
location 

T sac saturation temperature 

T wo outer wall temperature of the boiling test tube 
V voltage across the cartridge heater 

V s voltage reading by AC-DC true RMS converter 

a thermal diffusivity 

/5 volumetric thermal expansion coefficient 

6 uncertainty in measurement and calibration 

6 superheat (Two-Tsat) 

/i dynamic viscosity 

v kinematic viscosity 

p density 

a surface tension 
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I . INTRODUCTION 



In an ongoing effort to eliminate the use of 
chlorof luorocarbons (CFCs) , which have been linked to the 
depletion of the earth's protective ozone layer, the United 
States Navy is looking into using recently developed 
alternative refrigerants. There are 2 principal alternative 
refrigerants: hydrofluorocarbons (HFCs), and 
hydrochlorof luorocarbons (HCFCs) . HFCs contain no chlorine and 
have zero ozone depletion potential (ODP) . HCFCs do contain 
chlorine, but with the attached hydrogen atom to the CFC 
structure, the molecule is less stable and the ozone damaging 
chlorine is released in the lower atmosphere prior to reaching 
the ozone layer. The ODP of HCFC-124 is about 2%, 
significantly less than the 70% ODP for CFC- 114. Figure 1.1, 
provided by Baehr [Ref. 1], displays the advantages of HCFC- 
124 over CFC- 114 with respect to both ozone depletion 
potential (ODP) and global warming potential (GWP) . 

The US Navy is currently using CFC-114 (R-114) in all of 
its larger capacity centrifugal air conditioning units (>100 
tons) . The search for a 'drop-in' replacement has identified 
HCFC-124 (R-124) as a possible alternative to R-114. 'Drop-in' 
here means operating at the same evaporating temperature as 
existing systems. Therefore the operating parameters of R-124 
under similar conditions as present R-114 chiller systems 
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Figure 1.1 Halocarbon Ozone Depletion Potentials and Global Warming Potentials 

relative to R 11 (OOP = 1, GWP = l) 



must be taken into account. For example, the shell must be 
capable of withstanding higher pressures (11 psig for R-124 at 
2.2 °C compared to 0 psig for R-114 at 2.2 °C) . The shutdown 
system pressure is also a significant factor when the system 
is subjected to higher temperatures from a surrounding engine 
room or perhaps a casualty condition. A temperature of 100 °F 
would result in a system vapor pressure of 31 psig for R-114 
compared to 68 psig for R-124. If the heat exchanger shell 
can handle the increased system pressure, hardware 
modifications should be limited to alterations such as 
compressor impeller size and changes to sealing materials to 
maintain leak- tight integrity. 

If hardware modifications can be made, then the other 
major consideration is the heat transfer characteristics of R- 
124 compared to those of R-114. This is not such a simple 
comparison as the two fluids may exhibit varying boiling 
characteristics for each type of enhanced surface. In 
addition, the oil used for each fluid is different since each 
must be miscible in its particular refrigerant and this may 
significantly affect heat transfer performance. 

To evaluate and compare the heat transfer performance of 
R-124 with that of R-114, the following thesis objectives were 
established : 

1. Install a new fiberglass - strengthened evaporator and 
condenser section into existing apparatus to withstand 
the higher saturation pressure associated with R-124. 
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2. Collect pool boiling heat transfer data of pure R-114 and 

R- 114/oil mixtures over a range of increasing and 

decreasing heat fluxes for the newly manufactured tubes 
and check their repeatability with existing R-114 data 
using the newly modified apparatus. 

3. Modify existing data reduction program to account for the 
thermophysical properties of R-124. 

4. Collect pool boiling heat transfer data of pure R-124 and 

R- 124/oil mixtures over a range of increasing and 

decreasing heat fluxes using a smooth tube and four 
enhanced surface tubes. 

5. Study effect of oil on boiling hysteresis for a smooth 
and porous coated tube. 

6. Compare single tube heat transfer characteristics of R- 
124 with previously obtained R-114 data taken by Sugiyama 
[Ref. 2] . 
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II. POOL BOILING OF REFRIGERANTS 



A. POOL BOILING CURVE 

The pool boiling curve, as provided by Bar-Cohen [Ref. 3], 
is shown in Figure 2.1. The applicable regions within the 
curve where single tube pool boiling heat transfer tests of R- 
114 and R-124 refrigerants were conducted are reviewed. 




Figure 2.1 Pool Boiling Curve 
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1. Natural Convection 

This region (a-b) occurs during the increasing heat 
flux phase where there are no vapor bubbles being produced at 
the tube surface. In this region, heat transfer is due solely 
to single-phase free convection effects. Several correlations 
are available to predict heat transfer performance in this 
region, including the Churchill and Chu [Ref. 4] correlation 
shown below in equation (1) . 

h=k/D 0 [ . 6+ . 3 87 (g-/3-e-D 0 3 /u - a) 1/6 / [1+ (0 . 559 /Pr) 9/l6 ] 8/27 ] 2 ( 1 ) 



2. Incipient Boiling 

The transition from the natural convection region to 
the nucleate boiling region is the subject of numerous 
studies. The superheat (T wall -T sat ) at which bubble inception 
(point b) occurs is referred to as the onset of nucleate 
boiling (ONB) or incipient boiling. Considerable reduction in 
tube wall temperature occurs at this point, resulting in a 
considerable reduction in superheat noted as AT ex in Figure 
2.1. This 'jump' in performance only occurs during 
increasing heat flux and is referred to as the thermal 
overshoot. A hysteresis effect occurs when decreasing heat 
flux data (from point c to point a) are compared to the 
initial ascending path (a-b-c) . Several factors believed to 
influence the thermal overshoot at the ONB for highly- wetting 
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liquids were studied by Bar-Cohen [Ref. 3]. Probable 
influential factors included: vibration, surface tension, 
enhanced commercial surfaces, pressure, time between test 
runs, oil mixture, temperature attained at end of decreasing 
heat flux phase of the cycle, and variation with cool-down 
procedure. Because of these numerous influences on the boiling 
incipience superheat, past studies, including Bar-Cohen [Ref. 
3], reveal highly variable and widely distributed superheats 
at ONB . Furthermore, the process of first nucleation is 
inherently stochastic, leading to variations in behavior from 
one day to the next . 

3. Nucleate Boiling Region 

Once nucleate boiling commences, during increasing 
heat flux (point c) , sustained improvement in heat transfer is 
realized throughout the entire fully-developed nucleate 
boiling region (c-d) . As vapor bubbles form at nucleation 
sites and separate from the surface, dramatic increases in 
heat transfer are attained. There are several factors 
affecting where this region is located on the pool boiling 
curve, the most important of which is the number of active 
sites on the surface, as may occur on different boiling 
surfaces. Thome [Ref. 5] provides detailed insight into 
enhanced boiling heat transfer from various commercially 
available surfaces. For plain tubes, several correlations have 
been proposed to predict boiling performance. Two such 
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correlations were provided by Stephan and Abdelsalam [Ref. 6] 
and Chongrungreong and Sauer [Ref. 7] . 

B. SINGLE TUBE BOILING EXPERIMENTS AT 

NAVAL POSTGRADUATE SCHOOL 

Heat transfer measurements were collected by McManus 
[Ref. 8] for single tube pool boiling of R-114 at a saturation 
temperature of 13.8°C with up to 6% oil concentration. The 
boiling tubes in these tests were heated internally by warm 
water. HIGH FLUX, TURBO- B and smooth tubes were tested in this 
manner. McManus [Ref. 8] reported enhancements (compared to 
the smooth tube) in pure R-114 as high as 14.6 and 6.4 for the 
HIGH FLUX and TURBO- B tubes respectively, while performance 
decreased to enhancements of 7.0 and 4.9 respectively with a 
6% oil concentration. 

Pool boiling heat transfer coefficients in pure R-114 and 
R- 114/oil mixtures up to 10% (by weight) were measured by 
Wanniarachchi et al [Ref. 9] for electrically-heated smooth 
and enhanced tubes. These tests were run at a saturation 
temperature of 2.2°C. At 30 kW/m 2 Wanniarachchi et al [Ref. 9] 
reported enhancement results for the HIGH FLUX, Thermoexcel - E , 
Thermoexcel -HE, and GEWA-T tubes of 9.1, 8.2, 6.8, and 4.4 
respectively. Also reported were the enhancements for these 
tubes with 3%, and 10% oil concentration. 

Sugiyama [Ref. 2] conducted extensive pool boiling tests 
with R-114 and R- 114/oil mixtures using a large variety of 
enhanced surface tubes. For the re-entrant cavity tubes with 
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pure R-114, Sugiyama [Ref. 2] reported enhancements of 6.4 and 
6.2 for the HIGH FLUX and TURBO-B tubes respectively at 35 
kW/rri . With 10% oil these enhancements were reduced to 3 . 8 and 
5.6 respectively. This degradation in performance was 
partially attributed to the clogging of surface pores with 
oil. For GEWA-K 26 and 40 fpi tubes, enhancements of 2.4 and 
3.1 were reported with pure R-114 at 35 kW/rri. With 10% oil 
these same tubes were surprisingly reported to have 
enhancements of 3.2 and 4.6 respectively. This indicated an 
improved heat transfer performance with addition of oil in the 
case of finned tubes, believed to be caused by foaming action 
created within the pool that was most prevalent at higher oil 
concentrations (3-10%) and heat fluxes (>30 kW/m 2 ). 

The effect of oil on the incipience of nucleate boiling 
for pure R-114 and R- 114/oil mixtures was studied by Memory 
and Marto [Ref. 10] with smooth, finned and re-entrant cavity 
tubes. They observed a significantly lower heat flux at 
incipience for the HIGH FLUX and TURBO-B tubes compared with 
the smooth and finned tubes in pure R-114. However, consistent 
delay in the boiling incipience of the HIGH FLUX and TURBO-B 
tubes was exhibited in the presence of 3 and 10% oil 
concentrations. At these oil concentrations, boiling 
incipience for the smooth and finned tubes showed no 
dependence with oil . Also emphasized was the inherent 
randomness and variation connected with temperature overshoot 
associated with boiling incipience. 
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C. COMPARISON OF R-114 TO R-124 



1. Property Comparison of R-114/R-124 

Table 1 is provided to show direct comparison of the 
properties of R-114 and R-124 at T sat =2.2°C 



Table 1. R-114/R-124 PROPERTY COMPARISON 



Property 


R-114 


R-124 


Pul/ (kPa) 


96 


177 


P (liq) / (kg/m 2 ) 


1526 


1429 


p (vapor) / (kg/m 3 ) 


7.5 


11.3 


c p (liq) / (J/kg- K) 


932 


1060 


h fg / (kj/kg) 


134 . 7 


158 . 9 


k(liq) / (W/m-K) 


0 . 0701 


0 . 0744 


M (liq) / (g/m- s) 


0.4449 


0.3375 


cr/ (N/m) 


0 .0136 


0 . 0129 



2 . Need for Oil Mixtures 

The oil mixtures used are to simulate the actual 
environment in refrigeration systems which can collect lube 
oil within a unit's flooded evaporator to concentrations as 
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high as 10%. The oils used for R-124 and R-114 are both 
miscible within their respective refrigerants. However, with 
R-114, a York-C mineral oil was used, while with R-124, an 
alkylbenzene oil was used. 

3. Need for R-124 Data 

There is presently no data on single tube pool boiling 
heat transfer with R-124 or R-124/oil mixtures. Tests must 
therefore be conducted with R-124 to parallel R-114 results 
obtained to date. 
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III. DESCRIPTION OP EXPERIMENTAL APPARATUS 



A. SYSTEM DESCRIPTION 

Details of the relocation and early modifications to the 
original apparatus are given by Sugiyama [Ref. 2] and Lake 
[Ref. 11] . To accommodate the increased operating pressures 
associated with R-124, stronger evaporator and condenser units 
were installed. The data reduction program DRP71 was modified 
(DRPGB) to include the thermophysical properties of R-124. 

A schematic diagram representing the experimental 
apparatus is shown in Figure 3.1 and comprises the following 
components: 1) an evaporator tee used to boil the refrigerant; 
(2) a condenser tee used to condense the refrigerant; (3) a 
reservoir for liquid refrigerant; (4) a ref rigerant/oil 
subsystem; ( 5 ) a refrigeration cooling support system; (6) a 
vacuum pump; (7) a data acquisition and instrumentation 
system; (8) an aluminum/plexiglas framework housing items (1) 
to (4) . 

The general operation of the system was closed loop with 
vapor generated in the evaporator by the heated tube located 
within the refrigerant pool of the evaporator. The vapor 
passed through an aluminum 'L- shaped' tube to the condenser 
where the vapor was condensed by a chilled water/ethylene 
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Figure 3.1 Schematic of Apparatus 
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glycol mixture circulated through the condenser coil by an 8 
gpm turbine type positive displacement pump. The R-114 or R- 
124 condensate was then gravity fed back to the evaporator 
through a return line. R-502 and R-12 refrigeration systems 
maintained the water/ethylene glycol mixture in the sump 
between -12 and -18 degrees Celsius. Sump cooling was 
accomplished using a counter flow heat exchanger by 
recirculating the sump coolant (using an identical positive 
displacement pump) over the R-502 side of the heat exchanger. 
To further maintain this coolant temperature, the R-12 system 
was located directly within the water/ethylene glycol sump. 

The addition of oil to the refrigerant within the 
evaporator was accomplished by using a graduated cylinder to 
measure the precise quantity. This simply drained into the 
evaporator by gravity. 

Accurate horizontal positioning of the boiling tube was 
critical and accomplished by supporting each end of the tube 



by Teflon 


inserts 


. The 


tube was 


sealed using 0- rings. 


For 


safety, a 


relief 


valve 


on the 


aluminum pipe between 


the 



evaporator and condenser was set at a gage pressure of 49 psi. 

B. BOILING TEST SECTION 
1. Evaporator 

The evaporator was a Pyrex glass 'T' section coated 
with a continuous winding of fiberglass filaments impregnated 
with a modified epoxy resin. This upgraded evaporator shown in 
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Figure 3.2 was designed for a working pressure of 50 psi gage 
and initially hydrostatically safety tested to 75 psi gage. 
The original glass was designed for a working pressure of 30 
psi gage. The inside surface was the same as before and 
maintained the advantages of minimizing nucleate boiling at 
the surface as well as being corrosion resistant. The two ends 
of the evaporator were fitted with a cast-iron flange and 
gasket assembly as shown in Figure 3.2. They mated with 
aluminum flanges at each end of the evaporator, which in turn 
contained the thermocouple housings, oil entry connection and 
Teflon inserts containing the evaporator tube. 

2 . Evaporator tubes 

Figure 3.3 shows a schematic of a boiling tube. These were 
accurately held in place using Teflon inserts and sealed by 0- 
rings . A 4-stud strong-back secured the Teflon inserts inside 
the aluminum flanges. 

Five different tubes were tested within the apparatus. 
Dimensions of each tube are provided in Table 5 of Appendix C 
while an expanded view of the enhanced surfaces are shown in 
Figure 3.4. The tubes used were: 

1 . Smooth tube 

2. GEWA-K (19 fpi ) finned tube 

3. GEWA-K (26 fpi) finned tube 

4 . HIGH FLUX (porous coated) enhanced tube 

5. TURBO-B (deformed surface) enhanced tube 
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I 

a) Corning Pyrex Glass Evaporator (D x Q. = 
L = 178 mm, = 127 mm) 






d, 






MO 











insert 



b) Cast Iron Flange and Gasket (d^= 190mm, 
L 1 = 14mm, A = 21° ) 



402x51 mm 



d 0 = 210mm 



Figure 3.2 Sketch of Pyrex Glass Vessels and Flanges 




Figure 3.3 Schematic of boiling tube 
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Figure 3.4 Surface Profiles of Enhanced Tubes 
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The heater used within each boiling tube was nominally a 
1 kW, 240-Volt stainless - steel cartridge heater, measuring 
6.35 mm in outside diameter and 203.2 mm in length. The actual 
heated length of the heater was 190 mm. The HIGH FLUX and 
GEWA-K 26 fpi tubes were previously manufactured while the 
TURBO-B, GEWA-K 19 fpi and smooth tubes required assembly 
prior to testing with R-124 or R-114. The manufacturing 
process of these tubes is described below. 

The heater was placed inside a copper sleeve (Figure 3.5 
provides dimensions of smooth tube sleeve) with a tight 
mechanical fit of .001 inch clearance. The copper sleeve was 
machined with 6 grooves (1.27mm x 1.27mm) on its outer 
diameter of various lengths. Thermocouples were laid in each 
of these grooves and the edges of the grooves were peened over 
with a punch. Any resulting imperfections were smoothed out 
with 250 grit emery paper. The thermocouple wires and heater 
leads all exited the copper sleeve at the same end. The copper 
sleeve containing the heater and thermocouples was then placed 
on 2 aluminum V-blocks and the outer surface of the sleeve was 
covered with a flux and solder paste (50/50 tin/lead) ; great 
care was taken to ensure complete coverage of the thermocouple 
wires lying in the grooves. The thermocouple and heater leads 
were then passed through the enhanced boiling tube (the outer 
diameter of the copper sleeve was machined to fit tightly into 
each respective tube) . To prevent shorting between the heater 
and outer tube shell, a high temperature RTV silicone rubber 
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(DIMENSIONS IN mm) 




Figure 3.5 Schematic of copper sleeve 
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was applied to the heater and thermocouple leads at the point 
where they leave the sleeve. The copper sleeve was then 
tapped gently into the tube and positioned accurately so as to 
line up precisely with the 8 inch active length of the 
enhanced surface. This enhanced surface length was different 
from the actual heated length which was only 7.5 inches. 

The entire assembly was then placed horizontally in an 
oven and heated to between 400°F and 500°F. This was the 
temperature 'window' at which the flux melted, but the 
thermocouple leads did not. The temperature of the tube was 
monitored by attaching the thermocouples to a simple 
temperature indicator. After stabilizing the oven temperature 
between these temperatures and ensuring that the solder had 
melted, the tube was taken out and quenched with a damp cloth 
at one end (as close to the end as possible) where the 
thermocouple and heater leads protrude. The oven was then 
placed vertically and the tube reinserted from below keeping 
the lower end quenched while holding both tube and cloth with 
a pair of channel locks. When the temperature increased to 
>400°F , the solder again began to melt. Additional 40/60 tin- 
lead solder was then applied by hand from the end which 
protrudes from the top of the oven. Remelting and addition of 
extra solder in this way ensured a uniform circumferential 
layer of solder and helped to prevent formation of air 
pockets. The tube was then removed from the oven and quenched 
over its whole length with a damp cloth while in a vertical 
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position. It was then left to cool in air for 24 hours. This 
procedure differed from earlier manufactured tubes in that the 
previous copper sleeves had 4 thermocouples leaving each end 
of the sleeve. Trying to repeat this previous procedure 
accurately led to large discrepancies in the thermocouple 
readings (>20°C) as heater power was increased. It was felt 
this was due to air pockets created in the fabrication 
procedure (the tube could never be placed vertically in the 
oven with leads protruding from both ends) . The above 
procedure was much simpler and led to excellent uniformity in 
the thermocouple values at all heater powers. 

C . CONDENSER SECTION 

The condenser was also a fiberglass coated Corning Pyrex 
glass Tee identical in size and shape to the evaporator shown 
in Figure 3.2. The refrigerant was condensed on a 3/8 inch 
outside diameter copper tube that was coiled within the 
condenser. The condenser was mounted vertically. An O-ring 
sealed aluminum flange mounted at the top of the condenser 
contained not only the water ethylene/glycol coolant 
inlet/outlet, but also a 3/8 inch outside diameter copper tube 
vent, which was connected either to a vacuum pump via valve 
Vll or to a charging cylinder via valve V10. The bottom of the 
condenser was also capped with an O-ring sealed aluminum 
flange and contained the gravity feed return line to the 
evaporator via valve V5 . The refrigerant vapor from the 
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evaporator entered the condenser through an aluminum L- shaped 
tube, connected to the condenser by O-ring and gasketed 
flange. A pressure gage ranging from 30 inches mercury vacuum 
to 150 psig was mounted on the L- shaped tube; this was used 
only as a rough visual check of pressure within the facility. 

D. OIL ADDITION SECTION 

To provide oil addition into the refrigerant, a 
cylindrical aluminum reservoir, 6 inches in diameter and 6 
inches high, and a graduated cylinder 355 mm in length and 
25.4 mm in diameter is used as shown in Figure 3.1. The oil 
used with R-124 was a miscible alkylbenzene oil while the oil 
used with R-114 was a York-C mineral oil, also miscible. The 
graduated cylinder was filled from the reservoir via valve V2 , 
while the reservoir was replenished through valve V3 at the 
top. Oil entered the evaporator from the graduated cylinder by 
gravity feed via valve VI. Oil concentrations (by weight) of 
0,1, 2, 3, 6, and 10% were used for each of the tubes. 

E. COOLING SYSTEM 

1. Water /Ethylene- Glycol Coolant Tank 

At the base of the apparatus, a 5.35 cubic foot tank 
stored 30 gallons of ethylene -glycol/water mixture. 
Recirculation and discharge ports were connected to the tank 
together with a thermocouple to monitor sump temperature. 
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2 . Refrigeration Plants 

To provide the cooling of the water/ethylene glycol 
mixture, a 1/2 ton R-502 and 1/4 ton R-12 refrigeration system 
were installed near the sump as seen in Figure 3.1. Both 
plants contained a hermetically sealed compressor assembly, 
an-air cooled condenser, receiver, filter-dryer, pressure 
regulator, temperature control switch, and thermostatic 
expansion valve. The R-502 plant (Figure 3.6) used a counter 
current heat exchanger while the R-12 plant used a coiled 
evaporator placed within the sump. The sump temperature was 
maintained at approximately -17 °C when both these systems 
were used. Each plant was controlled by a temperature control 
switch and thermostatic expansion valve. 

3 . Pump and Control Valve 

An 8 gpm, positive displacement pump was installed 
with a 1 inch diameter PVC pipe connected to the 
water/ethylene glycol sump. The discharge of the pump was 
piped to the condenser via valve VC as shown in Figure 3.1. A 
by-pass valve V9 was placed upstream of the control valve VC 
to avoid overloading the pump when control valve VC was closed 
and provides for mixing of the return coolant from the 
condenser. At high heat fluxes, the by-pass valve needed to be 
closed down to increase flow of coolant through the condenser. 
An identical pump discharged the ethylene glycol/water coolant 
through the R-502 counter- current heat exchanger. 
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Figure 3.6 



R-502 Refrigerant plant 
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F. REFRIGERANT RESERVOIR 

The reservoir used to store the R-114 or R-124 was made 
from an aluminum cylinder, 9 inches in diameter and 10 inches 
in height. A sight glass attached on the reservoir allowed for 
the liquid level of refrigerant to be monitored. The reservoir 
was connected to the ' L' shaped tube through valve V7 and to 
the liquid line into the evaporator through valve V6 as shown 
in Figure 3.1. 

G. FRAME 

To house the entire assembly, an aluminum frame was 
constructed measuring 1070.1 mm x 509.8 mm x 610.1 mm. With a 
view to keeping the apparatus in an isothermal atmosphere, as 
well as for safety considerations, the sides were covered with 
12 . 7 mm thick plexiglas with hinges mounted on the 2 sides to 
provide easy access to the assembly. The top of the frame was 
covered by plywood while the bottom was covered by an aluminum 
plate. The entire frame rested on the coolant sump tank. The 
front side of plexiglas had small holes to accommodate the 
valve stems of valves VI through V8 , so they can be operated 
while maintaining the integrity of the enclosed frame. 

H. INSTRUMENTATION 

1 . Power Measurement 

The boiling tube heater was powered by a 240 volt AC 
source, adjusted to 0-220 volts and 0-5 amps by a variac 
control. Power input to the heater was measured with an AC 
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Figure 3.7 



Data acquisition/control unit 
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current inductive sensor and a voltage sensor. The output was 
a proportional signal in volts having been converted through 
an AC-DC true R.M.S. converter. Both the AC current sensor and 
the AC-DC true R.M.S. converter were connected to the data 
acquisition/control unit as shown in Figure 3.7. 

2 . Temperature Measurement 

The primary means of obtaining temperatures was the 
use of 30 gage copper- constantan thermocouples. Specific 
temperatures recorded were from 6 thermocouples in the boiling 
tube (8 in the earlier manufactured tubes) , 3 thermocouples in 
the evaporator pool (housed within separate sealed 
thermocouple wells within the evaporator pool), and one 
thermocouple to monitor the water/ethylene glycol sump tank. 
A Hewlett-Packard 3497A data acquisition system controlled by 
a Hewlett-Packard 9826 computer read these thermocouple 
measurements . 

I. DATA ACQUISITION AND REDUCTION PROGRAM 

Data collection and reduction program DRP71 , used by 
Sugiyama [Ref. 2) for use with R-114 was modified for use with 
R - 124 . The thermophysical properties of saturated R-124 were 
obtained from Ref Prop software [Ref. 12] and were used to 
generate equations for use in the calculations of the 
reduction program. Appendix A contains plots of each property 
versus temperature and the respective equation. The modified 
program was renamed DRPGB . 
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The Hewlett-Packard 9826 computer reduced all the 
thermocouple, current, and R.M.S. voltage readings using the 
DRPGB program. The program was run through the keyboard 
allowing regulation of the data acquisition unit. A printout 
of the reduced data was provided by a Hewlett-Packard Inkjet 
printer. Table 2 lists the channel assignments for the sensor 
inputs to the data acquisition system. 



J. PROCEDURE FROM KEYBOARD 

The following procedure outlines the procedure used to 
collect and process data from the Hewlett Packard data 
acquisition system: 



1. Selected the 'taking data' option after initial 

loading of program DRPGB. 

2. Selected the refrigerant being tested, tube heating 

mode, (ie electrically or water heated) arid 

thermocouple type from DRPGB. 

3. Assigned name to file of data to be processed and 

stored . 

4. Identified and entered any defective thermocouples by 

channel # 1-8, noting that the new tubes have only 6 
thermocouples and channels 7 & 8 were always 

' defective' . 

5. Selected test tube type. 

6. Set desired saturation temperature (2.22°C) for 

present tests. 

7. Set desired heat flux desired using variac adjustment. 

8. Attained desired saturation temperature by adjusting 
flow of coolant through condenser coils with control 
valve VC. 
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9. Once saturation temperature was achieved, waited 5 
minutes for steady state conditions prior to taking 
data . 

10. Prompted data acquisition unit to scan all the 
channels listed in Table 2. All channel readings were 
made in volts and stored in user specified fields. 

11. Calculated temperature and power from these voltages. 

12. Calculated heat transfer rate from cartridge heater. 

13. Calculated the average wall temperature of the test 
tube and calculated wall superheat (Twall- Tsat) . 

14. Calculated physical properties of R-124 using property 
correlations from RefProp [Ref. 12] at film temperature 
(Twall +Tsat ) /2 . 

15. Calculated the natural convection heat - transfer 
coefficient of R-124 from the unheated ends of the 
test tube. 

16. Computed heat losses from unheated ends of the test 
tube . 

17. Calculated heat flux from heated length of the test 
tube . 

18. Calculated the heat - transfer coefficient from the 
heated length of the test tube. 

19 . Stored heat flux and wall superheat for each data set 
in user specified fields. 

20. Plotted data on available software. 



Sample calculations of the above procedure are given in 
Appendix C. 
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TABLE 2. HP 3497A CHANNEL ASSIGNMENTS 



Channel 


Channel Assignments 


00-07 


tube wall temperatures 


08-10 


pool liquid temperature 


11 


sump temperature 


20 


RMS voltage 


21 


current sensor 
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IV. EXPERIMENTAL PROCEDURE 



A. ASSEMBLY PREPARATION 

1. Vacuum test of the apparatus 

After placing the desired boiler tube m position 
within the evaporator, the system was evacuated down to a 
pressure of 29 inch Hg . vacuum using a portable vacuum pump 
via valves Vll and V8 in Figure 3.1. The system was then 
isolated for approximately 30 minutes to check for leak- tight 
integrity (indicated by any significant vacuum loss). If a 
leak was detected, a positive air pressure (20 psig) was 
placed on the system via valves V10 and V8 . A soap/water 
solution was applied on the sealing areas to identify any 
leaks (if there is a leak, bubbles are seen) . If identified, 
the leak source was corrected and the system integrity 
retested. 

2. Charging apparatus with refrigerant 

After system integrity was ensured, the evaporator was 
filled with refrigerant to a level of 20 mm above the test 
tube. The refrigerant was drained by gravity from the 
reservoir into the evaporator through valves V6 and V5 . Valve 
V7 was also open to equalize pressure in the system. The 
reservoir was then isolated from the system with the closure 
of valves V5 through V7 . 
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3. Degassing and Data Acquisition Channel Check 

Prior to connecting all thermocouple leads, the 
resistance between the heater leads and outer tube wall was 
checked with a simple multimeter. Readings were usually 
greater than 1 Mft . Each channel's output was checked using 
program SETUP71. Any problem thermocouples were either fixed 
or excluded from the calculations in program DRPGB . Before 
energizing the heater and starting the tests, the evaporator 
pool was slowly cooled to a saturation pressure of 2.22 °C by 
allowing a low coolant flow through the condenser. Once 
completed, the heating surface flux was set to approximately 
90 kW/m 2 to degas the refrigerant; any non- condensible gases 
were collected in the condenser and vented off using the 
vacuum pump for approximately one minute. The apparatus was 
then secured for a period of approximately 8 hours prior to 
operation to allow the surface to cool and become fully wetted 
in the refrigerant. 

B . OPERATIONAL PROCEDURE 

The following procedure was used to operate the facility 
(this follows closely the procedure used by Sugiyama [Ref .2] ) : 

1. The R-502 and R-12 refrigeration units were utilized 
until the temperature of the water- ethylene glycol sump 
reached a minimum of - 8 °C before starting tests (with 
continued cooling to -17 °C) . 

2. The data acquisition unit, computer and variac control 
panel were switched on. 
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3. The computer program SETUP was loaded and run. All data 
acquisition channels were rechecked. The refrigerant 
pool temperature was then cooled down slowly to 2.2 °C 
by circulating coolant through the condenser. 

4. The data acquisition/reduction program DRPGB was loaded 
and run . 

5. The desired heat flux setting for the test tube was 
input to the program DRPGB and the vanac was adjusted 
to attain that heat flux. Starting heat flux was 
between 500 and 600 W/m 2 , 

6. Control valve VC was adjusted to control coolant flow 
through the condenser to maintain a constant saturation 
temperature of 2.22 °C + o.l °C. 

7. At each heat flux setting, the thermocouple values were 
allowed to stabilize for about 5 minutes prior to 
taking data. 

8. The following data were then taken and used in the 
acquisition/reduction program DRPGB: pool temperature; 
test tube wall temperatures; sump temperature; current 
and voltage sensor readings. 

9. Data were taken at each desired heat flux setting and 

always at the same saturation temperature (the coolant 
flow in the condenser was adjusted to ensure this) . At 
each heat flux setting, program DRPGB generated a 
permanent printout listing: wall temperatures of the 

test tubes; refrigerant pool temperatures; sump 
temperature; wall superheat; heat - transfer coefficient 
and the corrected heat flux (ie. to account for heat- 
transfer lost from each end) . 

10. For each heat flux setting, steps 6 through 9 were 
repeated. The superheat (Twall-Tsat) and the respective 
heat flux were plotted on a log-log scale. 

11. Table 3 provides a summarized listing of all data runs. 
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TABLE 3. LISTING OF DATA RUNS 



Data File 


Refrigerant 


Tube Type 


Purpose 


1 D0821SM0 


R- 114 


Smooth 


repeatability 


j D0824SM0 


R- 114 


Smooth 


repeatability 


D1015SM0 


R- 124 


Smooth 


Data 


D1017SM0 


R- 124 


Smooth 


Data 


D122#SM0 


R- 124 


Smooth 


Data 


D1222SM0 


R- 124 


Smooth 


Data 


D1221SM0 


R '124 


Smooth 


Data/samp calc 


D1220SM0 


R- 124 


Smooth 


uncert analys 


D1219SM0 


R- 124 


Smooth 


Data 


D1217SM0 


R '124 


Smooth 


Data 


D1016SM4 .44* 


R- 124 


Smooth 


Data Tsat4.4°C 


D1114SM0 


R- 124 


Smooth 


Data 


D1114SM# 


R- 124 


Smooth 


Data 


D1114SM2 


R '124 


Smooth 


Data 


D0108SM3 


R '124 


Smooth 


Data 


D0107SM3 


R- 124 


Smooth 


Data 


D0106ASM3 


R '124 


Smooth 


Data 


D0106SM3 


R- 124 


Smooth 


Data 


D0105SM3 


R- 124 


Smooth 


Data 


D0104SM3 


R '124 


Smooth 


Data 


D1227SM3 


R- 124 


Smooth 


Data 


D1116SM3 


R '124 


Smooth 


Data 


D1116SM6 


R- 124 


Smooth 


Data 


D0112ASM10 


R- 124 


Smooth 


Data 


D0112SM10 


R- 124 


Smooth 


Data 


D0111ASM10 


R- 124 


Smooth 


Data 


D0111SM10 


R- 124 


Smooth 


Data 
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Continuation of Table 3 . 



Data File 


Refrigerant 


tube type 


Purpose 


DO 110SM10 


R- 124 


Smooth 


Data 


D0109ASM10 


R- 124 


Smooth 


Data 


DO 109SM10 


R- 124 


Smooth 


Data 


D1117SM10 


R- 124 


Smooth 


Data 


D1O1019O 


R- 114 


GK-19 


Data 


D1010190 


R- 114 


GK-19 


Data 


D0901190 


R- 114 


GK-19 


Data 


D0907193 


R- 114 


GK-19 


Data 


D09091910 


R- 114 


GK-19 


Data 


D1105190 


R- 124 


GK-19 


Data 


D1118190 


R- 124 


GK-19 


Data 


D111819# 


R- 124 


GK-19 


Data 


D1118192 


R- 124 


GK-19 


Data 


D1119193 


R- 124 


GK-19 


Data 


D1120196 


R- 124 


GK-19 


Data 


D11211910 


R- 124 


GK-19 


Data 


D0803260 


R- 114 


GK-26 


Repeatability 


D1021260 


R- 124 


GK-26 


Data 


D1122260 


R- 124 


GK-26 


Data 


D1122261 


R- 124 


GK-26 


Data 


D1122262 


R- 124 


GK-26 


Data 


D1123263 


R- 124 


GK-26 


Data 


D1123266 


R- 124 


GK-26 


Data 


D11252610 


R- 124 


GK-26 


Data 


D11242610 


R- 124 


GK-26 


Data 


D1004TB0 


R- 114 


TURBO -B 


Data 


1 D1006TB0 


R- 114 


TURBO - B 


1 Repeatability \ 
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Continuation of Table 3. 



Data File 


Refrigerant 


Tube 


Type 


Purpose 


D1002GYX0 


R- 114 


GEWA- 


YX 


Data 


D1110TB0 


R- 124 


TURBO 


-B 


Data 


D1110TB0 


R- 124 


TURBO 


-B 


Data 


D1111TB1 


R- 124 


TURBO 


-B 


Data 


D1111TB2 


R- 124 


TURBO 


>-B 


Data 


D1112TB3 


R- 124 


TURBO 


-B 


Data 


D1112TB6 


R- 124 


TURBO 


>-B 


Data 


D1113TB10 


R- 124 


TURBO 


-B 


Data 


D0119AHF0 


R- 124 


HIGH 


FLUX 


Data 


D0113BHF0 


R- 124 


HIGH 


FLUX 


Uncert Analvs 


D0118HF0 


R- 124 


HIGH 


FLUX 


Data 


DOliaHFO 


R- 124 


HIGH 


FLUX 


Data 


D0119AHFS 


R- 124 


HIGH 


FLUX 


Data 


D0113HF0 


R- 124 


HIGH 


FLUX 


Data 


D1108HF0 


R- 124 


HIGH 


FLUX 


Data 


D1126HF0 


R- 124 


HIGH 


FLUX 


Data 


D1126HF1 


R- 124 


HIGH 


FLUX 


Data 


D112 4HF# 


R- 124 


HIGH 


FLUX 


Data 


D0125HF3 


R- 124 


HIGH 


FLUX 


Data 


D012CHF3 


R- 124 


HIGH 


FLUX 


Data 


DO 12 1CHF3 


R- 124 


HIGH 


FLUX 


Data 


D0121KF3 


R- 124 


HIGH 


FLUX 


Data 


DO 12 0AHF3 


R- 124 


HIGH 


FLUX 


Data 


D0120HF3 


R- 124 


HIGH 


FLUX 


Data 


DO 119AHF3 


R- 124 


HIGH 


FLUX 


Data 


D1128HF3 


R- 124 


HIGH 


FLUX 


Data 


D1128HF6 


R- 124 


HIGH 


FLUX 


Data | 
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Continuation of Table 3 . 



Data File 


Refrigerant 


Tube 


d) 

s 


Purpose 


D1130HF10 


R-124 


HIGH 


FLUX 


Data 


D0128HF10 


R-124 


HIGH 


FLUX 


Data 


D0127CHF10 


R-124 


HIGH 


FLUX 


Data 


DO 12 7BHF10 


R-124 


HIGH 


FLUX 


Data 


D0127HF10 


R-124 


HIGH 


FLUX 


Data 


DO 12 6BHF10 


R- 124 


HIGH 


FLUX 


Data 


D012 6AHF10 


R-124 


HIGH 


FLUX 


Data 


D0126HF10 


R-124 


HIGH 


FLUX 


Data 



The tabulated data was filed using the following file name 
system: 

Example (D0126AHF10) : 

1. First letter D simply refers to data 

2. The following 4 characters are always the date; in this 
case Jan 26. 

3. The next letter represents order of the day test was 
completed on the respective date; 

no letter- 1st run 
A- 2nd run 
B- 3rd run 
etc . 

4. The following 2 letters represented the tube used; 

SM= smooth 
19= 19 fpi GEWA-K 
26= 26 fpi GEWA-K 
HF= HIGH FLUX 
TB= TURBO -B 

5 . The last 1 or 2 numbers represent the percentage of 
oil; in this case 10%. 

* this one run was pure R-124 tested at 4.4 °C vice 2.2 °C 
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V. RESULTS AND DISCUSSION 



A. REPRODUCTION OF EXPERIMENTAL DATA 

The first tests conducted were to verify the 
reproducibility of the single tube data taken by Sugiyama 
[Ref. 2] with R-114. A GEWA-K 26 fin per inch (fpi) tube 
comparison with R-114 is given in Figure 5.1 and shows very 
similar results with previous data. Uncertainty bands for both 
wall superheat and heat flux are given on the figure. 
Consequently, at low heat flux (< 1000 W/m 2 ) , there is more 
scatter due to the lower wall superheats and hence higher 
uncertainty. Overall, the agreement between the two data sets 
is satisfactory. 

In Figure 5.2, an increasing heat flux data run using 
a smooth tube is compared with Sugiyama' s [Ref. 2] data. It 
can be seen that in the low heat flux, natural convection 
region and the high heat flux, nucleate boiling region the 
data agreement is good. In the mid-range of heat flux, (3 
kW/m 2 < q < 15 kW/m 2 ) , the data varies somewhat, depending upon 
how the tube nucleates. Curve #1 nucleates at a relatively low 
wall superheat (12.5 °C) while curve #2 (both taken by 
Sugiyama [Ref.l]) nucleates at a high wall superheat of about 
35 °C (point C) . The present data (curve #3) exhibits partial 
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Tube (increasing flux) 
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nucleation (see point A) where 3 of the six thermocouples had 
a low wall superheat (< 14 °C) and the remaining 3 had a wall 
superheat (> 18 °C) . Observing the tube surface, it could be 
clearly seen that half the tube was boiling while the other 
half was in the natural convection mode. At point B, all six 
thermocouples were within 0.9 °C of each other, indicating 
complete nucleation over the entire surface. At this point 
(and above), agreement with the previous data is excellent. 
The discrepancies in this mid heat flux region demonstrate the 
random behavior of nucleate boiling incipience; this is 
further investigated later in this chapter. After all sites 
nucleated (Point B) , the standard deviation was 0.9 °C with an 
average thermocouple temperature of 12.9°C. The point of 
incipience and its random occurrence at different levels of 
superheat will be examined later in this chapter. 

Sugiyama [Ref .2] did not conduct tests on a 19 fpi GEWA- 
K tube. This is one of the most common tubes used in 
commercial evaporators. Therefore, base-line data for the 19 
fpi GEWA-K tube with R-114 and R-114/oil mixtures were taken 
during this thesis before evacuating and cleaning the 
apparatus, and filling with R-124. The R-114 data obtained was 
incorporated in the R-124/R-114 comparison of the 19 fpi tube 
( Figure 5.12) . 

Figure 5.3 shows the repeatability of decreasing heat 
flux data for the newly constructed 19 fpi GEWA-K tube using 
R-124. The good agreement is an indication of the consistency 
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Figure 5.3 Repeatability Comparison of R-124 Data with 0% Oil 
for 19 fpi GEWA-K Tube (decreasing flux) 
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of the data reduction program DRPGB , the apparatus and the 
operator. Decreasing heat flux runs provide better comparison 
as temperature overshoot and hysteresis effects are avoided. 

B. POOL BOILING HEAT -TRANSFER COEFFICIENT CORRELATIONS 

Among the correlations available to predict pool boiling 
heat transfer coefficients for a smooth tube, most can be 
expressed as heat flux as a function of wall superheat. The 
pure R-124 smooth tube data is fitted very closely by the 
equation : 

q=17 . 0 aT 3 - 0 (1) 

A correlation in terms of fluid physical properties and 
tube geometry is given by Chongrungreong and Sauer [Ref. 7] . 
For R-124 at 177 kPa, this reduces to: 

q=45 . 8 aT 2 ' 32 (2) 

Another simpler correlation based only on system pressure 
is also given by Chongrungreong and Sauer [Ref. 7] : 

q=76 . 9 aT 2 ' 22 (3) 

Stephan and Abdelsalem [Ref. 6] present a statistical 
correlation for four different groups of fluids. For 
refrigerants, their correlation is presented in the form: 

q= (CIaT) 3 • 92 (4) 

where Cl is given as a function of refrigerant type and system 
pressure. Unfortunately, as R-124 is a relatively new 
refrigerant, their graphical representation of Cl does not 
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include R-124. However, Che present data could be fitted to an 
equation of the form given in equation (4) : 

q= ( 1 . 34 AT ) 3 • 9 2 (5) 

Hence, at a pressure of 1 atm, Cl= 1.34. 

Figure 5.4 shows the present R-124 smooth tube data for 
decreasing heat flux along with the above correlations. At 
lower heat fluxes the prediction of the Chongrungreong and 
Sauer [Ref. 7] correlation (equation (2) above) exhibits 
considerably closer prediction than the other correlations. In 
the more complex nucleate boiling region, the Stephan and 
Abdelsalem [Ref. 6] correlation and the Chongrungreong and 
Sauer [Ref. 7] correlation (equation (3)) display a fair 
prediction of the slope. 

C. SMOOTH TUBE DATA IN REFRIGERANT/OIL MIXTURES 

Although tests were conducted with oil concentrations of 
0, 1,2,3, 6 , Sc 10%, only data with 0,3 and 10% have been plotted 
in most of the figures for clarity. Figures 5.5 and 5.6 show 
the heat transfer performance of pure R-124, R-124/3% oil, and 
R- 124/10% oil for a smooth tube during increasing and 
decreasing heat flux respectively. It can be seen that at the 
highest heat fluxes, as oil concentration is increased to 3%, 
the heat transfer is enhanced by a maximum of 17% over pure R- 
124. For further increases in oil concentrations to 10%, 
performance drops off. This drop-off in thermal performance 
has also been observed by Jensen and Jackman [Ref. 13] who 
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Figure 5.4 Comparison of Smooth Tube Performance for 
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Figure 5.5 Performance Comparison for Boiling R-124/ 0%, 3%, 
& 10% Mixtures for Smooth Tube (increasing flux) 
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Figure 5.6 Performance Comparison for Boiling R-124/ 0%, 3%, 
& 10% Mixtures for Smooth Tube (decreasing flux) 
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attributed the decline to a steeper oil concentration 
gradient, inhibiting bubble growth. 

When oil is mixed with R-124, two competing mechanisms 
occur: (1) Oil accumulates in the vicinity of the vapor bubble 
surface as the more volatile refrigerant is evaporated; 
(2) diffusion of this accumulated oil back into the bulk 
mixture establishes an equilibrium concentration gradient 
near the departing bubble. Mechanisms (1) and (2) should act 
to inhibit bubble growth and lead to an overall decrease 
in the heat transfer with oil addition compared with pure 
refrigerant. It therefore appears that there is an additional 
mechanism that has a significant impact on the thermal 
performance of pool boiling of refrigerants with small 
quantities of oil. For R-124/oil mixtures of 3% and greater, 
significant foaming in the evaporator was observed at mid to 
high heat fluxes (> 10000 kW/m 2 ). It is possible that a third 
mechanism which could be affecting heat transfer is this 
foaming that occurred within the bulk R- 124/oil mixture. 
Udomboresuwan and Mesler [Ref .14] reported enhancement in pool 
boiling in the presence of foam. They attributed this 
enhancement to the foam causing the liquid/vapor interface to 
be closer to the tube wall. This would create a thin liquid 
film between a vapor bubble and the heated tube surface, 
thereby providing large heat - transfer coefficients. Higher 
heat fluxes and greater oil concentrations tend to increase 
the foaming, which could explain why this enhancement is seen 
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more predominately at high heat fluxes. However, as oil 
concentration increases, eventually the oil concentration 
gradient becomes so steep that mechanisms (1) and (2) tend to 
dominate and heat transfer performance drops off (ie. the 
foaming enhancement effect is offset) . 

At low heat fluxes, the scatter in the data is attributed 
to uncertainty. For the increasing heat flux data (Fig 5.5), 
it appears that the addition of oil tends to delay the point 
at which boiling commences. Due to the statistical nature of 
incipience, however, one should be wary of drawing much 
conclusion from this figure. The effect of oil on incipience 
is discussed in greater detail later in this chapter. 

Figure 5.7 compares smooth tube decreasing heat flux data 
from this thesis for R-124 with Sugiyama [Ref. 2] R-114 data 
at oil concentrations of 0%, 3%, and 10%. At all levels of 
heat flux, the R-124 data consistently demonstrate improvement 
in heat transfer over R-114. This ranges from 50% at low heat 
flux to 15% at high heat flux. The heat flux from a boiling 
surface is proportional to bubble diameter (d b ) , bubble 
departure frequency (f b ) , number of active nucleation sites 
per unit area, vapor density (p) , and latent heat of 
vaporization (h fg ) . Barthau [Ref. 15] concluded, after 
determination of d b , f b , and active nucleation site density 
for a smooth copper tube in R-114 using optical techniques, 
that as saturation pressure was increased, d b decreased, while 
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Plgar* 5.7 Comparison of HCFC-124 and HCFC-124/OH Mixtures 
to CFC-114 and CFC-114/Oil Mixtures for Smooth Tube 



f b , nucleation site density and heat transfer increased. While 
conducting the tests at 2.2 °C the corresponding saturation 
pressure of 11 psig for R-124 is greater than the pressure of 
0 psig for R-114. This suggests that the nucleation site 
density was greater for R-124 than for R-114 at the same 
saturation temperature. Added to this is the higher latent 
heat of vaporization, vapor density, and thermal conductivity 
for R-124 at this saturation temperature. The R-114/R-124 
property comparison at 2.2 °C saturation temperature is 
summarized in Table 1. It is therefore reasonable to expect a 
higher heat transfer coefficient for R-124 at the tested 
saturation temperature of 2.2 °C. 

D. BOILING PERFORMANCE OF 19 AND 26 FPI GEWA-K TUBES 

Figures 5.8 and 5.9 show the heat transfer performance of 
a 19 fpi GEWA-K tube in pure R-124, R- 124/3% oil, and R- 
124/10% oil during increasing and decreasing heat flux 
respectively. There is a consistent increase in thermal 
performance for oil concentrations of 3% and 10% over pure R- 
124 at high heat fluxes. This increase is substantially 
greater than that obtained with the GEWA-K 26 fpi tube for the 
same oil mixtures, as shown in Figures 5.10 and 5.11 for 
increasing and decreasing flux respectively. Surprisingly for 
both tubes, there seems to be a significant influence of oil 
in the natural convection region (Figs 5.8 & 5.10). There 
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Figure 5.8 Performance Comparison for Boiling R-124/ 0%,3%,& 
10% Mixtures for 19 fpi GEWA-K Tube (increasing flux) 
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Figure 5.9 Performance Comparison for Boiling R- 124 /0%, 3%, 10% 
Mixtures for 19 fpi GEWA-K Tube (decreasing flux) 
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Figure 5.10 Performance Comparison for Boiling R-124/ 0%,3%, 
& 10% Oil Mixtures for 26 fpi GEWA-K Tube (increasing flux) 
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Figure 5.11 Performance Comparison for Boiling R-124/ 0%,3%, 
& 10% Oil Mixtures for 26 fpi GEWA-K Tube (decreasing flux) 
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tends to be a decrease in performance for increasing oil 
concentration . This was not so obvious with Sugiyama [Ref. 2] 
R- 114 data. Once again, the point of incipience seems to be 
slightly delayed with increasing oil percentage. The thermal 
performance of the finned tube is better than the smooth tube, 
evidenced by thermal enhancements of 2.1 to 2.9 with the 19 
fpi tube and 1.3 to 2.1 with the 26 fpi tube. However, it 
should be noted that the heat flux calculation does not 
account for the increased surface area of the finned tube, but 
rather uses a root diameter. Some of the enhancement is 
therefore due to the increased surface area of the fins (Figs. 
5.24-5.26 compare heat transfer coefficients with smooth 
tube) , with each fin behaving essentially as the copper 
surface of a smooth tube in terms of thermal performance. 

Figures 5.12 and 5.13 compare R-124 and R-114 decreasing 
heat flux data for the 19 and 26 fpi GEWA-K tubes respectively 
for 0%, 3%, and 10% oil concentrations. The 19 fpi tube shows 
considerable improvement over the entire heat flux range while 
the 26 fpi tube shows improvement only in the low heat flux 
region. The improvement of the lower fin density tube can be 
speculated to be due in large part to better circulation and 
mixing between the fins. For both tubes at low oil 
concentrations and high heat flux, the foaming effect provides 
improvement in heat transfer similar to that seen with the 
smooth tube . 
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Figure 5.12 Comparison of IICI-C-124 and 1 1CI-C- 1 24/Oil Mixtures 
to CI ; C-114 and CTC-1 14/Oil Mixtures for 19 fpi Tube 
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Pigure 5.13 Comparison ofHCFC-124 and HCpC-124/Oil Mixtures 
to CFC-114 and CFC-114/Oil Mixtures for 26 fpi Tube 



E. BOILING PERFORMANCE OF HIGH FLUX AND TURBO -B TUBES 

Figures 5.14 and 5.15 show the heat transfer performance 
of a HIGH FLUX tube in pure R-124, R-124/3% oil, and R-124/10% 
oil during increasing and decreasing heat flux respectively. 
Uncertainty analysis (Appendix D) was performed on the HIGH 
FLUX tube and typical uncertainty bands are also given. The 
uncertainty is largest at low wall superheat and low heat flux 
since the uncertainty is of the order of the measured values. 
It can be seen that the heat transfer performance decreases 
for all heat fluxes as oil concentration increases. This 
decrease in performance is due most likely to the oil becoming 
trapped within the porous coating. As the refrigerant is 
evaporated, a very steep concentration gradient is developed 
within the pores; oil diffusion is also retarded because of 
the presence of the interconnected pores; this is particularly 
severe at high heat fluxes with 10% oil mixture, where a very 
large drop-off in performance is seen. In the natural 
convection region, the effect of oil is not as pronounced as 
at high flux, exhibiting only minor reduction in heat transfer 
compared to the pure refrigerant. Any advantageous effect of 
foaming is negated due to the pores remaining impenetrable to 
the foam. Figure 5.16 compares R-124 and R-114 decreasing heat 
flux data for the HIGH FLUX tube for 0%,3%, and 10% oil 
concentrations. The trend of the data for both fluids is very 
similar. With such a high density of nucleation sites already 
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Figure 5.14 Performance Comparison for Boiling R-124/0%,3% 
& 10% Oil Mixtures for HIGH FLUX Tube (increasing flux) 
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Figure 5.15 Performance Comparison for Boiling R-124/0%,3% 
& 10% Oil Mixtures for HIGH FLUX Tube (decreasing flux) 
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Figur* 5.16 Comparison of HCFC-124 and HCFC-124/Oil Mixtures 
to CFC-114 and CFC-114/OH Mixtures for HIGH FLUX Tube 



present on the surface, there does not appear to be a 
significant thermal enhancement at higher operating pressures 
as seen with the finned and smooth tubes. 

Figures 5.17 and 5.18 show the heat transfer performance 
of a TURBO-B tube in pure R-124, R-124/3% oil, and R-124/10% 
oil mixtures during increasing and decreasing heat flux 
respectively. The results of oil addition are similar to the 
HIGH FLUX tube. The notable exception is that at high oil 
concentrations and high heat fluxes, the heat transfer 

performance of the TURBO-B tube does not drop off as rapidly 
as the HIGH FLUX tube. Again, only a small effect of oil in 
the natural convection region was observed. 

Figure 5.19 compares R-124 and R-114 decreasing heat flux 
data for the TURBO-B tube for 0%, 3%, and 10% oil 

concentrations. Again, this comparison is similar to the 
results already discussed with the HIGH FLUX tube. 

F. EFFECT OF R- 124/OIL MIXTURES ON INCIPIENCE OF NUCLEATE 

BOILING 

Separate tests were conducted to analyze the incipience of 
nucleate boiling or the onset of nucleate boiling (ONB) in 
pure R-124 and how this was affected by the presence of oil. 
These tests were taken in accordance with the general 

procedures outlined in the preceeding chapter. However, 

smaller heat flux increments were applied to the test tube 
during the natural convection region in order to obtain the 
most accurate wall superheat at the ONB. Also, during these 
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Figure 5.17 Performance Comparison for Boiling R- 124/0%, 3%, 
& 10% Oil Mixtures for TURBO-B Tube (increasing flux) 
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rigur* 5.19 Comparison of HCFC-124 and HCFC-124/Oil Mixtures 
to CFC-114 and CFC-114/Oil Mixtures for TURBO-B Tube 



incipience tests the refrigerant pool was cooled down at 
approximately half the rate of the prior increasing and 
decreasing heat flux runs. Each incipience test was conducted 
with a minimum of 8 hours between between re- testing to allow 
adequate time for the R-124 pool to settle and for the 
nucleation sites to de-activate. The maximum time allowed for 
the aging of the tube surface was 48 hours. Similar tests were 
repeated seven times on the smooth and HIGH FLUX tubes. Figure 
5.20 summarizes the data on a probability chart. Due to the 
statistical nature of ONB, You et al [Ref. 16] proposed this 
chart as the best means of presenting incipient data. The 
seven incipient wall superheats for each condition are given 
for the smooth and HIGH FLUX tube at 0%, 3%, and 10% oil 
concentrations as the percent probability of nucleation 
occurring. For example, the smooth tube in pure refrigerant 
never nucleated at a wall superheat less than 12°C. Also the 
data show that for a smooth tube in pure R-124, nucleation 
occurred four times at a wall superheat greater than 20 °C, 
and three times at a lower wall superheat; this corresponds to 
a probability of nucleation occurring at a superheat of 20 °C 
between 3/7 (43)% and 4/7 (57)%. As the superheat further 
increases, the probability increases also until at a superheat 
of 28°C there is a 100% probability of nucleation. Bar-Cohen 
[Ref. 3] also reported widely varying R-113 superheats when 
'departure from convection' or nucleation occurred. The 
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probability chart of Figure 5.20 can be made more accurate by 
plotting a greater number of incipient superheats for each 
tube . 

For the smooth tube, there is a distinct shift to higher 
superheats with increases in oil concentration. This 
represents a delay in incipience caused by the addition of 
oil. Memory and Marto [Ref. 10] conducted similar tests with 
R- 114 and R- 114/oil mixtures and reported the oil 
concentration to have no systematic effect on the point of 
incipient boiling for the smooth tube. Wanniarachchi et al 
[Ref. 9] reported an increase in the wall superheat needed to 
initiate boiling for the smooth tube in R-114 with 10% oil, 
noting that 3% oil caused no increase. Memory and Marto [Ref. 
10] precicted the effect of oil on superheat needed to 
nucleate from a vapor- filled cavity by incorporating the 
physical properties of the refrigerant/oil mixture into a 
nucleation parameter (N) given by Marto and Lepere [Ref. 17] : 

' ^sat/ Pvapor ' ^fg 

If the parameter were to increase, a delay in the incipient 
point with increasing oil concentrations would appear likely. 
However, it is not clear how the oil effects this parameter; 
with the smooth tube in R- 124/oil mixtures it resulted in a 
delay in incipience, while the HIGH FLUX tube exhibited no 
distinct delay in incipience. Also the HIGH FLUX tube has a 
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whole range of sites that can be easily activated while the 
smooth tube only tends to have very small sites which may move 
outside the range of sites that can be activated when a 
ref rigerant/oil mixture is used. 

In the case of pure refrigerant, it does appear that with 
the relatively few incipient tests conducted, the wall 
superheats at incipience for the smooth and HIGH FLUX tubes 
are larger on average for R-114 than for R-124. This initial 
result would support the approximately 50% greater nucleation 
parameter of R-114 over R-124. This subject requires further 
study . 

Figure 5.21 shows an expanded view of the smooth tube's 
heater sleeve rolled out to show the exact locations of each 
thermocouple (longitudinal and circumferential) imbedded on 
the surface of the sleeve. Several data sets were taken for 
the smooth tube and then repeated after reversing the position 
of the tube within the pool. This procedure was conducted to 
see if there were any noticeable trends in the thermocouple 
readings which may be an indication of local hot or cold 
regions within the pool. However, this reversal of the tube 
had no consistent effect on the thermocouple readings or the 
point at which nucleation occurred. This appears to underline 
the uniformity of thermocouple output within the pool. A more 
detailed survey of this data needs to be done. 
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Figure 5.21 Expanded View of Thermocouple Locations within 
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G. SUMMARY OF THE R-124 DATA 

Figures 5.22 and 5.23 summarize Che effect of all oil 
concentrations (0, 1, 2 , 3 , 6, &10%) in R-124 for the five tubes at 
heat fluxes of 25 and 90 KW/m 2 respectively; 25 kW/m 2 
represents a typical operating heat flux. For each tube, the 
ratio of the heat - transfer coefficient at a given oil 
concentration divided by that for the same tube in pure R-124 
is given. At low heat flux, the smooth tube and the 19 fpi 
tube give an increase in performance as oil is added, up to 
concentrations of 6%. The smooth tube does the best in this 
respect, giving an increase in heat transfer of about 23% 
(compared to a smooth tube in pure R-124) at 6% oil 
concentration. The HIGH FLUX and TURBO -B tubes show a decrease 
in performance as oil is added, while the 26 fpi tube shows no 
change. At the highest oil concentration (10%), performance 
for all tubes is worse than that at 6%, although the smooth 
tube still indicates better heat transfer than when no oil was 
added, which is rather surprising. Similar trends were 
obtained for R-114, although at 10% oil, all tubes had poorer 
heat transfer than with no oil . At higher heat flux, general 
effects are similar, although now the 19 fpi tube shows the 
best performance of about 35% at 6% oil concentration. 
Interestingly, between 3 and 6% oil concentrations, all tubes 
except the HIGH FLUX tube show improved heat transfer over the 
same tube in pure refrigerant. Although below 3% oil there is 
an unexplained 'dip' in the heat transfer coefficient ratio 
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with the TURBO-B tube. The HIGH FLUX surface has very poor 
performance at 10% oil concentration, being only 35% of that 
in the pure fluid. 

Figures 5.24, 5.25 and 5.26 compare the five tubes 
directly for oil concentrations of 0, 3, and 10% respectively. 
From these figures, enhancements (ie. improvement in heat 
transfer compared to the smooth tube) can be calculated at any 
desired heat flux. These are provided in Table 4 at heat 
fluxes 25 and 90 kw/m 2 . The advantages of enhanced surfaces 
are clearly apparent. The smooth tube always has the poorest 
performance (as expected) and is used as a base-line with 
which the other tubes are compared. Between the finned tubes, 
the 19 fpi tube consistently outperforms the 26 fpi tube at 
all heat fluxes and oil concentrations: enhancements range 
between 2. 1-2. 9 for the 19 fpi tube and 1.3 -2.1 for the 26 fpi 
tube. Indeed, at the highest heat fluxes at 10% oil 
concentration, the 19 fpi tube is the best tube tested 
(enhancement of 2.6) . The GEWA-K 19 fpi tube to smooth tube 
area ratio is 2.4, while the GEWA-K 26 fpi tube to smooth tube 
area ratio is 3.2. The HIGH FLUX and TURBO-B tubes are the 
best tubes over most conditions of heat flux and oil 
concentration, with both having fairly equal performance. The 
exception is the HIGH FLUX tube at high heat flux and oil 
concentration, where performance drops off rapidly, as already 
mentioned. Over practical heat flux ranges (15-30 kW/m 2 ) 
therefore, the porous coated and modified fin tubes are the 
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best for all oil concentrations. Although the TURBO-B and HIGH 
FLUX tube performance is degraded at 10% oil concentration, 
in-service oil separators should maintain a 3% or less 
concentration in an actual evaporator system. At this level, 
the TURBO-B and HIGH FLUX tube are superior at normal chiller 
heat flux between 10 and 30 KW/m 2 . 

A point of interest worth noting was that the Viton and 
neoprene O-rings used experienced considerable swelling of 
approximately 30% by volume when subjected to R-124. This was 
enough to prevent their re-use when the Teflon inserts were 
removed to change the tube. Also noted was that the 0- rings 
only swelled when taken out of the system. When undisturbed, 
they maintained a good seal (this might not be true if the 
tube were sliding past the 0-ring) . Also significant, after 
removed from use within the apparatus for 60 days, both types 
of 0-rings did shrink back to their original size. 



Table 4. ENHANCEMENT SUMMARY FOR R-124 



Heat Flux 


25 kW/m 2 


90 kW/m 2 


Oil Concentration 


0% 


3% 


10% 


0% 


3% 


10% 


Smooth 


1 . 0 


1.0 


1. S 


1.0 


1 . 0 


1.0 


19 fpi 


2.9 


2.7 


2 . 8 


2 . 1 


2 . 4 


2 . 6 


26 fpi 


2 . 1 


1.7 


1.9 


1.3 


l.f 


1 . 5 


HIGH FLUX 


6.2 


4.9 


4 . 8 


2.8 


2 . 3 


1 . 1 


TURBO - B 


6 . 0 


4.3 


3 . 7 


2 . 4 


2.2 


1 . 9 
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Figure 5 . 22 Effect of Oil Concentration on I leal Transfer 
Coefficient at a Heat Flux of 25 kW/m* 
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Figure 5.23 Effect of Oil Concentration on Ileal Transfer 
Coefficient at a Heat Flux of 90 kW/m* 
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Heat Flu*/( W/m A 2) 



VI. CONCLUSIONS 



1. For the smooth and finned tubes, significant improvements 

in heat transfer (>50%) were obtained in pure R-124 

compared with pure R-114 at the same saturation 

temperature. This was attributed to the higher saturation 
pressure of R-124 activating more nucleation sites and the 
higher latent heat of vaporization, vapor density, and 
thermal conductivity of R-124 over R-114. 

2. For the HIGH FLUX and TURBO-B surfaces, performance in the 
two refrigerants was similar, since the number of active 
nucleation sites is already high. 

3. For pure R-124, the HIGH FLUX tube gave the largest 

enhancement of any tube over the whole range of heat flux. 

4. For the smooth and finned tubes, significant improvements 
in the heat transfer (>50%) were also obtained in R- 
124/oil mixtures when compared with R-114/oil mixtures at 
the same saturation temperature. This is also assumed to 
be due to the higher saturation pressures associated with 
R-124 . 

5. Increasing the oil concentration in R-124 by up to 6% for 
the smooth and finned tubes improved heat transfer 
performance compared to the pure refrigerant. For example 
at a heat flux of at 90 kw/rrr , performance for the smooth 
and finned tubes increased by 18% and 35% respectively. At 
higher oil concentrations, performance deteriorated. 

6. Any increase in oil concentration in R-124 for the HIGH 
FLUX or TURBO-B tubes led to a decrease in performance 
such that at high heat fluxes and high oil concentrations, 
the HIGH FLUX tube gave only 3 5% of its performance in 
pure R-124 and the 19 fpi tube gave the best overall 
performance . 

7. At chiller design heat fluxes (10-30 kW/m 2 ) and oil 
concentrations of 3% maximum, the HIGH FLUX and TURBO-B 
surfaces still provide the best heat transfer which is 
about the same as with R-114. Therefore, similar or 
slightly better performance is expected with R-124 when 
used as a 'drop-in' replacement for R-114 in a refrigerant 
evaporator operated at the same saturation temperature. 

8. Addition of oil delays the point of incipience (occurs at 
a higher wall superheat) for the smooth tube. For the HIGH 
FLUX tube, oil addition has no measurable effect on the 
wall superheat at incipience. 
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VII. RECOMMENDATIONS 



1. The remaining enhanced tubes should be tested in R-124 at 
0%,3%, and 10% oil concentrations and compared to the five 
tubes above and their prior performance in R-114 and R- 
114/oil mixtures. 

2. Incipience tests performed with finned tubes and TURBO-B 
tubes should be conducted to compare the influence of oil 
on wall superheat. 

3 . The influence of the variation in time between test runs 
on boiling nucleation (ie. onset of nucleate boiling) must 
be investigated more thoroughly. 

4. The single tube apparatus should be modified to test an 
array of two horizontal boiling tubes and compare these 
results with the R-114 tube array data previously 
recorded . 

5. A prototype chiller presently using R-114 should be tested 
with R-124 to monitor the performance. 
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APPENDIX A. THERMOPHYSICAL PROPERTIES OF R-124 



The following thermophysical properties of saturated R-124 
are plotted versus temperature in degrees Celsius. They were 
generated from REFPROP [Ref. 12] . 



1 . 


Saturation Pressure / (kPa) 




2 . 


Liquid Density / (kg/m 3 ) 




3 . 


Latent Heat of Vaporization /(kJ/kg) 




4 . 


Liquid Specific Heat /(kJ/kg K) 




5 . 


Liquid Viscosity /(mP) note: 1 mP= IxlO' 4 kg/m s 




6 . 


Liquid Thermal Conductivity / (W/m K) 
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Saturation Pressure for R-124 
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. 94 6 (T) +9 . 2081xl0' z (T) z +6. 93b9xlO 4 (T) 



Liquid Density for R-124 




(£iu/§>i)/oqH 



K> 



II 

CL 
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1434. 8-2. 8619(T)-6. 7267x10 3 (T) 2 -7.2852xlO' 5 (T) 




Latent Heat of Vaporization for R-124 







-C 
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= 159 .93-. 4609 (T) -1 . 4 58X10' 3 (T) 2 -l . 3715X10' 5 (T) 



Liquid Specific Heat at Constant P for R-124 




(M 35!/r^l)/d3 
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Cp-1.054 2+2. 14X10 3 (T) +1 . 0709X10' 5 (T) 2 -6 . 4 72 lxlO' 8 (T) 

~1.4324x10' 9 (T)*+4.136x10' 11 (T) 5 



Liquid Viscosity for R-124 







oo r- vo vn ^ cn cs 



(dlU)/DSIA 
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\i- 3. 4593-4. 26xlO‘ 2 (T)+3. 94 58xlO' 4 (T) 2 -4. 193x10 6 (T) 3 +2 . 07x10 8 (T) 

(note: lmP= lxlO ' 4 kg/m s) 



Liquid Thermal Conductivity for R-124 




(Mui/AV)/* 



o 

rH 

X 

rH 

<J\ 

rH 

\D 

• 

r* 



II 

X 
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. 5436x10’* (T) -1 . 9545xl0' 7 (T) 2 +3 . 18 3 5xl0' 9 (T) 



APPENDIX B. REPRESENTATIVE DATA SET 



Date : 19 Nov 1992 

NOTE: Program name : 0PPS8 
Disk number * 01 

Date : 19 Nov 1992 

NOTE: Program name : 0RPG8 
Disk number - 00 
New file name: 01119193 



TC 


are defective 


at locat ions 


7 ; 


3 






Tube Number: 13 












Data Set 


Number * 1 


Sulk Oil X - 


3.0 








TIME: 


15:40:01 














TC No: 


1 2 


3 


4 


5 


6 7 


8 




Temp : 


3.23 3.38 


3.31 


3.20 


3.06 


3.32 -99.99 


-99. 


99 


Twa 


TliQd TliQd2 


Tvapr 


Psat 


T sump 








3.25 


2.14 2.14 


2.18 


10.94 


-13.3 








Thetab 


Htube Qdp 












1 .06 3 


S . 475E+02 S. 


879E+02 












Data Set 


Number - 2 


Bulk Oil X - 


3.0 








TIME: 


15:42:22 














TC No: 


1 2 


3 


4 


5 


6 7 


8 




Temp : 


3.58 3.80 


3.72 


3. SI 


3.45 


3.64 -99.93 


-99. 


99 


Twa 


TliQd TliQd2 


T vapr 


Psat 


T sump 








3. S3 


2.22 2.24 


2.25 


10.91 


-13.2 








Thetab 


Htube Qdp 












1 .377 


6.422E+02 8. 


945E+02 












Data Set 


Number • 3 


Bulk Oil X • 


3.0 








TIME: 


15:45:20 














TC No.- 


1 2 


3 


4 


5 


6 7 


8 




Temp : 


3.85 4.35 


4 . 36 


4.19 


3.37 


3.95 -99.99 


-99. 


99 


Twa 


TliQd TliQd2 


T vapr 


Psat 


T sump 








4.09 


2.09 2.21 


2. 17 


10.83 


-13.1 








Thetab 


Htube Qdp 












1.921 


6 . 523E+02 1. 


253E + 03 












Oata Set 


Number * 4 


Bulk Oil X - 


3.0 








TIME: 


15:47: 12 














o 

z 

o 


1 2 


3 


4 


5 


6 7 


8 




T emp : 


4.10 5.03 


5.02 


4.87 


4 . 46 


4.28 -99.99 


-99. 


99 


Twa 


TliQd T 1 i qd2 


T vapr 


Psat 


T sump 








4.62 


2.15 2.27 


2.27 


10.92 


-13.0 









Thetab Htube Qdp 

2.350 6 . 986E+02 1.642E+03 
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Data Set 


Number * 


5 


Bulk Oil X - 


3.0 






TIME: 


15:49:09 














TC No: 


I 2 




3 


4 


5 6 


7 


8 


T e^p : 


4.30 S. 


75 


5.76 


5.44 


4.58 4.56 


-99.99 


-39.99 


Twa 


Tliqd T 1 i q d2 


Tvapr 


P sat 


T sump 






5.07 


2.05 2 


.20 


2.23 


10.33 


-12.9 






T het ab 


Htube 


Qdp 












2.345 


7 . 332E*02 


2 . 229E *03 










Data Set 


Number 3 


6 


Bulk Oil X = 


3.0 






TIME: 


15:51:30 














TC No: 


1 2 




3 


4 


5 6 


7 


8 


Temp : 


4.40 6.1 


52 


6.60 


6.21 


5.12 4.73 


-99.93 


-99.99 


Twa 


T I iqd T 1 iqd2 


T vapr 


Psat 


T sump 






5.61 


2.12 2 


. 16 


2.28 


10.93 


-12.3 






Thetab 


H t ube 


Qdp 













3.331 


9.I54E+02 3. 


0495*03 












Data Set 


Number * 7 


8u 1 k Oil X - 


3.0 








TIME: 


15:54:01 














TC No: 


1 2 


3 


4 


5 


6 


7 


8 


Temp : 


4.75 7.39 


7.27 


6.75 


5.39 


5.22 


-99.99 


-99.99 


Twa 


T 1 iqd T 1 iqd2 


Tvapr 


Psat 


T sump 








6.11 


2.15 2.21 


2.28 


10.93 


-12.7 








Thetab 


Htube Qdp 












3.331 


1.0315*03 3. 


95 1 E *03 












Data Set 


Number » 3 


8ulk Oil X * 


3.0 








TIME: 


15:56:43 














TC No: 


1 2 


3 


4 


5 


6 


7 


8 


Temp : 


5.14 8.2 Q 


3.15 


7.38 


5.55 


5.59 


-99.99 


-99.99 


T wa 


T 1 iqd T 1 iqd2 


Tvapr 


Psat 


T sump 








6.66 


2.24 2.27 


2.35 


11.00 


-12.6 








Thetab 


Htube Qdp 












4.302 


1.2S3E*03 5. 


4 1 3E *03 












Data Set 


Number * 9 


Sulk Oil X - 


3.0 








TIME: 


15:59:50 














TC No : 


1 2 


3 


4 


5 


6 


7 


8 


T emp : 


4.78 6.39 


8. 17 


6.85 


4.88 


4.46 


-99.99 


-99.99 


Twa 


Tliqd T 1 1 q d2 


T vapr 


Psat 


T sump 








5.97 


2.29 2.24 


2.32 


10.97 


-12.5 








Thetab 


Htube Qdp 












3.650 


2.1516*03 7. 


849E*03 












Oata Set 


Number » 1 0 


8ulk Oil X » 


3.0 








TIME: 


16:02:33 














TC No: 


1 2 


3 


4 


5 


6 


7 


3 


Temp : 


4.83 5.18 


5.13 


4.84 


4.91 


4.53 


-99.99 


-99.99 


Twa 


Tliqd T 1 1 q d2 


Tvapr 


Psat 


T sump 








CD 

cn 


2.30 2.25 


2.23 


10.93 


-12.4 








Thetab 


Htube Qdp 












2.533 


4.1505*03 1. 


072E*04 












Data Set 


Number » 1 1 


Sulk Oil X - 


3.0 








TIME: 


16:17:38 














TC No: 


1 2 


3 


4 


5 


6 


7 


8 


T emp : 


5.14 5.66 


5.59 


5.22 


5.25 


4.77 


-99.99 


-99.99 


Twa 


T 1 iqd T 1 iqd2 


T vapr 


Psat 


T sump 








5.20 


2.20 2.12 


2.16 


10.32 


-11.9 








Thetab 


Htube Qdp 












3.044 


5 . 0935*03 1. 


5525*04 
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Oata Set Nu»oer * 12 Bulk Oil X * 3.8 



TIME: 


16:20:31 












No: 


1 2 3 


4 


5 


6 


7 


8 


T e*p 


5 . "9 6.34 6.23 


5.90 


6.01 


5.24 


-99.99 


-99.99 


Twa 


7 l i q d T I iqd2 T vapr 


P a a t 


T Sump 








5.93 


2.26 2.19 2.20 


10.96 


-11.9 








Thetab 


Htube Qdp 












3.630 


5 . 328E*03 2 . 297E*04 












Oata Set 


Number * 13 8uik Oil X * 


3.0 








TIME: 


16:23:15 












TC No: 


1 2 3 


4 


5 


6 


7 


3 


Temp : 


6.61 7.18 7.13 


6.69 


6.85 


5.99 


-99.99 


-99.99 


T wa 


Tiiqd Tliqd2 Tvapr 


P sa t 


T sump 








6.60 


2.31 2.26 2.28 


10.93 


-11.8 








The t ab 


Htube Qdp 












4.316 


7.794E+03 3.364E+04 












Oata Set 


Number - 14 8ulk Oil X - 


3.0 








TIME: 


16:24:49 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


T emp : 


7.40 7.84 7.85 


7.36 


7.58 


6.49 


-99.99 


-99.99 


T wa 


Tiiqd Tliqd2 Tvapr 


P sa t 


T sump 








7.22 


2.35 2.27 2.32 


10.97 


-11.8 








T he tab 


Htube Qdp 












4.895 


9 . 726E + 03 4 . 76 1 E +04 












Data Set 


Number - is 8uU Oil X - 


3.0 








TIME: 


16:29: IS 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


Temp : 


8.17 8.32 8.52 


8.11 


8.32 


7.09 


-99.99 


-99.99 


Twa 


Tiiqd Tliqd2 Tvapr 


Psat 


T 3U«P 








7.78 


2.29 2.23 2.28 


10.93 


-11.7 








Thetab 


Htube Qdp 












5.501 


1 .3006*04 7. 1 51 E*04 












Oata Set 


Number » 16 8ulk Oil X - 


3.0 








TIME: 


16:33:55 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


T em p : 


8.79 8.96 9.33 


8.87 


9.12 


7.59 


-99.99 


-99.99 


Twa 


Tiiqd Tliqd2 Tvapr 


Psat 


T sump 








8.37 


2.20 2.14 2.20 


10.86 


-11.4 








Thetab 


Htube Qdp 












6. 169 


1 . 553E*04 9 . 577E *04 












Oata Set 


Number - 17 8ulk Oil X • 


3.0 








TIME: 


16:36:37 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


Temp : 


7.44 7.74 7.92 


7.59 


7.77 


6.56 


-99.99 


-99.99 


T wa 


Tiiqd Tliqd2 Tvapr 


Psat 


T sump 








7.24 


2.18 2.11 2.17 


10.83 


-11.4 








Thetab 


Htube Qdp 












5.070 


1 . 230E*04 6.2356*04 












Data Set 


Number - 18 8uU Oil X - 


3.0 








TIME: 


16:40:00 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


Temp : 


6.63 7.03 7.12 


6.77 


6. 98 


5.96 


-99.99 


-99.99 


Twa 


Tiiqd T 1 iqd2 T vapr 


Psat 


T sump 








6.57 


2.24 2.18 2.24 


10.90 


-11.4 








Thetab 


Htube Qdp 












4.329 


9 . 6576 + 03 4 . 1 90 E * 0 4 
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Data Set 


, Number « 19 8ulk Oil X ■ 


3.0 








TI"E: 


IS: 43:34 












7C Nc : 


i 2 3 


i 


5 


5 


7 


3 


Temp : 


6.09 6.65 6.55 


5 . 2 ^ 


6 . 4 i 


5.54 


-99.99 


-99.99 


Twd 


T 1 1 qd T 1 iqd2 T vapr 


Psat 


T s^mp 








S. 1 4 


2.27 2.20 2.24 


10.90 


-11.5 








T het aD 


Htube Qdp 












3.699 


7.859E+03 3 . 06“*E +04 












Data Set 


Number * 20 BuH Oil X = 


3.0 








TIME: 


1 6 : 4S : 5 9 












TC No : 


1 2 3 


4 


5 


6 


7 


8 


T emp : 


5.62 6.13 6.13 


5.32 


5.92 


5.19 


-99.99 


-99.99 


Twa 


T 1 iqd T 1 i q d 2 T vapr 


Psat 


T sump 








5.70 


2.23 2.14 2.19 


10.35 


-11.5 








The t ab 


Htube Qdp 












3.50G 


6.754E + 03 2. 363E *0 4 












Data Set 


Number - 21 Bulk. Oil X = 


3.0 








TIME: 


16:43:34 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


T emp : 


5.17 5.73 5.71 


5.36 


5.40 


4.85 


-99.99 


-99.99 


Twa 


Tliqd Tliqd2 Tvapr 


Psat 


T sump 








5.30 


2.31 2.22 2.25 


10.91 


-11.7 








Thet ab 


Htube Odp 












3.045 


5 . 555E+03 1 . 692E+04 












Data Set 


Number * 22 Bulk Oil X - 


3.0 








TIME: 


16:51 : 59 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


Temp : 


4.69 5.20 5.15 


4.35 


4.92 


4 . 44 


-99.99 


-99.99 


T wa 


Tliqd Tliqd2 Tvapr 


Psat 


T sump 








4.32 


2.27 2.17 2.21 


10.87 


-11.8 








Thetab 


Htube Qdp 












2.613 


4 . 525E + 03 1 . 1 82E +04 












Data Set 


Number ■ 23 Bulk Oil X ■ 


3.0 








TIME: 


16:54:33 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


T emp : 


4.32 4.77 4.63 


4.45 


4.55 


4.21 


-99.99 


-99.99 


T-a 


TliQd Tliqd2 Tvapr 


Psat 


T sump 








4 . 45 


2.30 2.21 2.27 


10.92 


-11.9 








Thetab 


Htube Qdp 












2.136 


3 . 653E + 03 7 . 985E+03 












Data Set 


Number - 24 Sulk Cl 


1 X - 


3.0 








TIME: 


1 6 : 57 : 5S 












TC No: 


1 2 3 


4 


5 


6 


7 


3 


Temp 


3.94 4.31 4.23 


4.10 


4.11 


3.91 


-99.99 


-99.99 


Twa 


Tliqd T l iqd2 T vapr 


Psat 


T sumo 








4.08 


2.24 2.15 2.21 


10.37 


-12.0 








Thetab 


Htube Qdp 












1 .363 


2 . 839E+03 5 . 289E+03 












Data Set 


Number - 25 Bulk Oil X - 


3.0 








TIME: 


17:01:00 












TC No: 


1 2 3 


4 


5 


6 


7 


8 


T emp : 


3.56 3.91 3.34 


3.75 


3.74 


3.60 


-99.99 


-99.39 


Twa 


Tliqd T 1 iq d2 Tvapr 


Psat 


T sump 








3.72 


2.22 2.03 2.14 


10.30 


-12.1 








Thetab 


Htube Qdp 












1 .577 


2. 1 66E + 03 3. 41 9E+03 
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Oata Set 


; Numbe* - «.o Sw i k O 4 I * * 


3 0 






TIf«: 


17:03:15 










TC *o: 


1 2 3 


4 


S 


6 


7 8 


T«mo : 


3.56 3.83 3.77 


3.74 


3.69 


3.61 


-99.99 -39.99 


T-e 


Tliqd TUqdC 


0 sat 


T S^mp 






3. S3 


2.35 2.21 2.27 


10.92 


* i 2 . 2 






Thetab 


Htuba Gdp 










1 .420 


1 . 809E*03 2.S -0E*03 










Gala Set 


Number - 27 8uU Oil X - 


3.0 






T I HE : 


17:06:05 










TC No: 


1 2 3 


4 


S 


6 


7 8 


T emp : 


3.41 3.63 3.56 


3.52 


3.48 


3.4S 


-99.99 -99.99 


Twa 


T 1 iqd T1 iqd2 Tvapr 


Psat 


T sump 






3.S0 


2.35 2.20 2.24 


10.89 


-12.3 






Thetab 


Htuba Qdp 










1.280 


I .S7IE*03 1 .98IE*#3 










Oata Set 


Number • 28 8uU Oil X • 


3.0 






T I HE : 


17:08:04 










TC No: 


1 2 3 


4 


5 


6 


7 8 


Tamp : 


3.24 3.42 3.34 


3.34 


3.30 


3.32 


-99.99 -99.99 


Twa 


Tliqd Tliqd2 Tvapr 


Psat 


T sump 






3.32 


2.26 2.15 2.14 


10.80 


-12.4 






Thetab 


Htuba Qdp 










l.f 8S 


1 .1S9E+03 1 .373E+03 










Oata Sat 


Number • 29 8ulk Oil X • 


3.0 






TIHE: 


17: 1 0 : $S 










TC No: 


1 2 3 


4 


5 


6 


7 8 


Temp : 


3.15 3.23 3.18 


3.16 


3.12 


3.20 


-99.99 -99.99 


Twa 


Tliqd Tliqd2 Tvapr 


Psat 


T sump 






3.17 


2.28 2.25 2.23 


10.89 


-12. S 






The tab 


Htuba Qdp 










.937 


1 . 0S6E + 03 9. 893E+02 










Oata Sat 


Number • 30 BuU Oil X • 


3.0 






TIH€: 


17:13:34 










TC No:- 


1 2 3 


4 


5 


6 


7 8 


Tamp : 


3.20 3.40 3.36 


3.33 


3.29 


3.28 


-99.39 -99.99 


Twa 


T 1 1 q d T 1 iqd2 Tvaor 


Psat 


T Sump 






3.30 


2.21 2.08 2.10 


10.77 


-12.6 






T ha tab 


Htuba Qdo 










i .20S 


1 . 246E m03 1 .S0iE*#3 










Oata Sat 


Number • 31 BuU Oil X • 


3.0 






TIHE: 


17:16:00 










TC Ho: 


1 2 3 


4 


S 


6 


7 8 


Temp : 


3. IS 3.2S 3.20 


3.16 


3.13 


3.23 


-99.99 -99.99 


Twa 


Tliqd Tl;qd2 Tvapr 


Psat 


T Sump 






3.18 


2.30 2.27 2.28 


10.93 


-12.7 






Thetab 


Htuba Qdp 










.907 


9.922E*02 9.001 E*02 










Oata Set 


Number • 32 8uU Oil X • 


3.0 






TIHE: 


17:17:39 










TC Ho: 


1 2 3 


4 


5 


6 


7 8 


Temp : 


2.92 3.00 2.92 


2.89 


2.9i 


3.03 


-99.98 -99.99 


Twa 


Tliqd T 1 iqd2 Tvapr 


Psat 


T Sump 






2.94 


2.15 2.17 2.15 


10.81 


-12.7 






Thetab 


Htube Qdp 










.796 


8 . 888E *02 7.074E*02 











Oats Set Number - 33 8ulk Oil X • 3.0 



TIHE: 


1 7: l 9 


: 16 




TC No: 


1 


2 3 4 


5 6 7 8 


Temp ' 


2.92 


2.95 2.98 2.86 


2.87 2.99 -39.59 -99.99 


Twa 


Tliqd 


Tliqd2 Tvapr Psat 


T sump 


2.91 


2.17 


2.19 2.16 10.82 


-12.8 


Thetab 


Htube 


Odp 




.753 


7.74SE 


*02 5 . 829E*02 
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APPENDIX C. SAMPLE CALCULATIONS 



Data file number D1221SM0 was used Co conduct the 
calculations. The saturation temperature was 2.28 °C 
single smooth tube heat flux of 31290 W/m 2 . 

A. TEST-TUBE DIMENSIONS 

D q =. 015875 m 
Dj_= .0127 m 
D x =. 01245 m 
L= . 2032 m 
Lu=.0762 m 



B. MEASURED PARAMETERS 

V=141.95 volts 
1=2.34 amps 
T1=16 . 84 °C 
T2=17 . 15 °C 
T3=16 . 42 °C 
T4=16 . 22 °C 
T5=15 . 66 °C 
T6=15 . 48 °C 

T sat = 2 - 28 ° C 
k c =344 W/m-K 



sample 
with a 
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C. OUTER WALL TEMPERATURE OF THE BOILING TUBE 



p=TT • D 0 =7r •. 015875 m =.04987 m 

Ac = 7r(D 0 2 -D 1 2 ) /4 = 7r( ( .015875) 2 - ( .0127) 2 ) /4 
Ac=7 . 13xl0' 5 m 2 

Q h =VI=141 . 95 • 2 . 34= 332.16 W 

T avg =ETn/n=E 6 Tn/6 = 16 . 3 °C 

T wo= T avg-t Qh‘ (ln(D 0 /D 1 )/(27r-L-k c ) )i 

T wo =16 . 3 - [ 332.16- ( In ( . 015 875 / . 0 124 5 ) / ( 2 tt • L • k c ) ) j 

T wo =16 • 12 °C 

(9=T wo -T saC = l6 . 12-2,28-13 . 84 °C 

D. PROPERTIES OF R-124 AT FILM TEMPERATURE 

The thermophysical properties of R-124 were obtained from 
REFPROP [Ref. 12] and are shown in Appendix A. 

T f = ( T wo +T sat ) /2= (16. 12+2 .28) /2= 9.2 °C 

^ = 3 . 4593 - ( . 0426)T f + (3 . 94 85x10 ' 4 ) T f 2 - (4 . 19 3x10 ' 6 ) T f 3 + 

(2 . 0709x10'®) T f 4 
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^=3 .459 3- (.042 6) ■ 9 . 2 + ( 3 . 9485x10 ' 4 ) • (9.2) 2 - (4.193xl0' 6 ) • (9.2) 3 + 
( 2 . 0709x10 ’ 8 ) • ( 9 . 2 ) 4 
p=3 . 09 8x10 " 4 N-s/m 2 

p = 1434 . 8 - ( 2 . 8 6 19 ) T f - ( 6 . 7267xl0' 3 ) T f 2 - ( 7 . 2 852x10 ' 5 ) T f 3 
p = 1434 . 8 - (2.8619) -9.2- (6. 7267xl0' 3 ) • (9 .2) 2 - ( 7 . 2 852x10 " 5 ) • (9 . 2) 3 
p=1407 . 84 kg/m 3 



v=(i/ p = 3 . 09 8x10 ' 4 / 1407.84= 2.201xl0' 7 m 2 /s 

k=7.5191xl0' 2 - (3 .5436x1 O’ 4 ) T f - ( 1 . 9 545x10 ‘ 7 ) T f 2 + ( 3 . 18 3 5x10 ' 9 ) T f 3 
k=7 . 5191X10' 2 - ( 3 . 5436x10 * 4 ) • (9.2) - ( 1 . 9 545x10 ' 7 ) • (9.2) 2 + 

(3 . 183 5x10 ’ 9 ) • (9 .2) 3 
k=7 . 192xl0' 2 W/m-K 

c p =l .0 542+ (2 . 1405x10 ‘ 3 )T f + (1 . 0709x10 ' 5 ) T f 2 - (6.4721xl0' 8 )T f 3 - 
(1 . 4324x10 ' 9 ) T f 4 + (4 . 13 6x10 ' 11 ) T f 5 

c p =l . 0 542+ (2 . 14 05x10 ' 3 ) ■ (9.2) + (1.0709xl0' 5 ) • (9. 2) 2 - (6.4 721 
xlO* 8 ) • ( 9 . 2 ) 3 - ( 1.4324x10' 9 ) • (9 . 2 ) 4 + (4 . 136X10' 11 ) • (9. 2) 5 

c p =1074 . 74 J/kg ■ K 

«=k/ (p • c p ) = 7 . 192x1 O' 2 / (140 7.84 • 1074 . 74) 
a=4.753xl0' 8 m 2 /s 
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/3= - ( Ap/AT) /p = 2 . 09 4x10 ' 3 (1/K) 



Pr= v /a=A . 63 



E. HEAT- FLUX CALCULATION 

The average natural convection heat - transfer coefficient 
at the non-boiling ends of the test tube is calculated using 
the Churchill and Chu [Ref. 4] correlation: 

For simplification let 1 represent the following in the 
correlation : 

R= [ (g • /3 • D Q 3 • 9 • tanh (m- Lu) ) / ( u • a • Lu • m) ] 1/<6 
then 

h = k/D 0 [ . 6 + . 3 87 [ 1 / (1+ (0 . 559/Pr) 9/16)8/27] ] 2 
where 

m=[(h-p)/(k c -Ac)] 1/2 

and an iterative process is used in program DRPGB to compute 
the value of h by starting with the value for h of 190 
(W/m 2 'K). The resulting values for h and m were: 

h=265 . 43 (W/m 2 • K) 

m=23 .23 (1/m) 
therefore : 

Q f = (h ■ p ■ k c • Ac) - 1 / 2 • 6 ■ tanh (m •Lu)=7.44 W 
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F. HEAT FLUX THROUGH ACTIVE BOILING LENGTH 



Q=Q h - 2'Q f =332.16- (2-7. 44 ) =3 17. 28 W 



Ab=7T'D 0 - L = tt - 0 .01587 5 • 0 . 2032 = 1 . 0134x10 ' 2 m 2 
q=Q/Ab=317 . 28/1 . 0134x10 ' 2 = 3 1 , 3 0 8 W/m 2 
h=q/0=31 , 308/13 . 84=2262 . 1 W/m 2 • K 

The following results were produced by the data 
acquisition and reduction program DRPGB: 
q=31290 W/m 2 
0=13.84 °C 
h=2262 W/m 2 • °K 



TABLE 5. DIMENSIONS OF BOILING TUBES TESTED 



TUBE 


Dl mm 


D2 mm 


Di mm 


Do mm 


L mm 


Lu mm 


K c 

W/moK 


smooth 


12 .4 


15 . 9 


12 . 7 


15.9 


203.2 


76.2 


344 


GEWA-K 
19 fpi 


10 . 1 


12 . 7 


li.l 


12 . 7 


203.2 


76 . 2 


344 


GEWA-K 
26 fpi 


10 . 1 


12.7 


19 . 1 


12 . 7 


203.2 


76.2 


344 


TURBO -B 


11 . 6 


13 . a 


11 . 8 


15 . 8 


203 .2 


76.2 


398 


HIGH 

FLUX 


12.9 


15 . 8 


13.2 


15 . 8 


203.2 


76.2 


45 * 



* HIGH FLUX tube porous coating is on a copper/nickel base 
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APPENDIX D. UNCERTAINTY ANALYSIS 



Using methods of uncertainty analysis outlined by Kline 
and McClintlock [Ref. 18] and performed by Sugiyama [Ref. 2], 
four data points are analyzed. The data points chosen were for 
a Smooth tube and HIGH FLUX tube. For each tube, there was a 
data point at a high heat flux and one at a low heat flux. The 
following is a sample uncertainty analysis calculation from 
data file number D1220SM0 with a Smooth tube, 0% oil, at a 
heat flux of 86,000 W/m 2 . The results of the remaining sample 
points are listed in Table 6. Uncertainties are given as a 
percentage of the calculated parameter, i.e. a relative 
uncertainty . 

A. UNCERTAINTY IN HEAT - TRANSFER RATE 

I s =1.99 volts V s =9.325 volts 

Estimated uncertainties in the measured quantities I and V : 

6 1 = + 0.025 amps 
<5V= + 0.05 volts 



1=1 . 9182 • I s =3 . 82 amps 

V=25 • V s =233 . 13 volts Q h =VI 

where: 6= uncertainty in measurement 
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<5Q h /Q h = ( (6V/V s ) 2 + (6I/I s ) 2 ) 1/2 
6Q h /Q h = ( (.05/9. 325) 2 +(. 02 5/1.99 ) 2 ) 1/2 
6Q h /Q h = 1.37 % 

B. UNCERTAINTY IN SURFACE AREA 



Ab= 7T • Do • 


L 






D 0 = 15.88 


(mm) 


60 o = 


. 1 (mm) 


L= 203 . 2 


(mm) 


<5L= 


. 1 (mm) 



5Ab/Ab= ((5D o /0 o ) 2 + (6L/L) 2 ) 1/2 
6Ab/Ab= ( ( . 1/15 . 88 ) 2 + (. 1/203 .2) 2 ) 1/2 
5Ab/Ab= .63 % 

C. UNCERTAINTY IN WALL SUPERHEAT 

^T= T WQ -T sat T sat =2.23 °C 5T sat = .01 °C 

T wo= T avg ' Q h [ In ( ) / ( 2 • 7T • L • k c ) ] 

Tn= thermocouple readings 

T avg= STn/n 

Tl=21 . 63 °C 
T2=21 . 41 °C 
T3 =2 0 . 6 3 °C 
T4 = 2 0 . 94 °C 
T5=20 . 06 °C 
T6=20 . 40 °C 
T avg =20.85 °C 
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standard deviation= ( (E (Tn-T avg ) 2 ) /n) 1//2 
standard deviation= 0.55 °C 

The logarithmic term in the equation for Two can be 
neglected for the uncertainty analysis as it is very small in 
comparison to the standard deviation. 

T wo= T avg=2 0.8 5 °C 

<5T wo =standard deviation= . 55 °C 
AT=18 . 62 °C 

5 AT/ AT = ( (6T wo /AT) 2 + ( 6T sat /AT) 2 ) 1/2 
6AT/AT= ( (. 55/18. 62) 2 +(. 01/18. 62) 2 ) 1//2 
<5AT/AT=2 . 95 % 

D. UNCERTAINTY IN HEAT FLUX 

q= (Q h -2'Q f )/Ab 

Q h =VI Q h =890 . 56 W 6Q h =12.2 W 

assuming the same proportion in the uncertainty for Q f : 

q=86000 W/m 2 Ab= 0.010137 m 2 

Q h - 2 • Q f = 871.81 W 

Q f = 9.38 W 6Q f =0 . 12 8 W 

6q/q= [ (5Q h / (Q h -2Q f ) ) 2 + (2 5Q f / (Q h -2Q f ) ) 2 + (6Ab/Ab) 2 ] 1/2 
6q/q= [(12. 2/(890. 56-2-9.38)) 2 +(2-(. 12 8)/( 890. 56-2'9.38)) 2 + 

(6 . 39x10 ' 5 / .010137) 2 ] 1/2 
5q/q= 1.54 % 
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E. UNCERTAINTY IN BOILING HEAT -TRANSFER COEFFICIENT 

h=q/AT 

5h/h= [ (<5q/q) 2 + (6AT/AT) 2 ] 1/2 

6h/h= [(. 0154) 2 + (. 0295) 2 ] 1/2 

6h/h= 3.3 % 



TABLE 6. UNCERTAINTY ANALYSIS OF FOUR DATA POINTS 



1 Parameters 


Smooth tube 


Smooth tube 


High 


Flux 


High Flux 


i 

\ 


4019 W/m 2 


86000 W/m 2 


2244 


W/m 2 


88860 W/m 2 


6Q h /Q h (%) 


3 . 52 


1.37 


8 . 11 


1 .36 


5Ab/Ab (%) 


. 63 


.63 


.634 


. 634 


3AT/AT ( % ) 


5.20 


2.95 


17.4 


12.5 


6q/q (%) 


1.29 


17.4 


8 .41 


1 . 5 


6h/h (%) 


5 . 4 


3 . 3 


19 . 3 


12 . 6 



All uncertainty analysis calculations were conducted 
during decreasing heat flux runs. The wall superheat 
uncertainty (6AT/AT) for the HIGH FLUX tube was significantly 
greater than that of the smooth tube due to a greater 
deviation in the thermocouple temperature readings. 
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APPENDIX E. SETUP PROGRAM 



The following setup program was used in the preparation of 
the apparatus prior to commencing testing. The program was 
written in Hewlett-Packard Basic 5.0 for both the Hewlett- 
Packard 9300 and 9852A series data acquisition/control unit. 
The setup program performs the following: 

1. Monitor the Sump temperature 

2. Monitor the evaporator liquid temperature 

3. Measurement and readout of all thermocouple channels. 

4. Measurement and readout of the power supplied to the 

tube cartridge heater. 

5. Measurement and readout of the power supplied to the 
auxiliary heaters if used. 



104 



1 » 
2 ' 

4 I 

5 » 
10 
20 

1 1 . 
21 
22 
23 
25 

27 

28 

30 

31 

32 

33 

34 

35 

36 

37 

38 

40 

41 

42 

43 

44 

45 
45 

47 

48 

491 

50 

51 

53 

54 
60 
70 
80 
50 
100 
1 1 0 
1 20 
130 
1 40 
141 
1 42 
150 



PROGRAM : SETUP 
DATE: AU6UST 3,1991 

PROGRAMMER: LT DEAN SUGIYAMA 
MODIFIED BY LANNIE LAKE JAN 22, 1992 
COM /Cc/ C ( 7 > 



DATA 0.10086091 ,25727.94369,-767345.9295,78025595.81 
. 66 1 92E + 1 3 
DATA 3. 94078EH4 
READ C< • > 

ON KEY 1,15 60T0 
PRINTER IS 1 
PRINT 
PRINT 

PRINT USIN6 ' 4X ,* 

PRINT USIN6 “6X , ■ 

PRINT USING * 6X , * 

PRINT USIN6 “ 6X .* 

PRINT USIN6 *6X 
PRINT USIN6 “6X ,* 

PRINT USIN6 ' 6 X ,* 

PRINT USIN6 ‘ 4X 
BEEP 

INPUT Ido 
IF Ido^5 THEN Ido 
IF I do=0 THEN 50 
IF Ido-1 THEN 155 
IF Ido-2 THEN 173 



-9247486589 ,5.975336^ 
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‘SELECT OPTION-- 
•0-MONITOR SUMP'*- 
- 1 -MONITOR LIQUID* “ * 

* 2 -CHECK THERMOCOUPLES*" * 

* 3 -CHECK MAIN HEATER* “ 
-4-CHECK AUX HEA T ERS* ** 
‘5-EXIT PROGRAM" * * 

•NOTE: KEY 1 - ESCAPE* * “ 



IF Ido-3 THEN 195 
IF Ido-4 THEN 135 
IF Ido-5 THEN 231 
PRINT 



PRINT 

PRINT ‘SUMP TEMPERATURE DE5 C * 
PRINT 

OUTPUT 709 i * AR AF 11 ALII UR5 * 
OUTPUT 709 < "AS SA* 

Su«-0 

FOR J-1 TO 5 
ENTER 709 i E 
Suw-Suw+E 
NEXT J 
Edve-Sum/S 
Tenp*FNT vsv( Eave ) 

PRINT USING "4X .MDD.OD* ;Te^p 

BEEP 

PRINT 

ulAIT 5 
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GOTO 50 



• 5 1 

i s: 

! 55 PRINT 

156 PRINT “LIQUID TEMPERATURE DEG C“ 

! 53 PRINT 

159 OUTPUT 709 5 * AP AF08 AL09 LR5' 

1 60 Sum-0 

161 FOR I-l TO 2 

162 OUTPUT 70S i 'AS SA” 

163 ENTER 709 ; E 

1 64 Sum=Sum+E 

155 NE<T I 

166 Eave*Sum/2 

167 Tewp«FNT w-sv( Eave ) 

168 PRINT USING “4X .MCO.OO'sTenp 

169 BEEP 

170 WAIT 5 

171 GOTO 155 
172' 

173 PRINT 

176 PRINT “CHANNEL TEMPERATURE DEG C* 

177 OUTPUT 709 i * AR AF00 ALU UR5 * 

178 FOR 1-1 TO 12 

179 OUTPUT 709; “AS SA “ 

180 Sum-0 

181 FOR J-1 TO 5 

182 ENTER 7 0 9 ; E 

1 83 Sum«Sum + E 

184 NEXT J 

585 Ed ve-Sum/5 

186 Temp »F nTv* v ( Ed v© > 

187 PRINT TA6(3);I;TAe( l5);Temp 

188 NEXT I 

189 BEEP 

190 WAIT 5 

191 GOTO 173 
1941 

195 PRINT 

196 OUTPUT 709 ; ‘ AR AF20 Al22 VR5“ 

197 FOR 1-1 TO 3 

198 OUTPUT 709 i '* AS SA* 

1 99 S u m = 0 

200 FOR J-1 TO 5 

201 ENTER 709;E 

202 Sum -Sum *E 
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233 

234 

205 

209 

210 

21 I 

212 

2 ! S 

217 

:i a 

219 

220 

221 

222 

223 

224 

225 

226 

227 

223 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

24’ 

242 

250 



NEXT J 

IF 1*1 Tr*EN Vo 1 1 * Sun/ 5 
IF 1*2 AND I do*3 THEN 

PRINT -make SURE V0LTA6E BOX IS SET TO MAIN HEATERS' 

Amp*Sum/5 

ENO IF 

IF 1*3 AND Ido*4 ThEN 

PRINT "HAKE SURE VOLTAGE BOX IS SE T TO Aux HEATERS' 

A mp *Su«/5 
ENO IF 
NE < T I 

Amp*ABS( Amp *1.9182) 

Volt =A63( Vo 1 t *25 ) 

Po*er*Vo 1 1 * Amp 
Resistence*Volt /Amp 
PRINT 
BEEP 

PRINT " VOLTAGE ( V ) CURRENT ( A ) RES I STENCE ( ohms > POUER<U)‘ 
PRINT 

PRINT USING ' 1 X ,5( MDDOD.OD ,4X >' i V o 1 1 ,Amp .Res i s tence .Power 

WAIT 5 

60T0 195 

BEEP 

PRINT 

PRINT ’THAT'S ALL FOLKS' * 

ENO 

DEF F NT v s v ( V ) 

COM /Cc/ C ( 7 > 

T«C < 0 ) 

FOR 1*1 TO 7 
T*T+C( I >*V‘ I 
NEXT I 

T- T+8.526397E-2* T *< 3 . ?S 1 I 99E-Z-T*5 . 0639255E-5 t 

RETURN T 

FNENO 
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APPENDIX F. PROGRAM DRPGB 



The data acquisition and reduction program DRPGB is 
written in Hewlett-Packard Basic 5.0 for the HP 9300 series 
computer and listed on the following pages. 
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10 

20 

30 

4 1 
50 
50 
70 
80 
90 
100 
1 10 
120 
1 30 
1 40 

150 

151 
160 
170 
180 

1 90 
Et( 1 
200 
210 
220 
230 
240 
250 
250 
270 
280 

2 90 
300 ' 
310 
320 ' 
330 
340 

350 ' 

351 
360 
370 
380 
390 
400 
4 1 0 
420 

430 

431 

432 

433 
440 
450 
460 
470 
480 
4 90 
500 
510 
520 
530 
540 i 
550 i 



i FILE NAME: 0RP68 
OATE: October 19, 1984 

REVISED: July 11, 1991 BY Oaan Sugiyama 

i REVISED: Aug IB, 1992 BY 6aorge Bertsch (including R - 1 2 4 properties) 

COM /I dp/ I dp 
PRINTER IS 1 
CALL Select 

INPUT "WANT TO SELECT ANOTHER OPTION < 1 -Y ,0-N ) ^ ’ *ael 

IF Isel-1 THEN GOTO 70 

BEEP 

BEEP 

PRINTER IS 1 

PRINT " OATA COLLECTION/REPROCESSING COMPLETEO “ 

ENO 

SUB Main 

COM /Iprop/ Ift 

COM /I dp/ I dp 

COM /Cc/ C( 7 > , lea 1 

COM /Wil/ 02 ,0 1 ,Do ,L ,Lu ,Kcu 

OIM Emf ( 12 > ,T( 12 > ,01 a( 1 3 > ,02a( 1 3 > ,Oia< 1 3 ) ,0oa< 13 ) ,La( 13 > ,Lua( 1 3 ) ,Kcua( 13), 
3) ,Tn$<4 >C 15] 

OATA 0.100B6091 ,25727.94359,-767345.8295,78025595.81 
OATA -9247486589 , 6 . 976BBE+1 1 ,-2.66!92E+13,3.9407BE+!4 
REAO C( • ) 

i OATA "Smooth’ , 'High F lux ’ , ’ Thermoaxe 1 -E " , “ Thermoex e 1 -HE ’ 

OATA Smooth, High Flux ,Turbo-B .High Flux Mod,Turpo-B Mod 
REAO TnS( • > 

PRINTER IS 701 
BEEP 

IF Idp-4 THEN 'PRINTER IS 1 
IF Idp-4 THEN GOTO 2660 

INPUT ’ENTER MONTH, OATE AND TIME ( MM : 00 : HH : MM : SS ) " ,Oat e$ 

OUTPUT 709 i ’TO’ ;OataS 
OUTPUT 709i ’TD’ 

ENTER 709 1 0at a$ 

PRINT 

PRINT M Month, date and time :";Oate$ 

PRINT ’ Oat a : ' ,OATE*( TIMEOATE ) 

PRINT 

PRINT USIN6 ' 1 0X , ’ "NOTE : Program name : ORPGB ’ " " 

BEEP 

INPUT 'ENTER DISK NUMBER’, On 

PRINT USING *16X,”0isk number - ”,22'iOn 

BEEP 

INPUT ’ENTER INPUT MOOE ( 0-3054A , 1 -FILE )’ ,1m 

BEEP 

If t-0 

INPUT ’SELECT FLUID ( 0- 1 1 4 , 1 - 1 24 ) ’ f I f t 
BEEP 

INPUT ’SELECT HEATIN6 MOOE (0-ELEC; 1 -WATER ) ’ , I hm 
BEEP 

INPUT ’ENTER THERMOCOUPLE TYPE < 0-NEW , 1 -OLG ) " , leal 

IF Im-0 THEN 

BEEP 

INPUT ’GIVE A NAME FOR THE RAW OATA FILE' ,02_f l la§ 

PRINT USING ’ 1 5 X , ’ " Nau file name: ” , 1 4A ’ ; 02_f l la* 

5 i za 1 -20 

CREATE BOAT D2_f l let ,Si zel 
ASSIGN ®File2 TO 02_filaS 

OUMMY FILE UNTIL Nrun KNOWN 
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560 01 _f 1 la*-*0UMMY • 

570 CREATE BOAT 01 _f l US .Size 1 
580 ASSI6N OFilal TO Ol.filaS 
590 OUTPUT ®Fi l a 1 t Oa t at 
600 IF Ihn-0 THEN 
610 BEEP 

620 INPUT 'ENTER NUMBER OF OEFECTIUE TCS (0-OEFAULT, 4 MAX.J'.Idtc 
630 IF Idtc-0 THEN 
640 Ldtcl-0 

650 Ldtc2-0 

651 Ldtc3-0 

652 Ldt c 4-0 

660 PRINT USIN6 ' 1 6X , * * No defecttvc TCs axisf- 

670 ENO IF 

680 IF Idtc-1 THEN 

690 8EEP 

700 INPUT 'ENTER OEFECTIVE TC LOCATION <l-8)",Ldtc1 

710 PRINT USING '16X,**TC is defective at location "“,00'iLdtcl 

720 Ldtc2-0 

721 Ldtc3*0 

722 Ldtc4-0 
730 ENO IF 

740 IF Idtc-2 THEN 
750 8EEP 

760 INPUT 'ENTER DEFECT I UE TC LOCATIONS < 1 -8 >" .Ldtd ,Ldtc2 

770 PRINT USIN6 " 1 6 X , * * TC are defective at locations " * ,00 , 4X , 00 * { Ldtc 1 ,Ldt c2 

771 Ldt c 3-0 

772 Ldt c 4-0 

780 ENO IF 

781 IF Idtc-3 THEN 

782 SEEP 

783 INPUT 'ENTER DEFECTIUE TC LOCATION < 1 -8 ) " ,Ldt c 1 .Ldt c2 ,Ldtc3 

784 PRINT USIN6 '1GX.-"TC is defactive at location * ■ ,3< 4X ,00 ) “ tLdtc 1 ,Ldtc2 ,Ld 
tc3 

7B5 Ldtc4-0 

7B6 ENO IF 

787 IF Idtc-4 THEN 

788 8EE D 

789 INPUT 'ENTER OEFECTIUE TC LOCATION ( 1 -8 ) “ ,Ldt c 1 ,Ldt c 2 ,Ldt c3 ,Ldt c 4 

790 PRINT USIN6 -IBX.-'TC is dafective at location “ ’ , 4< 4X ,00 ) " t Ldt c I ,Ldt c2 ,Ld 

tc3 ,Ldt c4 

792 ENO IF 

793 IF Idtc>4 THEN 

800 BEEP 

B 1 0 PRINTER IS 1 

B20 8EEP 

830 PRINT 'INUALIO ENTRY- 

840 PRINTER IS 701 

850 60T0 610 

B60 ENO IF 

870 ENO IF 

880 OUTPUT CFtlal iLdtcl ,Ldtc2 ,Ldtc3 ,Ldtc4 
B90 1 In-1 option 
900 ELSE 

910 8EEP 

920 INPUT -6IUE THE NAME OF THE EX ISTIN6 OATA r lLE* ,02_f lie# 

930 PRINT USIN6 -16X.*“01d file nane: , 1 4 A " i 02_f l 1 eS 

940 ASSIGN #Fila2 TO 02_flleS 
950 ENTER «Ftla2iNrun 

960 ENTER •File2iOoldt 

970 PRINT USIN6 '16X,''This data sat taxen on : " * , 1 4A ' i Oo l dS 

980 ENTER 0F i la2 iLdtc 1 ,Ldtc2 .Ldt c3 ,^dt c4 

990 IF Ldt c 1 >0 OR Ldtc2>0 OR Ldtc3>0 OR Ldtc4>0 THEN 

1000 PRINT USING ' 1 6X , ' ' Thernocoup les were defective at locat ions i " " ,4( 30 ,4X ) ' i 
^dtcl .Ldt c2 ,Ldtc3 ,Ldtc4 
'010 ENO IF 
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1020 ENTER ©F neZiltt 
1030 END IF 
1040 Idtc-0 

1050 IF Ldt c 1 >0 THEN Idtc-Idtc+i 

1 060 IF Ldtc2>0 THEN Idtc-Idtc+1 

1051 IF Ldtc3>0 THEN Idtc-Idtc+1 

1 052 IF Ldtc4>0 THEN Idtc-Idtc+1 

10701 IF In-0 AND Ihn-1 THEN 1595 
1 060 BEEP 

1090 INPUT ‘WANT TD CREATE A PLOT FILE' 5 ( 0-N , 1 -Y ) " , I p lot 
1100 IF Iplot-1 THEN 

1110 BEEP 

1120 INPUT “GIVE NAME FOR PLOT FILE" ,P_ f i leS 

1 130 CREATE BOAT P_f i lef ,4 

1 U0 ASSIGN ©Plot TO P.fiieS 

1 150 END IF 

1160 IF Ihn- I THEN 

1170 BEEP 

1180 INPUT "WANT TD CREATE Uo FILE 7 < 0-N , 1 -Y ) " , Iu f 
1190 IF Iuf-I THEN 
1200 BEEP 

1210 INPUT "ENTER Uo FILE NAME" ,Uf i leS 
1220 CREATE BOAT Ufile*,4 
1230 ASSIGN ©Ufile TD Ufile* 

1240 END IF 
1250 BEEP 

1260 INPUT "WANT TO CREATE Re FILE 7 ( 0-N , 1 -Y > " , Ire 
1270 IF Ire-1 THEN 
1 2B0 BEEP 

1290 INPUT "ENTER Re FILE NAME " ,Re f i leS 

1300 CREATE BOAT RefileS.10 

1310 ASSIGN ©Ref i le TO RefileS 

1320 END IF 

1330 END IF 

1340 PRINTER IS 1 

1350 IF In-0 THEN 



1360 


BEEP 










1370 


PRINT 


USING 


"AX , 


""Select tube nunder""" 


13B0 


IF Ihn-0 THEN 






1390 


PRINT 


USING 


* 6X , 


" "0 


Snoot h 4 inch Re f * " “ 


1 400 


PRINT 


USING 


"6X . 


•" 1 


Snooth 4 inch Cu ( Press/S 1 1 de > " * ' 


1410 


PRINT 


USING 


* 6X . 


""2 


Soft Solder 4 inch Cu""' 


1420 


PRINT 


USING 


*6X . 


""3 


Soft Solder 4 inch HIGH FLUX “ * " 


1430 


PRINT 


USING 


’ 6X . 


""4 


Uieland Hard B inch*"" 


1 440 


PRINT 


USING 


*6X , 


""5 


HIGH FLUX B inch""" 


1450 


PRINT 


USING 


•6X , 


""6 


GEUIA-K 40 Fins/ in" * " 


1460 


PRINT 


USING 


*6X , 


" "7 


GEUA-K 26 Fins/ in* * " 


1 470 


PRINT 


USING 


’ 6X , 


""B 


GEUA-T 19 Fi ns/ i n" * * 


1 480 


PRINT 


USING 


•6X , 


""9 


GEWA-T DR GEUA-TY 26 Fins/in‘” 


1490 


PRINT 


USING 


’ 6X , 


""10 


THERMDEXCEL-E* * " 


1500 


PRINT 


USING 


*6X , 


""11 


THERMDEXCEL-HE ' * " 


1510 


PRINT 


USING 


•6X , 


""12 


TURBD-B" • " 


1520 


PRINT 


USING 


•8X . 


" " 13 


GEWA-K 19 Fins/ l n' ’ • 


1530 


ELSE 










1 540 


PRINT 


USING 


* 6X , 


"•0 


Snooth tube""* 


1550 


PRINT 


USING 


"6X , 


" • 1 


High Flux 


1560 


PRINT 


USING 


*6X , 


"•2 


Turdo-B”* 


1570 


PRINT 


USING 


‘6X , 


•’3 


Hi Qh Flux Mod*"’ 


1 5B0 


PRINT 


USING 


*6X , 


" "4 


Turdo-B Mod" ” 


1590 


END IF 








1600 


INPUT 


Itt 









1610 OUTPUT ©Filel iltt 

1620 END IF 

1630 PRINTER IS 701 

1640 IF 1 1 1 < 1 0 THEN PRINT USIN6 "16X,’"Tube Nunber: •",D"iItt 

’650 IF 1 1 1 >9 THEN »RINT USIN6 - 16X . ’"Tube NunDer: “ " ,DD" ; Itt 



m 



1 660 IF I nm- 1 TntN PRINT uSiNo Ioa, tuoe lype- , on ; ( n-» 

7 670 SEEP 

1580 INPUT 'ENTER OUTPUT UERSION ( 0-LON6 . 1 -SHORT ,2-NONE )', Iov 
1680 0EEP 

1700 INPUT 'SELECT ( 0-LIQ , 1 -UAP ,2-< LIQ+UAP >/2 ) ' , I lqv 
17101 

1720 DIMENSIONS FD TESTEO TU0ES 
1730' ELECTRIC HEATED MODE 

17401 Di-Diameter at thermocouple positions 

1750 OATA .01 1 1 125 , .01 1 1 125 , .01 1 1 125,. 01 29540 ..012446 , . 0 1 29540 ,. 0 1 00965 
1760 DATA .0100965 , .01 1 57 , . 0 1 1 57 , . 0 1 1 57 , . 0 I 1 57 , , 0 1 1 57 , .0100965 
1 770 READ 01 a( • > 

1780 Dl-Dla(Itt) 

1790' 

1800' 02-Diametar of test section to the base of fins 
1810 DATA .015875 . .015875 , .015875 . .015824 , .01 5875 . . 01 5824 . .0 1 270 
1820 DATA .01 27 ,.01 38 ,.01 38,. 01 38,. 01 38.. 01 38 ..0127 
1 030 READ D2a( ♦ ) 

18401 

1850! 0i*Inside diameter of unanhanced ends 

1060 DATA .0127 . .0127, .0127..0132 . .0127, .0132 . .01 1 1 12S..0I 1 1 125 
1870 DATA .01 18 . .01 1 8 . .01 18 . .011 0 ,.01 1 0 . .01 1 1 125 
1 880 REAO DiaC • ) 

1090' 

1900' Oo-Outside diameter of unenhanced ends 

1910 DATA .015075 ,.015075 , .015075 , .015024 , .015075, .015024 ,.01 270 ,.01270 
1920 DATA .01331 ,.01331 ,.01331 ,.01331 , .01 50 . -01 27 
1 930 REAO 0oa( * ) 

19401 

1950' L-Length of enhanced surface 

1 960 DATA .1016, .1016,. 1016 ,. 1016 ,.2032,. 2032 ,.2032 ,.2032 ,.2032 ,.2032 ,.2032 ,.20 
32,. 2032, .2032 
1 970 READ La< ♦ ) 

19801 

1990' Lu*Langth of unenhancad surface at the' ends 

2000 DATA .0254 , .0254 , .0254 , .0254 , .0762 . .0762 , .0762 , .0762 ,.0762 , .0762 , .0762 , .07 
62 ,.0762 ,.0762 
2010 READ Lua( • ) 

2020' 



2030' 

2040 

2050 

2060 

2070* 

2000 ' 

2090' 

2100 

2M0 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 



Kcu-Thermal Conductivity of tuba 

DATA 398 ,344 ,344 . 45,344 .45 ,344 ,344 .398,398 ,398 ,398 ,398 ,344 
READ Kcua( * ) 

IF Ihm*1 THEN 

Data statements for uatar heating moda 

DATA 0.015075 ,0.015075 ,0.0169,0.0138,0.0169,0,0,0,0,0 
READ 02a( • ) 

DATA 0.0127 ,0.0127 ,0.01 45 ,0.01 27 ,0.01 45 ,0 ,0 ,0 ,0,0 
READ 0 ia( • ) 

DATA 0.015075 ,0.015875 ,0.0169,0.015875 ,0.0169,0,0,0,0,0 
READ 0oa( ♦ > 

DATA 0.3048 ,0.3048 ,0.3048 ,0.3048 ,0.3048 ,0 ,0 ,0 ,0 ,0 
READ La( ♦ ) 

OATA 0.0254 ,0.0254 ,0.0254 ,0.0254 ,0.0254 ,0 ,0 ,0 .0 ,0 
READ Lua< * ) 

DATA 398,45,398,45,398,0,0,0,0,0 
READ Kcua( « ) 

END IF 
02-02a( Itt ) 

D l *D l a ( Itt ) 

Do*Doa( Itt) 

L-La< Itt > 

Lu-Lua< Itt ) 

Kcu-Kcua (Itt) 

Xn- . 8 
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2300 rr- . 3 

2310 IF Itt-0 THEN Cf-1 .70E+9 

2320 IF ltt>0 THEN C f-3 . 7037E* 1 0 

2330 A-PI * ( 0o*2-0i * 2 )/4 

2340 P-P I *0o 

2350 IF Ihm-1 THEN 

2350 8EEP 

2370 INPUT * TU8E INITIATION MOOE. ( 1 -HOT WATER ,2-STEAM ,3-COLO WATER)', Itim 
2380 IF Itim-1 THEN PRINT USING 'IGX/'Tube Initiate: Hot Water'" 

2390 IF It im-2 THEN PRINT USING '16X,"Tube Initiate: Steam"' 

2400 IF Itim-3 THEN PRINT USING '1GX,"Tube Initiate: Cold Water'" 

2410 INPUT 'TEMP/VEL MOOE: ( 0-T-CONST ,V-0EC i 1 -T-OEC .U-CONST ; 2-T-INC .U-CONST )' , 
I tv 

2420 IF I tv-0 THEN PRINT USING ' 1 6X , ' ' Temp /Ve 1 Mode: Constant /Decrees ing" ' 

2430 IF Itv-1 THEN PRINT USING ’ 1 GX , "Temp/Ve 1 Mode: Decreasing/Constant"* 

2440 IF Itv-2 THEN PRINT USING ' 1 5X , ' ' Temp/ Ve 1 Mode: Increasing/Constant'" 

2450 INPUT 'WANT TO RUN WILSON PLOT? < 1 -Y ,0-N ) ' , Iwi 1 

24S0 IF Ihm-1 ANO Iuil-0 THEN 

2470 IF Itt-0 THEN Ci-.032 

2480 IF Itt-1 OR Itt-3 THEN Ci-.0S9 

2490 IF Itt-2 OR Itt-4 THEN C1-.052 

2500 BEEP 

2510 INPUT 'ENTER Cl ( OEF : WH- . 032 ,HF- . 059 ,TB- . 0G2 ) ' ,C i 
2520 PRINT USING ' 1 6X , " S l eder-Tat e "' 

2530 PRINT USING '1SX," Constant - * ' . 2 . 40* ; C l 
2540 ENO IF 
2550 ENO IF 

2560 IF Ihm-1 ANO Im-1 ANO Iuil-1 THEN 

2570 IF Itt-0 THEN Cl-. 032 

2580 IF Itt-I OR Itt-3 THEN Ci-.0S9 

2590 IF Itt-2 OR Itt-4 THEN Ci-.0G2 

2G00 ASSIGN ©File2 TO ♦ 

2G10 CALL W i 1 son< C f ,Ci ) 

2G20 ASSIGN ©File2 TO 02_fileS 

2530 ENTER ©F 1 1 ©2 5 Nrun ,0o 1 dS ,Ldt c I ,Ldt c2 ,Ldtc3 ,Ldt c4 , 1 1 1 
2540 ENO IF 
2550 Nsub-0 

2550 IF I dp-4 THEN Ihm-1 

2570 IF Ihm-1 THEN NsuD-8 

2580 Ntc-G 

2590 IF Ihm-0 THEN Ntc-12 

2700 J-l 

2710 Sx-0 

2720 Sy-0 

2730 Sxs-0 

2740 Sxy-0 

2’ 7 50 Repeat: i 

2750 IF Im-0 THEN 

2780 Ido-2 

2 ’’90 ON KEY 1,15 RECOUER 2750 
2800 PRINTER IS 1 

2810 PRINT USING *4X, "SELECT OPTION'" 

2820 PRINT USIN6 * SX , ' * 0-TAKE OATA"' 

2830 IF Ihm-0 THEN PRINT USING 'GX,''1-SET HEAT FLUX'" 

2840 IF Ihm-1 THEN PRINT USING 'GX," 1-SET WATER FLOW RATE*" 

2850 PRINT USING 'SX, "2-SET Tsat'" 

2860 PRINT USING *4X,"N0TE: KEY 1 - ESCAPE"' 

2870 8EEP 

2880 INPUT Ido 

2890 IF Ido>2 THEN Ido-2 

2900 IF Ido-0 THEN 4440 

2910' 

29201 LOOP TO SET HEAT FLUX OR FLOWMETER SETTIN6 
2930 IF I do- I THEN 

2940 IF Ihm-0 THEN 
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2950 OUTPUT 709»*AR AF20 AL2l UR5 * 

2960 BEEP 

2970 INPUT 'ENTER DESIRED Qdp'.Dqdo 

29B0 PRINT USING *4X, "DESIRED Odp ACTUAL Qdp*" 

2990 Err- 1000 

3000 FOR 1-1 TO 2 

3010 OUTPUT 709 i ’AS SA' 

3020 Sun-0 
3030 FOR Ji-1 TO 5 
3040 ENTER 709 i E 
3050 Sun-Sun+E 
3060 NEXT Ji 

3070 IF 1-1 THEN Uolt-Sun/5 
3080 IF 1-2 THEN Anp-Sun/5 

3090 NEXT I 

3091 PRINT *Ua : * ,Uolt 

3092 PRINT “ I 9 : “ ,Anp 
3100 Ano-ABS( Ano« 1 . 9 1 82 ) 

31 10 Uolt-ABS( Uolt*25 ) 

3120 Aqdp-Uol t • Anp/ ( P I *D2«L ) 

3130 IF ABS( Aqdo-Oqdp >>Err THEN 

3140 IF Aqdp>Dqdp THEN 

3150 BEEP 4000,. 2 

3160 BEEP 4000,. 2 

3170 BEEP 4000,. 2 

3180 ELSE 

3190 BEEP 250 ,.2 

3200 BEEP 250 ,.2 

3210 BEEP 250, .2 

3220 END IF 

3230 PRINT USING ' 4X ,MZ . 3DE ,2X ,MZ . 3DE * iDqdo ,Aqdp 

3240 WAIT 2 

3250 GOTO 3000 

3260 ELSE 

3270 BEEP 

3280 POINT USING * 4X ,MZ . 30E ,2X ,MZ . 3DE * ; Dqdp , Aqdp 

3290 Err -500 

3300 WAIT 2 

3310 GOTO 3000 

3320 END IF 

3330 ELSE 

3340 BEEP 

3350 INPUT 'ENTER FLOWMETER SETTING*, Fna 

3360 60T0 2B10 

3370 END IF 

33B0 END IF 

3390* 

3400) LOOP TO SET Taat 
3410 IF Ido-2 THEN 
3420 IF Ikdt-1 THEN 3470 
3430 BEEP 

3440 INPUT 'ENTER DESIRED Tsat’.Otld 

3450' PRINT USING *4X," DTeat ATaat Rato Tv Rate"* 

34S0 Ikdt-1 

3470 Oldl-0 

3480 Old2-0 

3490 Nn-1 

3500 Nr 9-Nn MOD 15 

3510 Nn-Nn+l 

3520 IF Nr 9 - 1 THEN 

3530 IF Ihn-0 THEN PRINT USING *4X," Teat T 1 d 1 Tld2 Tv Taunp 

3540 IF Ihn-1 THEN PRINT USIN6 *4X," Toat T 1 d 1 Tld2 Tv Taunp 
t Tpile Tout'** 

3550 END IF 

3560 IF Ihn-0 THEN OUTPUT 709i'AR AF08 ALM UR5 * 

3570 IF Ihn-1 THEN OUTPUT 709; "AR AF0 AL5 UR5 " 



T mlo 
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3580 FOR 1-1 TO 6 

3590 IF Ihn-0 AND I > 4 THEN 38G0 

3G00 Sun-0 

3610 OUTPUT 709 i “AS SA" 

3620 FOR Ji-1 TO 20 
3G30 ENTER 709iEliq 
3640 Sun-Sun+El iQ 
3G50 NEXT Ji 
3GG0 Eliq-Sun/20 
3G70 Tl d-FNTvs v( E 1 iq ) 

3G80 IF 1-1 THEN Tldl-Tld 

3G90 IF 1-2 THEN Tld2-Tld 

3700 IF 1-3 THEN Tv-Tld 

3710 IF 1-4 THEN Tsunp-Tld 

3720 IF 1-5 THEN Tinlet-Tld 

3730 IF I-G THEN Tout-Tld 

3740 NEXT I 

3750 IF Ihn-1 THEN 

37G0 OUTPUT 709 i " AR AF00 AL00 UR5 * 

3770 OUTPUT 709 t ’AS SA“ 

3780 5un-0 

3790 FOR Kk-l TO 20 

3800 ENTER 709 ; E 

3810 Sun-Sun+E 

3820 NEXT Kk 

3830 Enf (7 >-A0$< Sun/20 ) 

3840 Tp:le-Enf(7)/3.9GE-4 

3850 END IF 

38G0 At 1 d-( T Idl +Tld2 >* .5 

3870 IF A9S ( At 1 d-Dt 1 d ) > . 2 THEN 

3380 IF Atld>Dtld THEN 

3890 BEEP 4000 ,.2 

3900 8EEP 4000 ,.2 

3910 8EEP 4000, .2 

3920 ELSE 

3930 BEEP 250 , .2 

3940 8EEP 250 ,.2 

3950 8EEP 250,. 2 

39G0 END IF 

3970 ErM-Atld-Oldl 

3980 Oldl-Atld 

3990 Err2-Tv-01d2 

4000 0 1 d2-Tv 

4010 IF T 1 d 1 >1 00. THEN 4060 

4020 IF Ihn-0 THEN PRINT USING “ 4X ,5( P1DDD .DD ,2X )’ iDtld.Tldl ,Tld2 ,Tv , Tsunp 

4030 IF Ihn-1 AND Idp-0 THEN PRINT USING * 4X , 7{ HDD . DD ,2 X ) * i D 1 1 d ,T 1 d 1 , T 1 d2 ,T v , Ts 

unp ,T inlet ,Tpi le 

4040 IF Ihn-1 AND Idp-4 THEN PRINT USIN6 ' 4X ,5( MOD .DD ,2X ) ,3< H3D.DD ,2X )" iDt id ,T1 
dl ,Tld2 ,Tv .Taunp ,T inlet , Tp; le .Tout 



4050 


WAIT 2 




4060 


SOTO 3500 




4070 


ELSE 




4080 


IF ABS< At ld-Dt Id )> . 1 


1 THEN 


4090 


IF A t 1 d>Dt 1 d THEN 




4100 


8EEP 3000,. 2 




41 10 


BEEP 3000 2 




4120 


ELSE 




4130 


SEEP 800,. 2 




41 40 


8EEP 800 , .2 




4 150 


END IF 




4160 


Err 1 -At ld-Oldl 




4 170 


Oldl-Atld 




4180 


Err2-Tv-01d2 




4190 


01 d2*Tv 




4200 


IF Ihn-0 THEN PRINT 


USIN6 " 4X ,5( MDDD . DD ,2X )’ ;Dtld,Tldl ,Tld2 ,Tv .Tsunp 
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4 z 1 0 IP Ihn-1 THEN PRINT USING ‘ 4X ,5< MDD.DO ,2A > ,3< H30 . DD , i x " ; 0 1 1 d , ~ 1 d 1 , T i d 2 , T v 
, Tsunp ,T iruai .ipiie , i out 
4120 WAIT 2 
4230 60TO 3500 
4240 ELSE 
4250 BEEP 

4260 Err)- Atld-Oldl 
4270 0 1 dl -At Id 

4280 Err2-T v-01 d2 
4290 0Ld2-Tv 

4300 IF Ihn-0 THEN PRINT USING * 4X ,5( ttDDD . D0,2X)"iDtld,Tldl ,Tld2,Tv,T suno 

4310 IF Ih n-1 THEN PRINT USING * 4X ,8 ( MOO . DD ,2X )*;DtId,Tldl . T 1 d 2 ,Tv .Tsuno ,T i nle t 

,Tpi 1 e .Tout 
4320 WAIT 2 
4330 GOTO 3500 
4340 END IF 
4350 END IF 
4360 END IF 

4370' ERROR TRAP FDR Ido OUT OF BOUNDS 

4380 IF I do>2 THEN 

4390 BEEP 

4400 GOTO 2810 

4410 END IF 

4420 1 

44301 TAKE OATA IF In-0 LOOP 
4440 IF Ikol-1 THEN 4480 
4450 8EEP 

4460 INPUT ‘ ENTER BULK OIL X* .Boo 
4470 Ikol-1 

4480 IF Ihn-0 THEN OUTPUT 709i‘AR AF00 AL M UR5 " 

4490 IF Ihm-1 THEN OUTPUT 709j'AR AF 0 ALS UR5* 

4500 IF Ihm-0 THEN Ntc-12 

4510 FOR I-l TO Ntc 
4520 OUTPUT 709! *AS SA' 

4530 5un-0 

4540 FOR Ji-1 TO 20 ’ 

4550 ENTER 709;E 
4560 Sun-Sun*E 

4570 IF I- ( 9-Nsuo ) OR I-U0-Nsub) THEN Et ( J i-l )-£ 

4580 NEXT Ji 
4590 Kdl-0 

4600 IF I-( 9-Nsub ) OR I •( 1 0-Nsub ) THEN 
4610 Eave-Sun/20 
4620 Sun-0. 

4630 FOR Jk-0 TO 19 

4640 IF ABS(Et( Jk >-Eave X5.0E-6 THEN 

4650 Sun-Sun+Et ( Jk ) 

4660 ELSE 
4670 Kdl-Kdi+1 
4680 END IF 
4690 NEXT Jk 

4700 IF I-( 9-Naub ) OR I-M0-Naub> THEN PRINT USING *4X,”Kdl - "‘.DO’iKdl 
4710 IF Kdl> 1 0 THEN 
4720 BEEP 
4-30 BEEP 

4740 PRINT USIN6 *4X,*'Too ciucpi scattering in data - repeat data aat*"“ 

4750 GOTO 2800 
4760 ENO IF 
4770 END IF 

4780 Enf ( I >-Sun/( 20-Kdl ) 

4* , 90 NEXT I 

4800 IF lhn-1 THEN 

4810 OUTPUT 709 i * AR AF00 AL00 UR5 ’ 

4820 OUTPUT 709 i ‘AS SA - 

4830 Sun-0 

4840 FOR Kk-1 TO 20 
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4352 EN’E^ ' , 05;E 

4552 S^n-Sun+E 

4870 NEXT Kk 

4880 En f ( 7 )-ABS( Sun )/20 

4890 ENO IF 

4900 IF Ihn-0 THEN 

4910 OUTPUT 709 i * AR AF20 AL2 1 UR5 m 

4920 FOP 1-1 TO 2 

4930 OUTPUT 709 j 'AS SA' 

4940 Sun-0 
4950 FOP Ji-1 TO 5 
4960 ENTER 709 j E 
4970 Sun-Sun+E 
4980 NEXT Ji 

4990 IF 1-1 THEN Vr-Sun/5 
5000 IF 1-2 THEN Ir-Sun/5 
5010 NEXT I 
5020 ENO IF 
5030 ELSE 

5040 IF Ihn-0 THEN ENTER ©F i le2 ; Boo , To Id* ,En f (•) ,Vr , Ir 
5050 IF Ihn-1 THEN ENTER ©Fi Ie2 ;Bop ,Told$ .Er^f ( • > ,Fna 
5060 ENO IF 
50701 

5080' CONVERT enf’S TO TEMP , VOLT .CURRENT 

5090 Tua-0 

5100 FOR 1-1 TO Ntc 

5110 IF I atc>0 THEN 

5120 IF I -Ldtc 1 OR I-Ldtc2 OR I-Ldtc3 OR I-Ldtc4 THEN 

5130 T< I > — 99.99 

5140 60T0 5240 

5150 ENO IF 

5160 ENO IF 

5170 IF 1 1 1 < 4 ANO Ihn-0 THEN 

5180 IF I > 4 ANO I < 9 THEN 

5190 T(I> — 99.99 

5200 GOTO 5240 

5210 ENO IF 

5220 ENO IF 

5230 T ( I >-FNT v a v ( En f ( I ) ) 

5240 NEXT I 
5250 IF Itt <4 THEN 

5251' UAPNIN6: 4 BAO THERMOCOUPLES ARE NOT POSSIBLE FOR ITT<4 

5260 FOR I-l TO 4 

5270 IF I-Ldtc 1 OR I-Ldtc2 THEN 

5280 Twa-Twa 

5290 ELSE 

5300 Tua-Twa + T( I ) 

5310 ENO IF 
5320 NEXT I 
5330 Tua-Twa/( 4-Idtc ) 

5340 ELSE 

5350 IF Ihn-1 THEN 5450 
5360 FOP 1-1 TO 8 

5370 IF I-Ldtc 1 OR I-Ldtc2 OR I-Ldtc3 OR I-Ldtc4 THEN 

5380 Twa-Tua 

5390 ELSE 

5400 Tua-Tya+T(I) 

5410 ENO IF 
5420 NEXT I 
5430 Tw-Twa/ ( 8-1 dt c ) 

5440 ENO IF 
5450 T ld-T ( 9-Noub ) 

5460 Tld2-T( 10-Naub ) 

5470 TIda-< TId+Tld2 )*.5 
5480 Tv-T( 1 1-Naub ) 

5490 IF 1 1 1<3 ANO Ihn-0 THEN 
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5500 Tld2--99.99 

55 l 0 T v< T< 1 0 )♦ T< ) ! > )/2 

5520 END IF 

5530 Tsunp-T ( \ 2-Naub > 

5540 IF Ihn-0 THEN 5570 
5550 Tinlat-H 13-Naub ) 

55G0 Tou t -T< M -N sub > 

5570 IF Ihn-0 THEN 
5590 Anp-A8S< Ir*l .9182 ) 

5590 Uolt-A8SCUr >*25 
5600 Q-Uo 1 1 *Anp 
5G10 END IF 

5G20 IF Itt-0 AND Ihn-0 THEN 
5G30 Kcu-FNKcu< T y > 

5640 ELSE 

5G50 Kcu-Kcua< 1 1 1 ) 

5GG0 END IF 
5670 i 

5680' FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
5G90 IF Ihn-0 THEN Ty-Ty-Q«L06( D2/D 1 >/ < 2 *P I *Kcu *L > 

5700 IF 1 1 q v — 0 THEN Taat-Tlda 

5710 IF Ilqv-1 THEN Tsat-^ Tlda+Tv >♦ . 5 

5720 IF Ilqv-2 THEN Tsat-Tv 

5730 IF I hn- ! THEN 

5740 Tavg-Tinlet 

5750 6r ad- 37. 9853+. 104388*Tavg 

57G0 T drop-A8S( Enf ( 7 ) >*1 .E*6/< I0*6rad> 

5770 Tavgc-T tnlet-Tdrop* . 5 

5780 IF A8S< Tavg-Ta vgc >> . 0 1 THEN 

5790 Tavg-( Tavg+Tavgc )* .5 

5800 GOTO 5750 

5810 END IF 

5820' 

5830' COMPUTE WATER PROPERTIES 
5340 IF Ihn-1 THEN 
5850 K w-FNK y < Ta vg ) 

5860 Muwa-FNMuy ( Tavg ) 

5870 Cpy-FNCpy ( Tavg ) 

5880 Pry -FNPry ( Ta vg ) 

5990 Rhoy-FNRhpy ( Ta vg ) 

5900 Ty i -T avg 
5910' 

5S20 1 Compute MDDT 

5930 Mdot-3.9657E-3*Fns*< 3 .6 I 9S5E-3-Fns*< 8 . 82006E-6-Fma*< 1 . 23688E-7-Fna *4 . 3 1 897 

E- 1 0 > ) > 

5940' Mdot-Mdo t •< 1 . 0365 -T inlet *< 1 . 96644E-3-Tinlet*5. 252E-6 > >/ 1 . 0037 
5950 Kdt-0 

5960 Q-Mdot *Cpy*Tdrop 

5970 Lntd-Tdrop/LD6< < T mlet-Taat )/<Tinlat-Tdrop-Taat > > 

5980 Uo-G/(PI*Do*L*Lntd> 

5990 Ry-Do*LD6(Do/Di >/<2.*Kcu> 

6000 Ty-Taat+Fr*Lntd 

6010 Uy-Mdot/< Rhoy*PI*Dl "2/4 ) 

6020 Rey-Rhoy # Uw*Di/Muya 

6030 Hi-Ci*Ku/Dl*Rey".8*PryM 1/3. >•( Muya/ FNM uu( T ui > )* . 14 
6040 Ty ic-Tavg-Q/ ( P I *Do # L*Hi ) 

6050 IF A8S< Ty i-Tuic >> . 0 1 THEN 

6060 Ty l • ( Ty l ♦Ty ic > * . 5 

6070 60TD 6030 

6080 END IF 

6090 Tyi-( Tyi *Ty ic >• . 5 

6100 Ho-1 /< 1 /Uo-Do/< Di*Hi )-Ry ) 

6110 END IF 

6120 END IF 

6130 IF I hn- l THEN 

6140 ThetaP-Q/ ( Ho*P I *Do*L > 



118 



6150 Tw-Taat+Tneiao 
6 1 62 ElSE 

6)70 Thetab-Tw-Taat 
6180 ENO IF 
6190 IF Thetab<0 THEN 
6200 8EEP 

6210 INPUT "TUALKTSAT (0-CONTINUE, l-ENOT.Iev 

6220 IF Iev-0 THEN GOTO 27G0 

6230 IF lev-! THEN 7220 

6240 ENO IF 

G250I 

6260 * COMPUTE VARIOUS PROPERTIES 
6270 T f i In- ( Tw+T aat )*.5 
6280 Rho-FNRho( Tf t In ) 

6290 Mu-FNMy ( Tfiln) 

6300 K-FNK ( Tf l In ) 

63)0 Cp-FNCp< Tf l In ) 

6320 8eta-FN8eta( Tfiln ) 

6330 Hf g-FNHf g< Tsat ) 

6340 Ni-Mu/Rho 
6350 Alpha-K/( Rho*Cp ) 

6360 Pr-N i / A 1 phe 
6370 Psa t-FNPsa t < Taat ) 

6380' 

638) PRINT 'Rho: " ,Rho 

6382 PRINT ‘Mu: * ,Hu 

6383 PRINT * K : ’ ,K 

6384 PRINT * Cp : * ,Cp 

6385 PRINT "Tf iln: “ , Tfiln 

6386 PRINT ’Hfg: “ ,Hfg 

6387 PRINT 'Psat:* ,Psat 

6388 PRINT 'Taat : * ,Tsat 

6390' COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
64001 FOR UNENHANCEO ENO(S) 

6410 Hoar- 1 90 

6420 Fe-< Hbar*P/( Kcu*A ) ) A . 5*Lu 
6430 Tanh-FNTanh< Fe ) 

6440 The t a-Tne t ab • Tanh/F e 

6450 X x- < 9.8) *8eta»Thetao«0o A 3*Tanh/(Fe»Ni*Alpha ) ) A . 166667 
6460 Y y-< ) ♦< . 559/ Pr ) A ( 9/ 1 6 > ) “ ( 8/27 ) 

6470 Hbarc-K /Oo # ( . 6+ . 387* X x / Yy ) A 2 

6480 IF A8S( ( Hbar-Hbarc )/Hbarc )> .001 THEN 

6490 Hbar-( Hbar+Hbarc ) • . 5 

6500 GOTO 6420 

65)0 END IF 

651) PRINT ’ NAT CONV h: ' ,Hbar 
65201 

6530' COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENO(S) 

6540 Ql-< Hbar» P*Kcu*A ) A .5 # Thetab«Tanh 
6550 Qc-Q-2 «Q 1 
6560 As-PI *02*L 
65701 

6580' COMPUTE ACTUAL HEAT FLUX ANO 80ILIN6 COEFFICIENT 

6590 Qdp-Qc/Aa 

6600 Htube-Qap/Thet ab 

66 10 Caf-(Cp«Thetab/Hfg )/( Qdp/( Mu*Hfg )•( .01 4/( 9.81 *Rho ) A .5 ) A < ) /3. )*Pr A ) .7 ) 
6620' 

6630' RECORD TIME OF OATA TAKIN6 
6640 IF In-0 THEN 
6650 OUTPUT 709 i * TD" 

6660 ENTER 709;Told$ 

6670 ENO IF 
6680' 

6690' OUTPUT OATA TO PRINTER 
6700 PRINTER IS 70) 

67)0 IF Iov-0 THEN 
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6720 


PRINT 




6730' 


PRINT USING 


- 1 0X , 


4A * i J . Sop .To 1 dS 




6731 


PRINT USIN6 


- 1 0X 


op 






6732 


PRINT ’ 




6740 


IF Ihn-0 THEN 


6750 


PRINT USIN6 


• 10X 


8 


. . . 




6760 


PRINT USING 


•10X 


7 ) ,T( 8 > 




6770 


PRINT USIN6 


• 10X 


6780 


PRINT USIN6 


- 1 0X , 


ld2 , T 


v ,Psat .Tsunp 




6790 


PRINT USIN6 


• 1 0X 


6800 


PRINT USIN6 


• 10X , 


6810 


ELSE 




6820 


PRINT USING 


■ 10X , 




Ho*-- 




6830 


PRINT USING 


■ 10X , 


nlet .Tdrop .Tnetab , 


,Qap .<■ 



“•Data Set Nunoe'* - * ‘ .ODD ,2X , * * S* i * Oi. X • '*,DD.O,5x,l 
, * ‘Data Set NunDer - *' ,DOO ,2X ‘Bulk Oil X - * , ,DD.O'iJ,B 
TIME:* ,TIME*( TIMEDATE ) 

, * " TC No: I 2 3 4 5 S 7 



Qdp * * ’ 



6840 END IF 

6850 ENO IF 

6860 IF Iov-1 THEN 

6870 IF J-l THEN 

6380 PRINT 

6890 IF Ihn-0 THEN 

That ab ' ’ " 
THETA8 * * ' 

6940 ENO IF 
6950 IF Ihn-0 THEN 

6960 PRINT USIN6 " 1 2X ,30 ,4X ,00 t 2X ,M00 . 00 ,3( 1 X ,MZ . 30E ) * i J ,8op ,Tsat .Htube ,Qdp .The 
tab 

6970 ELSE 

6980 PRINT USIN6 * 1 2X ,30 ,4X ,00 ,2X ,MDO . 00 ,31 1 X ,nZ . 30E ) * ;Fns ,8oo ,Tsat , Htube ,Qdo ,T 



6900 

6910 


PRINT USIN6 ' 10X RUN No 

ELSE 


Oils 


Ts at 


Htube 


Odp 


6920 

6930 


PRINT USIN6 -10X FMS 

ENO IF 


OILS 


TSAT 


HTUBE 


QDP 



he t ao 

6990 ENO IF 
7000 ENO IF 
7010 IF In-0 THEN 
7020 8EEP 

7030 INPUT 'OK TO STORE THIS OATA SET C 1 - V , 0-N ) 7 ’ ,0k. 

7040 ENO IF 

7050 IF Ok - 1 OR In-1 THEN J-J + 1 
7060 IF Ok- 1 AND In-0 THEN 

7070 IF Ihn-0 THEN OUTPUT ®F 1 1 e 1 i Boo , To 1 dS ,En f ( • > ,Ur , Ir 
7080 IF Ihn-I THEN OUTPUT ©F 1 1 e 1 « 8oo . To 1 d$ , Enf ( • ) ,Fns 
7090 ENO IF 

7100 IF Iuf-1 THEN OUTPUT ©UfiletUw.Uo 
7110 IF Ire-1 THEN OUTPUT ©Re f 1 1 e « Fns ,Re w 

7120 IF (In-1 OR Ok- 1 ) ANO Iolot-I THEN OUTPUT ©P lot iQdp ,The tab 
7 1 30 IF In-0 THEN 
7140 8EEP 

7150 INPUT ’WILL THERE BE ANOTHER RUN ( 1 -Y , 0-N ) 7 * ,6o_on 
7160 Nrun-J 

7170 IF 6o_on-0 THEN 7220 
7180 IF 6o_on<>0 THEN Repeat 
7190 ELSE 

7200 IF J<Nrun*l THEN Repeat 
7210 ENO IF 
7220 IF In-0 THEN 
7230 8EEP 

7240 PRINT USIN6 ‘10X,’ - NOTE: -",ZZ.’" data runs were stored in file "-.10A"iJ 
1 ,02_f lief 

7250 ASSIGN ©Filel TO * 
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7260 OUTPUT 0F 1 1 e2 1 Nrun- 1 

7270 ASSIGN ©Fi lei TO Ol.flieS 

7280 ENTER OF 1 1 e 1 ; Oat e* ,Ldt c 1 ,Ldtc2,Itt 

7290 OUTPUT •File2iOateS,Ldtcl ,Ldtc2,Itt 

7300 FOR 1-1 TO Nrun-1 

7310 IF Ihn-0 THEN 

7320 ENTER ®Fi le I jBod , T oldS , Emf < •) ,Ur , Ir 
7330 OUTPUT IFile2iBop , ToldS .Enf ( • ) ,Ur t Ir 
7340 ELSE 

7350 ENTER CF i le 1 ;0op , ToldS ,E*f (•) ,Fw S 
7360 OUTPUT ®Fi le2 iBop .ToldS ,Enf(* ) .Fms 
7370 ENO IF 
7380 NEXT I 

7390 ASSIGN ©Filel TO • 

7400 PURGE “0UMMY“ 

7410 END IF 
7420 BEEP 
7430 PRINT 

7440 IF Ip lot- 1 THEN PRINT USING " 1 0X , " - NOTE : X-Y pairs were stored in 

plot data file “ “ , 1 0A “ i J- I ,P.f l le* 

7450 ASSIGN 0FiIe2 TO • 

7460 ASSIGN OPiot TO • 

7470 IF Iuf-1 THEN ASSIGN ©Ufile TO * 

7480 IF Ire-1 THEN ASSIGN GRefile TO • 

7490 CALL Stats 
7500 0EEP 

7510 INPUT 'LIKE TO PLOT DATA < 1 -Y , 0-N P ' ,0k 
7520 IF Ok - 1 THEN CALL Plot 
7530 SUBENO 
7540) 

7550' CURUE FITS OF PROPERTY FUNCTIONS 
7560 OEF FNKcu(T) 

75701 OFHC COPPER 250 TO 300 K 
7580 Tk-T+273 . 1 5 'C TO K 

7590 K-434-.112*Tk . 

7600 RETURN K 
7610 FNENO 

7620 OEF FNMu<T ) 

7 62 1 COM /Iprop/ Ift 

7622 IF Ift-0 THEN 

7630' 170 TO 360 K CURUE FIT OF VISCOUSITY 
7640 Tk-T+273 . 1 5 ' C TO K 

7650 Mu-EXP(-4.4636+< 1011 .47/Tk ) )«1 .0E-3 

7651 ENO IF 

7652 IF Ift-1 THEN 
76531 223 TO 373 K 

7657 Mu-3. 4593-4. 26E-2*T+3.9485E-4*T‘ 2-4. 1 93E-6 *T “ 3+2 . 0709E-8*T * 4 

7658 Mu-Mu* 1 . 0E-4 

7659 ENO IF 

7660 RETURN Mu 
7670 FNENO 

7680 OEF FNCp< T ) 

7681 COM /Iprop/ Ift 

7682 IF Ift-0 THEN 

76901 180 TO 400 K CURUE FIT OF Cp 
7700 Tk-T+273 . 15 ' C TO K 

7710 Cp-. 40 188+1 . 65007E-3*Tk + l . 5 1 494E-5*Tk “2 -6 . 67053E- 1 0*Tk ‘ 3 

7720 Cp-Cp • 1 000 

7721 ENO IF 

7722 IF Ift-1 THEN 

7723' 223 TO 373 K CURUE FIT Cp(Rl24) 

7725 Cp- 1 .0542+2. 1405E-3*T+ I . 0709E-5*T “ 2-6 . 472 1 E-8* T * 3- 1 . 4324E-9*T - 4*4 . 136E-1 1 * 
T"5 

7726 Cp-Cp*1000 

7727 ENO IF 
7730 RETURN Cp 
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7^40 FNEND 

7750 OEF FNRho( T ) 

7751 COM / 1 prop / Ift 

7752 IF Ift-0 THEN 

7750 Tk-T + 273 . 1 5 1C TO K 
7770 X-l-< 1 .8*Tk/753.9S ) l K TO R 

7780 Ro-36 .32*61 . 1 464 U*X * ( 1 /3 )♦ 1 6 . 41 00 1 5*X* 1 7 . 476838*X * .5+ 1 . 119828*X*2 

7790 Ro-Ro/. 062428 

7791 END IF 

7792 IF Ift-1 THEN 

7793' 223 TO 373 K (R-124 ) 

77 96 Ro-1434.8-2.861 9*T-6 .7267E-3*T-2-7.2852E-5*T*3 

7797 END IF 

7800 RETURN Ro 

7810 FNEND 

7820 DEF FNPr(T) 

7830 Pt-FnCd ( T )*FNhu(T )/FNK(T > 

7840 RETURN Pr 
7850 FNEND 

7860 DEF FNK(T) 

7861 CON /Iorop/ Ift 

7862 IF Ift-0 THEN 
78701 T <360 K WITH T IN C 

7880 K-.071-. 000261 *T 

7881 END IF 

7882 IF Ift-1 THEN 

7883' 223 TO 373 K (R-124 ) 

7885 K-7.5l9lE-2-3.5436E-4«T-1 . 9S45E-7*T “ 2+3 . 1 03SE-9*T*3 

7886 END IF 
7890 RETURN K 
7900 FNEND 

7910 DEF FNT anh( X ) 

7920 P-EXP(X) 

7930 0-1 /P 

7940 Tanh- ( P-Q ) / ( P + Q ) 

7950 RETURN Tanh 

7960 FNEND 

7970 DEF FNTvsv(U) 

7980 CDn /Cc/ C(7),Ical 
7990 T-C( 0 ) 

8000 FDR 1-1 TO 7 
8010 T-T*C< I )*U* I 
8020 NEXT I 
8030 IF leal - 1 THEN 

8040 T-T-6.7422934E-2 + TM 9. 0277043E-3-T •( -9 . 32599 1 7E-S ) ) 

8050 ELSE 

8060 T-T+8.626897E-2+T«(3.76199E-3-T*5.0689259E-5 ) 

8070 END IF 
8080 RETURN T 
8090 FNEND 
8100 DEF FN8«ta(T ) 

81 1 0 Rop-FNRho( T ♦ . 1 ) 

8120 Row-FNRho( T-. 1 ) 

8130 8ata--2 /( Rod ♦R ow )• < Rod -R ow )/ .2 

8140 RETURN Seta 
8150 FNEND 

8160 DEF FNHf g( T ) 

8161 COn /Iprop/ Ift 

8162 IF Ift-0 THEN 

8 170 Hfg-I .3741 344E + S-TM 3.309436 IE-2 *T*1 .21 65 1 43 > 

8171 END IF 

8172 IF Ift-1 THEN 

8173' 223 TO 373 K CURUE FIT (R124) 

8176 Hfg-159.93-.4609*T-1 . 458E-3*T *2- 1 .3715E-5*T*3 

8177 Hfg-Hfg* 1000 

8178 END IF 
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8 1 80 RETURN Hfg 
819(3 FNEND 

8200 DEF FNPsat ( Tc > 

8201 COM / Iprop/ Ift 

8202 IF Ift-0 THEN 

8210* 0 TO 80 deg F CURVE FIT OF Psat 
8220 T f- 1 . 8*Tc+32 

8230 Pa-5 . 945525+T f * ( . 1 S352082+T f • ( I . 4840963E-3 + T f *9 . S 1 5067 1 E-6 ) ) 

8240 Pg-Pa-14.7 

8241 END IF 

8242 IF Ift-1 THEN 

8243 1 223 TO 373 K (R-124) 

8245 Pa-( 162.68+5.946 *Tc+9.2081E-2*Tc-2+6.9359E-4*Tc*3 >*.14504 

8247 Pg-Pa-14.7 

8248 END IF 

82S0 IF Pg>0 THEN 1 +-PSI6, — in Hg 

8250 Psat-Pg 
8270 ELSE 

8280 Psat-Pg*29.92/ 14.7 

8290 END IF 
8300 RETURN Psat 
8310 FNEND 

8320 DEF FNHsnoot h< X ,8op , I sat > 

8330 DIM A( 5 ) ,8 ( 5 ) ,C( 5 > ,D( 5 > 

8340 DATA .20526 ,.25322 ,.319048 ,. 55322 7990 9 , 1 .00258 

8350 DATA . 745 1 5 , . 72992 , . 73 1 89 , . 7 1 225 , . 68472 , . 6 4 1 97 
8360 DATA . 4 1 092 , . 1 7726 ,.25 M2 . .54806 , . 8 1 9 1 6 , 1 . 0845 
8370 DATA . 7 1 403 . . 729 1 3 , .72565 , .6966 91 , .665867 ,.61889 
8380 READ A( • > ,B( * > ,C< * > ,D< • > 

8390 IF 8od< 6 THEN I-8op 
8400 IF 8 od- 6 THEN 1-4 
8410 IF 8op- 1 0 THEN 1-5 
8420 IF Isa t- 1 THEN 
8430 Hs-EXP< A( I >+8( I >*L06< X > > 

8440 ELSE 

84S0 Hs-EXP(C( I > + D( I )*L0G( X )> 

8460 END IF 
8470 RETURN Hs 
8480 FNEND 
8490 DEF FNPol y ( X ) 

8500 COM /Cply/ A( 1 0 , 1 0 ) ,C( 1 0 ) ,8 ( 5 ) ,Nop , Iprn t ,0po , I log , I f n , I j oi n ,N j o i n 
8510 XI-X 
8520 Po 1 y-8 ( 0 > 

8530 FOR 1-1 TD Nop 

8540 IF I log- 1 THEN X1-L06<X> 

8550 Poly-Poly+8( I >#X1*I 
8560 NEXT I 

8570 IF I log- 1 THEN Poly-EXP( Poly > 

8580 RETURN Poly 
8590 FNEND 
8600 SUB D oly 

8610 DIM R( 10 > ,S< 10 > ,Sy( 1 2 > ,S*< 12 > ,Xx< 100 ) ,Y y ( 1 00 ) 

8620 COM /Cply/ A( 1 0 , 1 0 ) ,CM 0 ) ,8 ( 5 ) ,N , Iprnt ,0po , I log , I f n , I j oi n ,Njoi n 
8630 CDM /Xxyy/ Xp( 25 ) ,Yp( 25 > 

8640 FOR 1-0 TD 4 
8650 8 < I >-0 

8660 NEXT I 
8670 8EEP 

8680 INPUT ‘SELECT ( 0-FILE , 1 -KEY8DARD , 2-PROGRAM >“, In 
8690 In-In+1 
8700 8EEP 

8710 input -enter numser of x-y pairs-, N p 

8720 IF In- 1 THEN 
8730 8EEP 

8740 INPUT "ENTER DATA FILE NAME * , 0 f l lei 

8750 BEEP 
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8750 

8770 

8780 

8780 

8800 

8310 

8820 

8830 

3340 

8850 

3850 

8870 

8880 

8830 

8300 

8910 

8920 

8930 

8940 

8950 

8950 

8970 

8980 

8990 

9000 

9010 

9020 

9030 

9040 

9050 

9060 

9070 

9080 

9090 

9100 

91 10 

9120 

9130 

9140 

9150 

9150 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

9240 

9250 

9260 

9270 

9280 

9290 

9300 

9310 

9320 

9330 

9340 

9350 

9360 

9370 

9380 

9390 

9400 

9410 



INPIT ‘LIKE TO EXCLUDE DAT* =>MR5 ( -v ,0-N )*> les 
IF led- 1 THEN 

8EE P 

INPUT 'ENTER NUM8ER OF PAIRS TO BE EXCLUOED' ,Ipe* 
ENO IF 

ASSI6N CFile TO O.filet 

ElSE 

8EEP 

INPUT * UANT TO CREATE A DATA FILE < 1-Y ,0-N >7' ,Yea 

IF Yea-1 THEN 

8EEP 

INPUT '6IUE A NAME FOR OATA FILE' ,D_f l leS 
CREATE 8DAT D_file$,5 
ASSI5N ©File TD D_file« 

ENO IF 
END IF 
8EEP 

INPUT 'ENTER THE ORDER OF P DLYNDHI AL " ,N 

POR 1-0 TO N*2 

Sy< I )-0 

Sx< I )-0 

NEXT I 

IF led- 1 ANO In-1 THEN 
FOR 1-1 TO Ipex 
ENTER •FileiX.Y 
NEXT I 
ENO IF 

FDR 1-1 TD Np 
IF In- 1 THEN 

IF Dpo-2 THEN ENTER ® p ileiX,Y 
IF Odo<2 THEN ENTER OFxieiY.X 
IF Ddo-1 THEN Y-Y/X 
IF Ilog-1 THEN 
IF Ddo-2 THEN Xt-X/Y 
X-l06( X ) 

IF Odo-2 THEN Y -LDG ( X t ) 

IF Ddo<2 THEN Y«LD6< Y ) 

END IF 
END IF 

IF In-2 THEN 
8EEP 

input -enter next x-y pair-,x,y 

IF Yes-1 THEN OUTPUT ©FiIe;X,Y 
ENO IF 

IF In<3 THEN 
Kx( I )-X 
Y y ( I )-Y 
ELSE 

X-Xp( 1-1 ) 

Y«Yp< 1-1 ) 

ENO IF 
R< 0 >- Y 

Sy < 0 >«Sy( 0 )♦ Y 

s< i )-x 

Sx< 1 )-Sx ( 1 > + X 
FOR J-1 TD N 
RC J )-R( J-1 > • X 
Sy< J >-Sy ( J > + R( J ) 

NEXT J 

FOR J-2 TD N*2 
S( J )-S< J-1 )*X 
Sx< J )-Sx( J )^S( J ) 

NEXT J 
NEXT I 

IF Yes-1 AND In-2 THEN 
8EEP 
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9420 PRINT USIN6 *12X,Ub," *“* pairs were siorea in 'lie , ;Np,u_nleS 

3430 ENO IF 
9440 Sx(0)-Np 
9450 FOR 1-0 TO N 
9460 C(I)-Sy(I) 

9470 FOR J-0 TO N 
94B0 A( I ,J >-Sx( Id ) 

9490 NEXT J 
9500 NEXT I 
9510 FOR 1-0 TO N-1 
9520 CALL Divided ) 

9530 CALL Subtract ( 1*1 > 

9540 NEXT I 

9550 B(N)-C(N)/A(N,N) 

9560 FOR 1-0 TO N-1 
9570 B < N— 1 — I >-C( N- 1 -I ) 

95B0 FOR J-0 TO I 

9590 B ( N- I - 1 >-B < N- 1 - 1 )- A( N- 1 - 1 ,N-J >*B(N-J ) 

9600 NEXT J 

9610 B( N-1 -I )-B( N-1 -I )/A( N-l-I ,N-1 -I > 

9620 NEXT I 

9630 ' PRINTER IS 701 

9640 i PR INI B( • > 

9650 i PRINTER IS 705 
9660 IF Iprnt-0 THEN 

9670 PRINT USIN6 " 1 2X , ' "EXPONENT COEFFICIENT'"' 

96B0 FOR 1-0 TO N 

9690 PRINT USIN6 " 1 5X ,00 ,5X ,M0.70E“ 1 1 ,B( I ) 

9700 NEXT I 
9710 PRINT - - 

9720 PRINT USIN6 '12X."0ATA POINT X Y Y( CALCULATED ) OISCR 

EPANCY'" 

9730 FOR 1-1 TO Np 
9740 Yc-B ( 0 ) 

9 7 50 FOR J-l TO N 
9760 Yc-Yc + B< J >«Xx< I >" J 
9770 NEXT J 
97B0 0- Y y ( I )- Y c 

9790 PRINT USIN6 ' 1 5X ,30 ,4X ,4( MO . SDE , 1 X > ' > I ,Xxd > , Yy d > ,Yc ,0 

9B00 NEXT I 

9B10 END IF 

9B20 ASSIGN ©File TO • 

9B30 SUBENO 

9B40 SUB Oivide(M) 

9B50 COM /Cp 1 y / A( 1 0 , 1 0 ) ,C( 1 0 ) ,B( 5 ) ,N , Iprnt ,0po , I I oq ,1 fn , I join ,N join 
9B60 FOR I-M TO N 
9B70 Ao-A( I , M ) 

9BB0 FOR J-M TO N 

9B90 A< I ,J >-A ( I ,J ) / Ao 

9900 NEXT J 

9910 C(I)-Cd)/Ao 

9920 NEXT I 

9930 SUBENO 

9940 SUB Subtree t ( K ) 

9950 COM /Cply/ A( 1 0 , 1 0 > ,C( 1 0 > ,B( 5 > ,N , Iprnt ,0po 1 1 Ioq , I f n , I join ,N join 

9960 FOR I-K TO N 

9970 FOR J — K — 1 TO N 

99B0 A( I ,J >-A( K- 1 ,J )-A( I ,J ) 

9990 NEXT J 

10000 C( I >-C ( K - 1 >-C( I ) 

10010 NEXT I 
10020 SUBENO 
10030 SUB P 1 in 

10040 COM /Cply/ A( 1 0 , 1 0 ) ,L( 1 0 ) ,B( 5 ) ,N , Iprnt ,0po , I log , I f n , I jo in ,N jo i n 
10050 COM / X x y y / Xx( 25 ) ,Yy ( 25 ) 

10060 PRINTER IS 705 
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10070 BEEP 

10080 INPUT 'UANT TO PLOT Uo vs Uw? U -Y ,0-N ) " , Iuo 
10090 IF Iuo-0 THEN 
10100 BEEP 

10110 INPUT 'SELECT <0-h/h0I sane tube , 1 -h< HF >/h< am > " , Irt 
10120 BEEP 

10130 INPUT 'SELECT h/h RATIO < 1 “FILE ,0*COMPUTED ) " , I hr at 
10140 IF Ihrat-0 THEN 
10150 BEEP 

10160 INPUT * UHICH Tsat ( 1 -6 . 7 ,0--2 . 2 > ' , I sat 
10170 ENO IF 
10180 Xm l n«0 
10190 Xmax-1 0 

10200 Xatep-2 

10210 IF Ir t -0 THEN 

10220 Ynin-0 

1 0230 Y max ■ 1 . 4 

10240 Yateo-.2 

10250 ELSE 

10260 Yrnin-0 

10270 Ymax-15 

10280 Ystep«5 

10290 ENO IF 

10300 ELSE 

10310 Opo-2 

10320 Ymin-0 

10330 Ymax-12 

10340 Yatep-3 

10350 Xmin-0 

10360 Xmax-4 

10370 X s t cp* 1 

10380 ENO IF 

10390 IF Ihrat-1 THEN 

10400 Ym i n«0 

10410 Ynax-18 

10420 Ystep«3 

10430 Xmin«0 

10440 Xmax-9 

10450 Xstep-2 

10460 ENO IF 

10470 BEEP 

10480 PRINT ' IN ; SP 1 s IP 2300.2200,8300,68001* 

10490 PRINT 'SC 0,100 ,0,100;TL 2,0i* 

10500 S f x« 1 00/ < Xmax -Xmi n > 

10510 Sfy-100/( Ynax-Ynin > 

10520 PRINT 'PU 0,0 PO' 

10530 FOR Xa«Xmin TO Xnax STEP Xatep 
10540 X«( Xa-Xmin )*S f x 
1 0550 PRINT 'PA' tX ,0i XT i ' 

10560 NEXT Xa 

10570 PRINT 'PA 100,0«PUi ' 

10580 PRINT 'PU PA 0,0 PO' 

10590 FOR Ya*Ymin TO Ymax STEP Yatep 
10600 Y-< Ya-Ymin >»Sfy 
10610 PRINT 'PA 0,' iY , ' YT ’ 

10620 NEXT Ya 

10630 PRINT 'PA 0,100 TL 0 2" 

10640 FOR Xa«Xmin TO Xnax STEP Xatep 
1 0650 X«( Xa-Xmin >*Sfx 
10660 PRINT 'PA' |X ,100i XT- 
10670 NEXT Xa 

10680 PRINT *PA 100,100 PU PA 100,0 PO* 

10690 FOR Ya-Ymin TO Ymax STEP Yatep 

10700 Y-( Ya-Ymin)*Sfy 

10710 PRINT 'PO PA 1 00 , ' , Y , ' YT ' 

10720 NEXT Ya 
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'CP -2 .-1 iLB * jXai " 
*CP -I .5.-1 « LB ' ; Xa i " ' 



0730 PRINT 'PA 100,100 p U‘ 

0740 PRINT 'PA 0 3ft 1 .5 .2“ 

0750 FOR Xa-Xmn TO Xmax STEP Xatep 
07B0 X-< Xa-Xnin >*Sfx 
0770 PRINT * PA ' i X . * ,0 i ’ 

0780 IF Iuo-0 THEN PRINT 

0790 IF Iuo-1 THEN PRINT 

0800 NEXT Xa 
0810 PRINT "PU PA 0,0' 

0820 FOR Ya-Ymin TO Ywax STEP Yatep 
0830 IF ABS( Ya XI .E-5 THEN Ya-0 
0840 Y«( Ya-Ynin )*Sfy 
0B50 PRINT 'PA 0 ;Y ' 

08S0 IF Iuo-0 THEN PRINT * CP -4 ,-.25; LB’ i Ya » " " 

0870 IF Iuo-I THEN PRINT "CP -3 . - . 25 i LB " ; V a « ' " 

0880 NEXT Ya 

0890 XIabalS-'Oil Percent' 

0900 IF Iuo-0 THEN 

0910 IF Irt-0 THEN 

0920 Y 1 abe 1 $• ' h/h0X “ 

0930 ELSE 

0940 YlabelS-'h/hamooth' 

0950 END IF 

0980 PRINT 'SR 1 . 5 ,2 iPU PA 50,-10 CP" i-LEN( XlabalS ) ✓ 2 ; * 0 i LB * jXlabalS ; ' ' 

0970 PRINT 'PA -11,50 CP 0 ," ;-LEN( Yiabel* >/2*S/6 ; 'DI 0 , 1 <LB' i YlabeiSi' " 

0980 PRINT 'CP 0 ,0' 

0990 ELSE 

1000 PRINT ‘ SP0 i SP2* 

1010 PRINT 'SR 1 .2 ,2 . 4 { PU PA -B,35;0I 0,1iLBU;PR 1,0.5<LBoiPR - 1 , 0 . 5 ; L B (kU/m 

1020 PRINT 'PR -I ,0 . 5 ; SR 1 . 1 . 5 1 LB2 ; SR 1.5,2iPR .5,.5 iLB.iPR .5,0iLBK)' 

1030 PRINT 'PA 42 ,- 1 0 i DI 1,0;LBU;PR .4,-l ; LBw<PR 1 , . 5 1 LB < m/s >* 

1 040 PRINT ' SP0 ; SP 1 ' 

1050 END IF 
1080 Ion-0 
1070 BEEP 

1080 INPUT 'UANT TO PLOT DATA FROM A FILE ( 1 -Y , 0-N >*> " ,Okp 

1090 Icn-0 

1100 IF Okp-1 THEN 

1110 BEEP 

1120 INPUT 'ENTER THE NAME OF THE DATA FILE" ,0_f i let 
1130* IF Iuo-0 THEN 
1 140 BEEP 

1150 INPUT 'SELECT (0-LINEAR, 1 -L0S( X , Y ) " , I Ioq 
1160' END IF 

1 170 ASSI6N ©File TO 0_f ileS 
1 IB0 BEEP 

1190 INPUT 'ENTER THE BEGINNING RUN NUMBER' ,Md 
1200 BEEP 

1210 INPUT 'ENTER THE NUMBER OF X-Y PAIRS STORED' .Npairs 
1220 IF Iuo-0 AND Ihrat-0 THEN 
1230 BEEP 

1240 INPUT 'ENTER OESIRED HEAT FLUX', Q 
1250 END IF 
1280 BEEP 
1270 PRINTER IS 1 
1 2B0 PRINT USIN6 '4X 

1290 PRINT USING 

1300 PRINT USING 

’310 PRINT USIN6 '4X,"5 

1 320 PRINT USIN6 '4X , "6 

1330 PRINT USIN6 '4X ,"7 
1340 INPUT Sy* 

1350 PRINTER IS 705 
1380 PRINT 'PU DI' 

1370 IF Sym-1 THEN PRINT 



' 4X , 

*4X . 



Select a synbo 1 : " " 

1 Star 2 Pius siQn" " 

3 Circle 4 Sauare* ' ' 

5 Ronous"* 

6 Right -s l de-up triangle"' 

7 Up-side-down triangle'"* 



■Sh* 
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1380 IF Sy«*-2 THEN PRINT *SM+* 
1390 IF Syn-3 TnEN rRlNi sflo 
1400 Nn-4 

1410 IF I 1 oo" 1 THEN Nn-1 
1420 IF Md> 1 THEN 
1430 FOR I-l TO (Md-1 > 

1440 ENTER ©FileiXa ,Ya 

1450 NEXT I 

U60 END IF 

1470 IF Ihrat-0 THEN 

1480 QI-Q 

1490 IF I 1 oo- 1 THEN Q-LOS(Q) 
1500 ENO IF 

1510 FOR 1-1 TO Npaira 

1520 IF Iuo-0 AND Ihrat-0 THEN 

1530 ENTER ©FUejXa ,8( ♦ ) 

1540 Ya-B ( 0 ) 

1550 FOR K-l TO Nn 
1560 Y a-Ya+B( K >*Q A K 
1570 NEXT K 
1580 END IF 

1590 IF Iuo-1 OR Ihrat-1 THEN 
1600 ENTER 0F 1 1 a i Xa ,Ya 
1610 IF Iuo-1 THEN Ya-Ya/1000 
1620 ENO IF 

1630 IF Iuo-0 AND Ihrat-0 THEN 

1640 IF 1 1 oo-l THEN Ya-EXP(Ya) 

1650 IF I log-0 THEN Ya-Ql/Ya 

1660 IF Irt-0 THEN 

1670 IF Xa-0 THEN 

1680 Yo-Ya 

1690 Ya-1 

1700 ELSE 

1710 Ya-Ya/Yo 

1720 ENO IF 

1730 ELSE 

1740 Ha w-FNHamoot h( Q ,Xa , I sat > 
1750 Ya-Ya/Han 
1760 END IF 
1770 END IF 
1 780 X x ( 1-1 >-Xa 



1790 
1800 
1810 
1820 
1830 
1840 
1850 
1 860 
1870 
1880 
1890 
1 900 
1910 
1920 
1930 
1 940 
1950 
1 960 
1970 
1 980 
1 990 
2000 
2010 
2020 
2030 



Yy( I-l >-Y a 

X -( Xa-Xn in )*5 f x 

Y-< Ya-Ynin )*Sfy 

IF Syn>3 THEN PRINT "SfT 

IF Syn<4 THEN PRINT *SR 1.4.2.4" 

PRINT *P A * t X ,Y 

IF Syn > 3 THEN PRINT *5R 1. 2,1.6’ 

IF Syn-4 THEN PRINT ’UC2 ,4 , 99, 0,-8, -4, 0,0, 8,4,0, i’ 
IF Syn-5 THEN PRINT • UC3 . 0 , 9 9 , -3 . -6 . -3 , 6 . 3 , 6 .3 . -6 i ' 
IF Sym-G THEN PRINT “UC0 ,5 . 3 ,99 ,3 ,-8 .-6 ,0 ,3 ,8 i ’ 

IF Sym-7 THEN PRINT *UC0 .-5 .3 ,99 ,-3 ,3 ,6 ,0 .-3 ,-8 i " 

NEXT I 
SEE P 

ASSI6N ©File TO • 

ENO IF 

PRINT - PU SfT 
BEEP 

INPUT ’WANT TO PLOT A POLYNOMIAL (1 -Y ,0-N )? * ,Okp 

IF Okp-1 THEN 

BEEP 

PRINTER IS 1 

PRINT USIN6 *4X ‘Select Una type:"" 

PRINT USIN6 "6X.""0 Solid Una""" 

PRINT USIN6 " 6X , " “ 1 Daahed * * * 

PRINT USING *6X , * *2 , , ,5 Longar lma - daa n""" 
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12040 INPUT Ipn 

rnlNTEn la 70a 
12050 BEEP 

12070 INPUT 'SELECT ( 0-LINEAR , 1 -L06< X , Y ) > * , 1 1 og 

12080 Iprnt-1 

12090 CALL Poly 

12100 IF Iuo-I THEN 

12110 BEEP 

12120 INPUT ‘DESIRE TO SET X Lower and Upper Limit ( ) -Y ,0-N )? ‘ . Ix 1 in 

12130 IF Ixlim-0 THEN 

12140 X 1 1-0 

12150 Xu 1 -7 

12150 ENO IF 

12170 IF Ixlim-1 THEN 

12180 BEEP 

12190 INPUT ‘ENTER X Lower Limit* ,X 1 1 
12200 BEEP 

12210 INPUT ‘ENTER X Upper Limit* ,Xui 
12220 ENO IF 

12230 ENO IF 

12240 FOR Xe-XIl TO Xul STEP Xstep/25 
12250 Icn-Icn+1 
12250 Ya-FNPoly(Xe ) 

12270 IF Iuo-1 THEN Ye-Ya/1000 

12280 Y-( Ye-Ynin )«S f y 

12290 X-< Xe-Xwin >*Sfx 

12300 IF Y<0 THEN Y-0 

12310 IF Y> 1 00 THEN GOTO 12410 

12320 Pu-0 

12330 IF Ipn* 1 THEN Idf-Icn MOO 2 
12340 IF Ipn-2 THEN Idf-Icn MOO 4 
12350 IF Ipn-3 THEN Idf-Icn MOO 8 
12360 IF Ipn-4 THEN Idf-Icn MOD 16 

12370 IF Ipn-5 THEN Idf-Icn MOO 32 

12380 IF Idf-1 THEN P u -1 

12390 IF Pu-0 THEN PRINT ‘PA* ,X , Y , *PO * 

12400 IF Pu-1 THEN PRINT " PA ‘ , X , Y , ‘ PU ‘ 

12410 NEXT Xa 
12420 PRINT ‘PU“ 

12430 60T0 1 1070 
12440 ENO IF 
12450 8EEP 

1 2460 INPUT " UANT TO QUIT ( 1 -Y ,0-N ) 7 ‘ , Iqui t 
12470 IF Iouit-I THEN 12490 
12480 GOTO 1 1070 
12490 PRINT ‘PU SP0‘ 

12500 SUBENO 
12510 SUB Steta 
12520 PRINTER IS 701 
12530 J-0 
12540 K-0 
12550 BEEP 

12550 IF Iplot-I THEN ASSIGN »Fila TO P_file$ 

12570 8EEP 

12580 INPUT ‘LAST RUN No 7 ( 0-QUIT ) ' ,Nn 

12590 IF Nn-0 THEN 12950 

12600 Nn-Nn-J 

12610 S x — 0 

12620 Sy-0 

12630 Sz-0 

12640 Sxa-0 

12650 Sy a-0 

12560 Szs-0 

12670 FOR 1-1 TO Nn 

12680 J-J+1 

12690 ENTER iFileiQ.T 
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2~00 h.Q/T 
2710 Sx-Sx+Q 
2720 Sxa-Sxa+Q'2 
2730 Sy-Sy+T 
2740 Sya-Sya+T‘2 
2750 Sz-Sz+H 
2760 Sza-Sza^H-2 
2770 NEXT I 
2780 Qave-Sx/Nn 
2790 Tave-Sy/Nn 
2800 Have-Sz/Nn 

2810 Sdevq-SQR< A8S< (Nn*Sxa-Sx‘2 >/<Nn*(Nn-1 > > ) ) 

2820 Sdevt-SQR< A8S< (Nn*Sya-5y*2 )/(Nn*( Nn-1 )))) 

2830 Sdevh-SQR< A8S< ( Nn*Sza-Sz "2 )/<NnMNn-1 )))) 

2840 Sh- 1 00*Sde vh/Ha ve 
2850 Sq- 1 00*Sdevq/Qave 
2860 St- 1 00*Sdevt/Tave 
2870 IF K - 1 THEN 12930 
2880 PRINT 

2890 PRINT USING 'IIX/'OATA FILE: ' ‘ . 1 4A’ ;Fi let 
2900 PRINT 

2910 PRINT USING * 1 1 X , " “RUN Htube SdevH Qdp SdevO Thetad SdevT*' 

2920 K-1 

2930 PRINT USING " I 1 X ,00 ,2< 2x ,D . 30E , 1 X ,30 . 20 > ,2X ,00 .30 , 1 x ,30 . 20* i J .Have ,Sh .Gave 

Sq .Tave ,St 

2940 GOTO 12570 

2950 ASSIGN ©F i l e 1 TO • 

2960 PRINTER IS 1 
2970 SUBENO 
2980 SU8 Coe f 

2990 COM /Cpiy/ A( 1 0 , 1 0 > ,C< 1 0 > ,B< 5 ) ,N , Iprnt ,0po , 1 1 og , I f n , I j oin ,Njo in 
3000 3EEP 

3010 INPUT *GIUE A NAME FOR CROSS-PLOT FILE’.CpfS 
3020 CREATE 8DAT Cpf$.2 
3030 ASSIGN ©File TO CpfS 
3040 8EEP 

30S0 INPUT 'SELECT ( 0-LlNEAR , 1 -L0S( X ,Y ) >• ,1 log 
3060 8EEP 

3070 INPUT ‘ENTER OIL PERCENT ( - 1 -STOP ) * ,8op 

3080 IF 8op<0 THEN 13120 

3090 CALL Poly 

3100 OUTPUT ©F i l e i 8op ,8< * > 

31 10 GOTO 13060 
3120 ASSIGN ©File TO • 

3’ 30 SU8EN0 

3140 SUB Ui Iaon< Cf ,Ci ) 

3150 COM /Uil/ 02 ,D i ,Do ,L ,Lu ,Kcu 
3160 Oin Emf( 12 ) 

31701 UL I SON PLOT SUBROUTINE OETERMINE CF AND Cl 
3180 BEEP 

3190 INPUT 'ENTER DATA FILE NAME*, File© 

3200 8EEP 

3210 PRINTER IS 1 

3220 PRINT USIN6 "4X. "Select option:*" 

3230 PRINT USING -4X," 0 Uary Cf and Cl"’ 

3240 PRINT USING * 4X , * ' 1 Fix Cf Ua^y Cl’" 

3250 PRINT USIN6 *4X,-- 2 Uary Cf Fix Ci’" 

3260 INPUT ‘ENTER OPTION * , Ic f i* 

3270 PRINTER IS 701 

3280 IF Icfix-0 THEN 13320 

3290 IF lcfix>0 THEN 8EEP 

3300 IF Icfix-I THEN INPUT ‘ENTER Cf.Caf 

3310 IF Icfix-2 THEN INPUT 'ENTER CI*,Ci 

3320 PRINTER IS 1 

3330 INPUT -Want To Uary Coef ?*>< 1 -Y ,0-N >‘ .Iccoef 
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13340 IF Iccoef-1 "HEN INPUT "ENTER COEFF“,R 

13350 PRINTER IS 701 

I33S0 IF Icfix-0 OR Icfix-2 THEN Cfa-.004 

13370 IF Icfix-1 THEN Cfa-Csf 

13300 Cia-Ci 

13390 Xn-.e 

13400 Fr-.3 

13410 Jj-0 

13420 Rr-3 . 

13430 IF Iccoaf-1 THEN Rr-R 

13440 PRINTER IS 1 

13450 PRINT Oo ,0i ,Kcu 

134G0 ASSI6N ©File TO FileS 

13470 ENTER ©File; Nrun , Dates, Ldtcl ,Ldtc2,Itt 

13480 Ru-0o*L06( Oo/Oi )/< 2*Kcu ) 

1 3490 Sx-0 

13500 Sy-0 

13510 Sxy-0 

13520 Sx2-0 

13530 Sy2-0 

13540 FOR 1-1 TO Nrun 

13550 ENTER ©File i0op , Told* ,Emf< •) .Fma 

13560' CONUERT EMF * S TO TEMPERATURE 

13570 FOR J-l TO 5 

13580 T ( J >-FNTvav<Emf< J ) ) 

13590 NEXT J 

13600 Taat ■ ( T( 1 ) + T(2))*.5 
13610 TavQ-T ( 5 ) 

13620 6rad-37.9853+. 104380*Tavg 
13630 Tdrop-Emf<7 )*1 .E+6/( 10.*6rad) 

13640 Tavgc-T<5)-Tdrop*.5 

13650 IF A8S< Tavg-Tavgc )> . 01 THEN 

13660 Tavg- ( Tavg+Tavgc ) • . 5 

13670 60T0 13620 

13680 ENO IF 

136901 

137001 Compute properties of water 
13710 Kw-FNKw(Tavg) 

13720 Muwa-FNMuw( Tavg ) 

13730 Cpw-FNCpw(Tavg) 

13740 Prw-FNPrw( Ta vg ) 

13750 Rhow-FNRhow ( Ta vg ) 

13760' 

137701 Compute properties of Freon-114 
1 3780 Lmt d-Tdrop/L06< ( T< 5 )-Taat >/<T(5)-Tdrop-Tsat > ) 

13790 IF Jj-0 THEN 

13800 Tu-Taat+Fr*Lmt d 

13810 Thetab-Tw-Tsat 

13820 Jj-I 

13830 END IF 

13840 T f-{ Tw + T sat >♦ .5 

13850 Rho-FNRho(Tf ) 

13860 Mu-FNMu<Tf ) 

I 3870 K-FNK( Tf ) 

13880 Cp-FNCp(Tf ) 

13890 Set a-FNBet a( T f ) 

13900 Hf g-FNHf g< T sat ) 

13910 Nt-Mu/Rho 
13920 Alpha-K/(Rho*Cp ) 

13930 Pr-Ni/Alpha 
1 3940 1 

1 3 950 1 Analysis 
139601 COMPUTE MOOT 
13970 A-PI*<0o*2-0i*2)/4 
13980 P-PI *0o 

13990 Mdot *3 . 9657E-3+Fms * (3.61 955E-3-Fma ♦ ( 8 . 82006E-6-Fma • ( 1 . 23688E-7-Fms * 4 . 3 1 897 
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E-10 > > > 

14000 Q-Mdot ♦Cow* Torop 

140101 CDMPUTE NATURAL-CDNUECTI UE HEAT-TRANSFER COEFFICIENT 
14020' FDR UNENHANCED END( S ) 

14030 Hbar-190 

14040 Fe*(Hbar*P/(Kcu*A) )* . 5*Lu 
14050 Tann-FNTanh< Fe ) 

14060 Theta-Thetab*Tanh/Fe 

14070 Xx-(9.81*9eta»Thetab*Do‘3*Tanh/(Fe*Ni*Alona>)'‘. 166567 
14080 Yy-( 1 ♦( . 559/Pr ) “ ( 9/ 1 6 ) > “ ( 8/ 27 > 

14090 Hbarc-K/Do*( ,5+ . 387*Xx/Yy )*2 

14100 IF A8S((Hbar-Hbarc)/Hbar>>.001 THEN 

14110 Hbar-( Hbar+Hbarc )• . 5 

14120 GOTO 14040 

14130 END IF 

1 4 U0 i 

14150' COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 

14160 Q1-(Hbar*P*Kcu»A)" .5*Thetao*Tanh 

14170 Qc-Q-2«Ql 

14180 As-PI •D2*L 

14190' COMPUTE ACTUAL HEAT FLUX 

14200 Odo*0c/As 

14210 IF Icfix-0 OR Icf ix> 1 THEN Csf-1/Cr<1 ./Rr ) 

14220 Thetab-C9f/Co*Hfg*(Qdo/(Mu*HfQ )♦( . 0 1 4/ ( 9 . 8 1 *Rho > > A . 5 ) A ( 1 /Rr ) *Pr * 1 . 7 

14230 Ho-Qdo/Thet ab 

14240 Onega-Ho/Cf 

14250 Uo-Q/(PI*Do«L*Lntd> 

14260 Uw-Mdot / ( Rhou*P I *Di ‘ Z/4 ) 

14270 Rew"Rhow*Uw*Di/Muwa 
14280 Twi-Tw*Q*Rw/ ( P I ♦Do*L > 

14290 Gana-Kw/Di«Reu“.8*Prw'‘<1/3. >♦< Muua/FNMuw( Tw i ) ) * . I 4 

14300' PRINTER IS 1 

14310 Yw-( 1 . /Uo-Rw >*Omega 

14320 Xu-0nega*0o/ ( Sana *0i ) 

14330 Sx -Sx + Xu 
1 4340 Sy-Sy+Yu 
14350 Sxy-Sxy ♦Yu*Xw 
14360 Sx2-5x2+Xw*Xy 
14370 Sy2-Sy2+Yu* Yu 
14380 NEXT I 
14390 ASSIGN ©File TO ♦ 

1 4400 M"(Sx«5y“Nrun*Sxy )/( Sx*Sx-Nrun*Sx2 ) 

14410 C-( Sy-Sx*M )/Nrun 

14420 IF Icflx-0 OR Icfix-3 OR Icflx-4 THEN 

14430 Cic-I/M 

14440 Cfc-I/C 

14450 END IF 

14460 IF Icf ix- 1 THEN 

14470 Cic-1/M 

14480 Cfc-Cf 

14490 END IF 

14500 IF Icfix-2 THEN 

14510 Cic-Cl 

14520 Cfc-l/C 

14530 END IF 

14540 IF ABS< ( Ci-Cic )/Cic >> .001 DR ABS( (Cf-Cfc >/Cfc >>.001 THEN 
14550 Ci-( Cl+Clc > ♦ . 5 
14560 Cf •( Cf +Cfc >♦ .5 
14570 PRINTER IS 1 

14580 PRINT USING "2X , * * Csf - * * ,MZ . 3DE ,2X , * * Ci - “ “ ,M2 . 3DE* i Cs f ,Ci 

14590 PRINTER IS 701 

14600 6DT0 13460 

14610 END IF 

14620 PRINT 

14630 PRINTER IS 701 

14640 PRINT USING Cf Ci**" 
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U65C PRINT USING " 8X , “ " ASSUMED , MZ , 2DE > 3X , MZ . 3DE" ;C c a , 2i a 

14663 PRINT USIN6 " 8X / " CALCULATED ,MZ . 3DE ,3X ,MZ . 3DE " 5 Ca f , 2 ; 

14670 PRINT 

146 80 Sum2-Sy2-2«M*Sxy-2*C*Sy+M'‘2*Sx2*2 t M*C«Sx + Nrun»C ,, ‘2 
14690 PRINT USIN6 *10X,*-Sun of Squares - " ’ , Z . 3DE * ; Sun2 
14700 PRINT USIN6 " 1 0X , * • Coe f f i c 1 ant - ’ * .0 . DDD" 1 Rr 
14710 SUBEND 
14720 DEF FNMuu(T) 

14730 A-247.8/(T+l33. 15 ) 

14740 Mu* 2 . 4E-5* 1 0 * A 
14750 RETURN Mu 
14760 FNEND 
14770 DEF FNCpu( T > 

14780 Cpw-4.21 120858-T*(2.26826E- 3- T*(4. 42361 E- 5+2. 71428 E-7*T)) 

14790 RETURN Cpu*1000 

14800 FNEND 

14810 DEF FNRhow(T) 

14820 Ro-999. 52946 + T*( . 0 1 26 9-T* ( 5 . 4825 1 3E-3-T * 1 .234 147E-5) > 

14830 RETURN Ro 
14840 FNEND 
14850 DEF FNPrw(T ) 

14 860 Pr u-FNCpw( T ) # FNHuw( T )/FNKu( T ) 

14870 RETURN Pru 
14880 FNEND 
1 4890 DEF FNKu(T ) 

14900 X-( T + 273 . 15 > /273 . 1 5 

14910 Kw— .92247+X*(2.8395-X*( 1 ,8007-X*< . S2577- . 07344 »X )) > 

14920 RETURN Ku 
14930 FNEND 
14940 SUB Plot 

14950 COM /Cply/ A( 10,10 ) ,C( 10) ,8(5) ,Nop , Iornt ,0po .Ilog.Ifn.Ijoin.Njoin 
14960 DIM Ba(3) 

14970 INTE6ER Ii 
14980 PRINTER IS 1 



14990 


Idv-0 








15000 


8EEP 








15010 


input 


“LIKE 


DEFAULT 


UALUES FOR PLOT ( 1 - 


15020 


Ooo-0 








15030 


8EEP 








1 5040 


PRINT 


USINS 


’4X, ""Select Option:""’ 


15050 


PRINT 


USING 


"6X .- -0 


q versus Pelta-T"" 


15060 


PRINT 


USIN6 


* 6 X , " " 1 


h versus delta-T"" 


15070 


PRINT 


USING 


* 6X , ' ’ 2 


h versus q’*“ 


15080 


INPUT 


Opo 






15090 


8EEP 








15100 


INPUT 


•SELECT UNITS 


(0-SI , 1 -ENGLISH > “ ,1. 



15110 PRINTER IS 705 
15120 IF IdvOI THEN 
15130 BEEP 

15140 INPUT “ENTER NUMBER OF CYCLES FOR X-AXIS", Cx 
15150 8EEP 

15160 INPUT “ENTER NUMBER OF CYCLES FOR Y-AXIS’, Cy 
15170 8EEP 

15180 INPUT -ENTER MIN X-UALUE (MULTIPLE OF 10)“,Xnm 
15190 8EEP 

15200 INPUT ‘ENTER MIN Y-UALUE (MULTIPLE OF 10)".Ynin 

15210 ELSE 

15220 IF Opo-0 THEN 

15230 Cy-3 

15240 Cx-3 

1 5250 Xni n- . 1 

15260 Yni n- 1 00 

15270 END IF 

15280 IF Opo-1 THEN 

15290 Cy-2 

15300 Cx-2 
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153’ 0 


X m ; n- . 1 




15320 


r n- 1 00 




15330 


END IF 




1 5340 


IF Opo-2 THEN 




15350 


IF Iun-0 THEN 




153B0 


Cy-3 




15370 


Cx-2 




15380 


Xn m- 1 000 




15390 


Ypu n- 1 00 




15400 


ELSE 




15410 


Cy-3 




15420 


Cx-3 




15430 


Xnin-1 00 




15440 


Yrun- 1 0 




1 5450 


END IF 




15460 


END IF 




15470 


END IF 




15480 


BEEP 




15490 


PRINT " IN s SP 1 s IP 2300 ,2200 ,8300 ,6800i" 


15500 


PRINT "SC 0 , 1 00 , 0 , 1 00 i TL 2,0s" 


15510 


Sfx-100/Cx 




15520 


Sfy-100/Cy 




15530 


SEEP 




1 5540 


INPUT 'WANT TO BY -PASS 


CASE 7 ( 1-Y , 0-N ) ’ , Ibyp 


15550 


IF Ibyp-I THEN 16790 




15560 


PRINT "PU 0,0 PD" 




15570 


Nn-9 




15580 


FOR I-l TO Cx ♦ 1 




15590 


Xat-Xmn* 1 0~( I- 1 ) 




15600 


IF I-Cx+1 THEN Nn-I 




15610 


FOR J-1 TO Nn 




15620 


IF J-1 THEN PRINT " TL 2 


0" 


15630 


IF J-2 THEN PRINT "TL 1 


0" 


15640 


Xa-Xat * J 




15650 


X-L6T(Xa/Xrun )•$ f x 




15660 


PRINT ’PA * i X ,0; X T i * 




15670 


NEXT J 




15680 


NEXT I 




15690 


PRINT ’PA 100 .01 PU i * 




15700 


PRINT "PU PA 0 ,0 PD" 




15710 


Nn-9 




15720 


FOR 1-1 TO Cy ♦ 1 




15730 


Yat-Ynin* 1 0 * ( I- 1 ) 




15^40 


IF I-Cy+1 THEN Nn- 1 




15750 


FOR J-1 TO Nn 




15760 


IF J-1 THEN PRINT "TL 2 


0" 


15770 


IF J-2 THEN PRINT "TL 1 


0" 


15780 


Ya-Yat • J 




15790 


Y-L6T( Ya/Ynin )*Sfy 




15800 


PRINT "PA 0 iY , "YT" 




15810 


NEXT J 




15820 


NEXT I 




15830 


PRINT "PA 0 ,100 TL 0 2" 




15840 


Nn-9 




15850 


FOP 1-1 TO Cx+1 




15860 


Xat-Xmn* 1 0 A < 1-1 ) 




15870 


IF I-Cx+1 THEN Nn- 1 




15880 


FOR J-1 TO Nn 




15890 


IF J-1 THEN PRINT "TL 0 


2" 


15900 


IF J>1 THEN PRINT "TL 0 


1 " 


15910 


Xa-Xat * J 




15920 


X-t6T( Xa/Xntn )*Sf x 




15930 


PRINT "PA * i X , " , 1001 XT- 




15940 


NEXT J 




15950 


NEXT I 




15960 


PRINT "PA 100 ,100 PU PA 


100 ,0 PD" 
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15970 Nn-9 

5390 P QR 1-1 TO Cyi-1 
15990 Vat -Y«m* 1 0* ( 1-1 ) 

16000 IF I-Cy+1 THEN Nn- 1 

16010 FOR J-l TO Nn 

16020 IF J-1 THEN PRINT *TL 0 2’ 

16030 IF J>1 THEN PRINT *TL 0 1" 

16040 Ya-Vat*J 

16050 Y-L6T ( Ya/Ywin )*Sf y 

16060 PRINT ’PO PA 100,' , Y , ’ YT ’ 

16070 NEXT J 
16080 NEXT I 

16090 PRINT *PA 100,100 PU’ 

16100 PRINT ‘PA 0 ,-2 SR 1 .5,2’ 

16110 Ii-LGT(Xrun ) 

16120 FDR 1-1 TO Cx+1 
16130 Xa-Xml n« 1 0 A ( 1-1 ) 

16U0 X-LST(Xa/Xnin)*Sfx 
16150 PRINT *PA ’ i X .’ ,0* * 

16160 IF It>-0 THEN PRINT * CP -2 , -2 i LB 1 0 \ PR -2 .2 \ LB " ; 1 1 ; * * 

16170 IF 1 1 <0 THEN PRINT ’CP -2 , -2 s LB 1 0 i PR 0.2iLB";Iu"" 

16180 Ii-Ii+I 

16190 NEXT I 

16200 PRINT ’PU PA 0 ,0’ 

16210 I i-LST ( Ynin > 

16220 Y 1 0-1 0 

16230 FOR 1-1 TO Cy+1 

16240 Ya-Ymin* 1 0* ( 1-1 > 

16250 Y-[_ST< Ya/Ynin >*$fy 
16260 PRINT ’PA 0 ,* i Y 

16270 PRINT ’CP -4 , - . 25 i LB 1 0 i PR -2 ,2 i LB * 1 1 1 j ‘ * 

16280 Ii-Ii+I 
16290 NEXT I 
16300 BEEP 

16310 INPUT ’WANT USE DEFAULT LABELS ( 1 -Y ,0-N >? ’ , I di 
16320 IF IdlOl THEN 
16330 BEEP 

16340 INPUT ’ENTER X-LABEL ’ ,X iaoelS 
16350 BEEP 

16360 INPUT ’ENTER Y-LABEL ’ ,Y label* 

16370 END IF 

16380 IF 0po<2 THEN 

16390 PRINT ’SR 1 ,2 1 PU PA 40,-Ui’ 

16400 PRINT ’LB( T j PR -1.6,3 PD PR 1.2,0 PUiPR . 5 , -4 j LBwo \ PR .5,1;’ 

16410 PRINT ’LB-TiPR .5,-1 ; LBsat i PR .5.1 s’ 

16420 IF Iun-0 THEN 
16430 PRINT ’LB > (K)’ 

16440 ELSE 

16450 PRINT ’LB) <F>’ 

16460 END IF 
16470 END IF 
16480 IF Dpo-2 THEN 
16490 IF Iun-0 THEN 

16500 PRINT ’SR 1 . 5 ,2 i PU PA 40,-MiLBq (U/wjSR I.I.SiPR 0.5,MLB2iSR 1 .5 ,2 *PR 
0.5.-1 i LB ) ’ 

16510 ELSE 

16520 PRINT ’SR 1 .5 ,2 iPU PA 34,-l4;LBq (Btu/hnPR .S,.SiLB.iPR .5,-. Si" 

16530 PRINT ’LBftiPR .5,1 i SR 1,1.5 iLB 2;SR 1.5,2iPR .5,-1 i LB > i ' 

16540 END IF 
16550 END IF 
16560 IF Dpo-0 THEN 
16570 IF Iun-0 THEN 

16580 PRINT ’SR 1 .5 ,2 iPU PA -12,40iDI 0,1iLBq (U/niPR -1,0.5iSR 1 . 1 . 5 t LB2 i SR 1 
. 5 ,2 i PR 1 , .5;LB )“ 

16590 ELSE 

16600 PRINT ’SR 1 . 5 ,2 i PU PA -1 2 .32 iDI 0,1;LBq (Btu/hriPR - . 5 , . 5 1 LB . i PR .5 , . 5 i * 
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’ 6S ' 0 PRINT "l6't;SR 1 .5« p R -I , .5i^B2;PP . 5i5R : .5,2;lB * 

16620 END IF 
15630 END IF 
16640 IF Ooo>0 THEN 
16650 IF Iun-0 THEN 

1 6650 1 PRINT * SR 1 .5 .2 i PU PA -12 ,3BiDI 0 , 1 « LBh (U/mPR -I ,.5i SR 1 .1 .SilB 2;SR I 
. 5 , 2 i PR .5, .Si- 

16670 PRINT 'SR I .2 ,2 . 4 i PU PA -12 .37 * D I 0,liLBhiPR 1,0.5iLBoiPR -l,0.5«LB (U/« 



16680 PRINT ‘PR - 1 , . 5 i SR 1,1.5iLB2iSR 1 .5 .2 i PR .5,.5;LB.iPR .5 f 0iL8K>" 

16690 ELSE 

16700 PRINT * SR 1 . 5 ,2 ; PU PA -12 .28 i DI 0,l»LBh <8tu/hriPR - . 5 . . 5 1 LB . i PR .5, .Si- 

16710 PRINT 'LBftiPR -1 ,.5 i SR l ,1 . 5 i L82 i SR l.5,2«PR .5,.5«L8.|PR .5,.5iLSF>" 

16720 END IF 

16730 END IF 

16740 IF Idl-0 THEN 

16750 PRINT ‘SR 1 . 5 ,2 i PU PA 50.-16 CP* i-LEN( XlabelS >/2; -0;L8- iXlabelSi * " 

16760 PRINT ‘PA -14,50 CP 0 , * i -LEN( YlabelS )/2*5/6i "DI 0,1 i L8 * ; Y label $ i " " 

16770 PRINT - CP 0.0 DI" 

16780 END IF 
16790 Ipn-0 
16800 Xll-1 . E+6 
16B10 Xul— 1 . E+6 
16B20 Icn-0 
16830 Ifn-0 



16840 
16850 
1 6860 
16870 
16880 
16890 
16900 
16910 
16920 
16930 
16940 
16950 
16960 
16970 
16980 
16990 
17000 
17010 
17020 
17030 
1 7040 
17050 
17060 
17070 
17080 
17090 
17100 



I J oin- 1 

8EEP 

INPUT ’WANT TD PLDT DATA FROM A FILE ( 1 - Y ,0-N P * ,Dk 

IF Dk-1 THEN 

BEEP 

INPUT ’ENTER THE NAME DF THE DATA FILE * ,D_f l 1 eS 

ASSIGN ©File TD D_fiieS 

8EEP 

BEEP 

INPUT 'ENTER THE BEGINNING RUN NUMBER * ,Md 
8EEP 

INPUT 'ENTER THE NUMBER DF X-Y PAIRS STDRED* .Npairs 
BEEP 

INPUT "CONNECT DATA WITH LINE ( 1 -Y ,0-N P ' , Ic 1 
8EEP 

PRINTER IS 1 

PRINT USING * 4X .* "Select 



PRINT USIN6 "6X , • " l 



PRINT USING * 6X , " * 3 
PRINT USIN6 * 6X , ■ " 5 
PRINT USING * 6X , * "6 
D R INT USIN6 ’ 6X -7 
INPUT Syr 
PRINTER IS 705 
PRINT -PU DI" 

IF Syn-I THEN PRINT "SM«" 
IF Syr-2 THEN PRINT -SM*" 



synool : " * 

Star 2 Plus 9 ign" " " 

Circle 4 Square"'* 

Ronbus " ’ " 

Ri gnt -3 l de-uo triangle**' 
Ud- 9 l de-doun triangle**' 



17110 IF Syr-3 THEN PRINT *SMo" 
17120 IF Md> 1 THEN 
17130 FDR I-l TO (Md-1 ) 



17140 ENTER ©FileiYa.Xa 

17150 NEXT I 

17160 END IF 

17170 FDR I-l TD Npairs 

17180 ENTER ©FileiYa.Xa 

17190 IF 1-1 THEN Qt-Ya 



17200 IF I-Nqairg THEN 32-Ya 
17210 IF Doo- 1 THEN Ya-Ya/Xa 
17220 IF Doo-2 THEN 



17230 O-Ya 
17240 Ya-Ya/Xa 
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r:se «a-o 

17260 £N0 IF 

17270 IF Xa<X 1 1 THEN Xll-Xa 
17280 IF Xa>Xul THEN Xul-Xa 
17290 IF Iun-1 THEN 

17300 IF 0 O o<2 THEN Xa-Xa* 1.8 

17310 IF Op o > 0 THEN Ya-Ya*.1761 

17320 IF Odo-0 THEN Ya-Ya*.3l7 

17330 IF Qpo-2 THEN Xa-Xa*. 317 

17340 ENO IF 



17350 X-L6T(Xa/Xmin)*Sfx 
17360 Y-L6T( Ya/Ymln )*Sfy 
17370 Kj-0 

17380 CALL Synb(X,Y,Sym,Icl ,Kj ) 

17390 60T0 17520 

17400 IF Syn>3 THEN PRINT ‘ SPT 

17410 IF Syn<4 THEN PRINT 'SR 1.4, 2. 4' 

17420 IF Icl-0 THEN 
17430 PRINT ' PA ' ,X ,Y ' 

17440 ELSE 

17450 PRINT ' P A ‘ ,X , Y , ' PO “ 

17460 ENO IF 

17470 IF Syw>3 THEN PRINT 'SR 1.2. 1.6“ 

17480 IF Syn-4 THEN PRINT “UC2 ,4,99 ,0 ,-8 ,-4 ,0,0,8 ,4 ,0» “ 
17490 IF Syn-5 THEN PRINT “ UC3 ,0 , 99 , -3 .-6 , -3 ,6 , 3 , 6 , 3 .-6 »' 
17500 IF Syn-6 THEN PRINT ' UC0 ,5 . 3 , 99 ,3 ,-8 ,-6 ,0 , 3 , 8 i “ 
17510 IF Sym-7 THEN PRINT “UC0 ,-S . 3 .99 ,-3 ,8 .6 ,0 .-3 ,-8i * 
17520 NEXT I 
17530 PRINT ‘PU“ 

17540 8EEP 

17550 INPUT “WANT TO LABEL 7 < 1 -Y ,0-N > “ , I lab 
17560 IF I lab- 1 THEN 
17570 PRINT “ SP0 < SP2 " 

17580 8EEP 

17530 IF K lab-0 THEN 
17600 X lab-5 
17610 Y lab-85 

17620 INPUT “ENTER INITIAL X ,Y LOCATI ONS “ ,X lab , Y 1 ab 
17630 Xt t-X lab-5 
17640 Ytt-Ylab+8 



17650 PRINT “SR 1,1 .5“ 

17660 PRINT *SM«PA“ ,Xtt ,Ytt ,'LB X Heat File- 

17670 Ytt-Ytt-3 

17680 PRINT “PA’ ,Xtt ,Ytt ,“LB Oil Flux Nane“ 

17690 IF Sy«-1 THEN PRINT ’S 
17700 IF Sym-2 THEN PRINT “SF1 + “ 

17710 IF Sym-3 THEN PRINT “SMo- 
17720 K 1 ab- 1 
17730 ENO IF 



17740 Kj-1 

17750 CALL Symb(Xlab ,Ylab .Syn.Icl ,Kj ) 

17760 PRINT “SR 1,1 .5 < SM* 

17770 IF Syn<4 THEN PRINT “PR 2,0“ 

17780 BEEP 

17790 INPUT “ENTER BOP“ ,Bop 

1 7800 IF Bop< 1 0 THEN PRINT “PR 2 , 0 i LB “ i Bop i“ * 
17810 IF Bop>9 THEN PRINT “PR . 5 , 0 1 LB “ « 8op i “ “ 
17820 Ihf-0 

17830 IF Q 1 >02 THEN Ihf-1 

17840 IF Ihf-0 THEN PRINT “PR 4,0iLBInc“ 

17850 IF Ihf-1 THEN PRINT “PR 4,0,LBOec' 

17860 PRINT “PR 2 , 0 i LB “ i 0_f i le* i “ “ 

17870 PRINT ’ SP0 1 SP 1 i SR 1 . S ,2 “ 

17880 Ylab-Ylab-5 
17890 ENO IF 
17900 BEEP 
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17910 ASSIGN TO * 

17920 Aii-xil/1 .2 



'"Select line type : 

”0 Solid line'** 

"l Dashed’” 

’’2,, ,5 Longer line - dash"'* 



17930 Xul-Xul»l .2 
17940' SOTO 8040 
17950 END IF 
17960 PRINT *PU SM’ 

17970 8EEP 

1 ^990 INPUT * WANT TD PLOT A POLYNOMIAL ( 1 -Y ,0-N )? ’ ,6o_on 
17990 IF 5o_on-J THEN 
18000 8EEP 
18010 PRINTER IS 1 
18020 PRINT USIN6 *4X , 

18030 PRINT USIN6 *6X , 

18040 PRINT USIN6 *8X , 

18050 PRINT USIN6 *6X , 

18060 INPUT Ipn 
19070 PRINTER IS 705 
18080 8EEP 

18090 INPUT ’SELECT ( 0-LIN , 1 -L06( X ,Y ) )* .1 Ioq 
18100 Iprnt-1 
18110 CALL Poly 
18120 IF Ifn-0 THEN 
18130 8EEP 

18140 INPUT ’ENTER NUM8ER OF FILES TO JOIN'?’ .Njoin 
18150 ENO IF 
18160 Ijoin-0 

18170 IF Ifn<Njoln THEN Ijoin-1 
18180 IF Ifn-0 OR Ijoin-1 THEN 
18190 FOR Ij-0 TO 3 
18200 8 a ( Ij )-8 a ( I j ) + 8( I j ) 

18210 NEXT Ij 
18220 Ifn-Ifn+1 
18230 END IF 
18240 IF Njoin-Ifn THEN 
18250 FOR Ij-0 TD 3 
18260 8 ( I j )-Bs( I j )/Njoin 
18270 Bs( I j )-0 
18280 NEXT Ij 
18290 Ifn-0 
18300 ElSE 
18310 6DT0 16850 
18320 ENO IF 
18330 8EEP 

18340 INPUT ’ENTER Y LOUER AND UPPER LIMITS’ ,Yl 1 ,Yul 
18350 FOR Xx-0 TO Cx STEP Cx/200 
18360 Xa-Xnin» 1 0’Xx 

18370 IF Xa<X 1 1 DR Xa>Xul THEN 18640 
18380 Icn-Icn+1 
18390 p u-0 

18400 IF Ion-1 THEN Idf-Icn MOD 2 

18410 IF Ion-2 THEN Idf-Icn MOO 4 

18420 IF Ion-3 THEN Idf-Icn MOO 8 

18430 IF Ipn-4 THEN Idf-Icn MOO 16 

18440 IF Ipn-5 THEN Idf-Icn MOD 28 

18450 IF Idf-I THEN Pu-1 

18460 IF Doo-0 THEN Ya-FNPo ly < Xa ) 

18470 IF Dpo-2 ANO I1 oq- 0 THEN Ya-Xa/FNPo 1 y ( Xa ) 

18480 IF Doo-2 ANO Ilog-1 THEN Ya-FNPo iy( Xa ) 

18490 IF Opo-1 THEN Ya-FNPo ly< Xa ) 

18500 IF Ya<Ymn THEN 18640 
18510 IF Ya < Y 1 1 OR Ya>Yul THEN 18640 
18520 IF Iun-1 THEN 

18530 IF 0po<2 THEN Xa-Xa«1.8 

18540 IF Oqo>0 THEN Ya- Ya • . 1761 

18550 IF Dpo-0 THEN Ya-Ya.,317 

18560 IF Ooo-2 THEN Xa-Xa».3l7 
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16570 EnD IF 

18530 *-t-87 (. Ta/ inun /*5r y 

18590 X-L5T(Xa/Xnin)*Sfx 

18600 IF Y<3 THEN Y-0 

18610 IF Y > 1 00 THEN 60T0 18640 

18620 IF Pu-0 THEN PRINT ‘PA* ,X , Y , ‘PO* 

18630 IF Pu-1 THEN PRINT ‘PA’ ,X ,Y , “PU“ 

18640 NEXT Xx 
18650 PRINT ‘PU‘ 

18660 60T0 16850 
18670 ENO IF 
18680 8EEP 

18690 INPUT ‘ WANT TO PLOT REILLY’S OATA? ( 1 -Y ,0-N >‘ , Ir ly 

18700 IF Opo-0 OR Opo-1 THEN 

18710 X 1 1-3 

18720 Xul-20 

18730 ENO IF 

18740 IF Opo-2 THEN 

18750 X 1 1-1 0000 

18760 Xul-100000 

18770 ENO IF 

18780 IF Irly-1 THEN 

18790 Yll-20 

18800 Yu 1 *70 

18810 8EEP 

18820 INPUT ‘ENTER LOWER ANO UPPER Y-LIMITS FOR PL0TTIN6" ,Y1I ,Yul 
18830 FOR Xx-0 TO Cx STEP Cx/200 
18840 Xa-Xnin* 10 A Xx 

18850 IF Xa<X 1 1 OR Xa>Xul THEN 18980 
18860 Xl-L06(Xa> 

18870 IF Opo-0 THEN Y 1--2 . 5402837E- 1 +X 1 • < 4 . 9720 1 5 1 -X 1 »2 . 5 1 34787E-1 ) 
18880 IF Opo-1 THEN Y 1--2 . 5402837E- 1 + X 1 • ( 3 . 9720 1 5 1 -X 1 «2 . 5 I 34787E- 1 > 
18890 IF Opo-2 THEN Y 1--3 . 707380 1 E- 1 *X I •( 8 . 72591 90E- 1 -X 1 *6 . 8826842E-3 ) 
189C0 Ya-EXP( Y 1 ) 

18910 Y«L6T< Ya/Ynin.>*Sfy 
18920 X-L6T( Xa/Xnin )*Sf x 
18930 Ipu-0 

18940 IF Y< Y 1 1 THEN Ipu-1 

18950 IF Y > Yu 1 THEN 60T0 18980 

18960 IF Iou-0 THEN PRINT "PA * ,X , Y . ‘ PO ’ 

18970 IF Ipu-1 THEN PRINT " PA " .X , Y , " PU " 

18980 NEXT Xx 
18990 PRINT *PU" 

19000 ENO IF 
19010 8EEP 

19020 INPUT "WANT TO PLOT ROHSENOW CORRELATION? ( 1 - Y ,0-N ) " , Iroha 

19030 IF Iroha- 1 THEN 

19040 YL1-15 

19050 Yu 1 -80 

19060 8EEP 

19070 INPUT ‘ENTER Taat (OaQ C)",Tsat 
19080 Cs f- . 0040 
19090 8EEP 

19100 INPUT ‘ENTER Caf ( OEF-0 . 004 ) ‘ ,Ca f 
19110 Tf-Tsat+2 

19120 FOR Xx-0 TO Cx STEP Cx/200 
19130 Xa-Xmin* 1 0 A Xx 

19140 IF Xa<X 1 1 OR Xa>Xul THEN 19420 
19150 X 1 -L06( Xa ) 

19160 IF Opo<2 THEN Tf -Taat +Xa/2 
19170 Rho-FNRho(Tf ) 

19180 K -FNK ( T f ) 

19190 fiu-FNHu ( T f ) 

19200 Cp-FNCp(Tf ) 

19210 Hf Q-FNHf q( T aat > 

19220 Nl-Hu/Rho 
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9230 -Co * Mu 

'92*i0 Onega*C s f g/Cs • . 0 1 ** 3.5i*Rnc' .9 

19250 IF Opo-0 THEN Ya-< Xa/Onega >"3 

19260 IF Opo-1 THEN Ya-< Xa/Onega ) *3/Xa 

19270 IF Opo-2 THEN Ya-Xa * ( 2 . /3 ) /Onega 

19280 IF Opo-2 THEN 

1 9290 Tfc-Tsat+Xa/Ya* ,5 

19300 IF A8S< T'-Tfc )>.0I THEN 

19310 T*-<Tf+Tfc >*.5 

19320 60T0 19170 

19330 ENO IF 

19340 ENO IF 

19350 Y«L6T( Ya/Ymin >*Sf y 
19360 X-L6T( Xa/Xnin )*Sf x 
19370 Ipu-0 

19380 IF Y < Y 1 1 THEN Ipu-1 
1 9390 IF Y> Yul THEN 1 9420 
19400 IF Ipu-0 THEN DRINT ■ PA • ,X , Y , * PO * 

19410 IF Ipu-1 THEN PRINT * PA * ,X , Y , “ PU " 

19420 NEXT Xx 
19430 PRINT 'PIT 
19440 END IF 
19450 BEEP 

19460 INPUT ’WANT TO QUIT < 1 -Y , 0-N ) * , Ip t 
19470 IF Iqt-1 THEN 19490 
19480 60T0 16800 
1 9490 PRINT *PU PA 0 ,0 SP0 ’ 

19500 SUBENO 

19510 SU8 Sy*D< X ,Y ,Syn . I c 1 ,Kj ) 

19520 IF Sym>3 THEN PRINT *SM* 

19530 IF Sym<4 THEN PRINT * SR 1. 4.2.4* 

19540 Yad-0 

19550 IF Kj-1 THEN Yad-.8 

19550 IF Icl-0 THEN 

19570 PRINT * P A * ,X , Y+Yad , * * 

19580 ELSE 

19590 PRINT *PA* ,X ( Y + Yad , *PO* 

19600 ENO IF 

19610 IF Sy*>3 THEN PRINT "SR 1.2. 1.6* 

19620 IF 5yn-4 THEN PRINT *002,4,99,0,-8.-4.0,0,8,4,05* 

19630 IF Syi^-S THEN PRINT *UC3 ,0 . 99 ,-3 ,-6 , -3 ,6 .3 ,6 ,3 ,-6 ; * 

19640 IF Syn-6 THEN PRINT *UC0 .5 . 3 . 99 ,3 ,-3 ,-6 .0 .3 ,8 i * 

19650 IF Syn-7 THEN PRINT * UC0 .-5 . 3 . 99 , -3 . 8 ,6 ,0 . -3 . -8 ! * 

19660 IF Kj-1 THEN PRINT ’SMiPR 0,-,8* 

19670 SUBEND 
19680 SUB Purg 
196 90 BEEP 

19700 INPUT ‘ENTER FILE NAME TO 8E OELETEO* .File* 

19710 PUR6E File* 

19720 GOTO 19690 
19730 SUBENO 
19740 SUB Tdcn 
19750 COM /Cc/ 0(7) , lea 1 
1 9760 OIM Enf < 1 ) 

19770 DATA 0.10086091,25727.94369,-767345.8295.78025595.81 
19780 OATA -9247486589, 6. 97688E*11 , -2. 56I92E+13. 3. 94078EM4 
19790 READ C< * ) 

19800 BEEP 

19B10 INPUT *6IVE A NAME FOR FILE TO 8E ORE ATEO * ,F 1 1 eS 
19820 BEEP 

19830 INPUT ’SELECT TUBE ( 0-WH , 1 -HF ,2-UT ) * , 1 1 X 
19840 8EEP 

19850 INPUT 'SELECT THERMOCOUPLE TYPE ( 0-NEU , 1 -OLO ) * . leal 
19860 IF 1 1 1 <2 THEN Oi-.0127 
19B70 CREATE BOAT Filet ,4 
19880 ASSI6N SFile TO File* 
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9890 OUTFIT : ie ; 1 1 1 

1S900 J-0 
19910 8EEP 

19920 INPUT ‘ENTER MONTH, OATE ANO TIME ( MM : 00 : HH : MM : SS > * .Oat eS 
19930 OUTPUT 709 i ‘TO* iDataS 
19940 OUTPUT 7 0 9 1 ’TO* 

19950 ENTER 709iOateS 
19950 PRINTER IS 1 
19970 PRINT 

19980 PRINT * Month, date and tine: *iOateS 

19990 PRINT 

20000 PRINT USIN6 *10XT* Fns Tin Tev Uw*2 

20010 IF K-0 THEN 
20020 PRINTER IS 701 
20030 PRINT 

20040 PRINT * Month, date and tine: * lOateS 

20050 IF Itt-0 THEN PRINT USING "10X,‘*Tube Type: 

20050 IF Itt-1 THEN PRINT USIN6 *10X,**Tube Type: 

20070 IF Itt-2 THEN PRINT USING *10X,**Tube Type: 

20080 PRINT 

20090 PRINT USING * 1 0X , * " Fns Tin Tev Uw*2 

20100 PRINTER IS 1 
201 10 K-1 
20120 ENO IF 
20130 8EEP 

20140 INPUT 'ENTER FLOUMETER REAOING ‘ ,Fms 
20150 OUTPUT 709; ‘AR AF0 AL4 UR 1 * 

20160 FOR L-0 TO 4 
20170 OUTPUT 709 i 'AS SA* 

20180 IF L>0 ANO L<4 THEN 20260 
20190 S-0 

20200 FOR 1-0 TO 9 
20210 ENTER 709<E 
20220 S-S+E 
20230 NEXT I 

20240 IF L-0 THEN En f < 0 )-A8S< S/ 1 0 > 

20250 IF L-4 THEN En f < 1 >-A8S< S/ 1 0 > 

20250 NEXT L 

20270 OUTPUT 709 i " AR AF00 AL00 URl* 

20280 OUTPUT 709 i *AS SA* 

20290 Etp-0 
20300 FOR 1-0 TO 9 
20310 ENTER 709<Et 
20320 Etp-Etp+Et 
20330 NEXT I 
20340 Etp-Etp/ 1 0 
20350 T in-FNTvsv( Enf ( 1 ) ) 

20360 Tev-FNT vs v ( Enf ( 0 ) > 

20370 Grad-37. 9853+. 104388*Tin 
20380 Mdot-3.9657E~3+Fns*< 3.61 9SSE-3-Fna* ( 8 . 82006E-6-Fna * < 1 . 23688E-7-Fns*4 . 3 1 897 
E-10) )) 

20390 Uu-Mdot / < 1000*PI*Oi'‘2)*4 
20400 Tdrop-Etp* 1 ,E+6/< 1 0*6rad ) 

20410 PRINT USIN6 * 10X ,3(00.00 ,4X ) ,1X ,2.00 ,4X , M2. 40* iFns .Tin, Tev ,Uw A 2 .Tdrop 
20420 BEEP 

20430 INPUT *UANT TO ACCEPT THIS OATA SET 7 ( 1 - Y ,0-N > * ,0k 

20440 J-J+1 

20450 IF Ok-0 THEN 

20460 J-J-t 

20470 GOTO 20130 

20480 ELSE 

20490 OUTPUT CFtleiFns ,Enf ( • ) ,Etp 
20500 PRINTER IS 701 

20510 PRINT USIN6 * 1 0X ,3( DO .00 ,4X ) . 1 X ,2 . 00 ,4X ,M2 . 40* iFna .Tin ,Tev ,Uw A 2 .Tdrop 
20520 PRINTER IS 1 
20530 8EEP 



T drop * ' ‘ 



Wieland Snooth* * * 
High Flux*** 
Tgrbo-B* * * 

T drop * ’ * 
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235^3 1%^’ -Ilw ’h£R.E 5E ANG t h£:; Gn~ A : r - - c r ^ » \ 3 - 

E 3553 I - 6 o__o ^ ■ ”hE\ £ 0 11 3 3 

20560 ENO IF 

20570 ASSI6N •File TO • 

205B0 PRINTER IS 701 
20590 PRINT 

20600 PRINT USING * 1 0X , ' "NDTE : ’ * . 22 . * * date eats ore stored in file "".ISA'iJ.F 
ilef 

20610 PRINTER IS 1 
20620 SU8END 
20630 SU8 Uoprt 
20640 PRINTER IS 1 
20650 BEEP 

20660 INPUT "Enter Uo File Nane’.FileS 
20670 BEEP 

206B0 INPUT "Number of Date Runs'.Nrun 

20690 INPUT "Oo You Went e Plot Fi le'?( 1 -Y ,0-N ) " , Iplot 

20700 8EEP 

20710 IF Io lot-1 THEN 

20720 INPUT "Give Plot File Name " .P_f i le* 

20730 CREATE BOAT P_file$,4 

20740 ASSIGN ©Plot TO P_fileS 

20750 END IF 

20760 PRINTER IS 701 

20770 PRINT 

20780 PRINT 

20790 PRINT USING "10X,"" Water Uel Uo""" 

20800 ASSIGN ©File TD Files 

20810 IF Iplot-1 THEN ASSIGN ©Filel TD P.fileS 
20820 FOP 1-1 TD Nrun 
20830 ENTER ©FilesUu .Uo 

20840 IF Iplot-1 THEN OUTPUT •Filel tUu.Uo 
20850 PRINT USING " 1 5X .0 . OD ,6X ,fl2 . 30E " \ Uw ,Uo 
20860 NEXT I 
20B70 ASSI6N ©File TD • 

20BB0 ASSIGN ©Filel TD ♦ 



20890 PRINT USIN6 " 1 0X . 


" "NOTE: 


" • .22 . 


,"" data sets ere stored 


in file " 


, 1 5A ’ 1 Nru 


n , F i 1 eS 












20900 IF Io lot- 1 THEN 












20910 PRINT USING "10X . 


" "NOTE: 


"" .22 , 


X-Y Pairs ere stored 


in file " " 


. 1 5A" ; Nru 



n ,P_f i le$ 

20920 ENO IF 
20930 PRINTER IS 1 
20940 SUBEND 
20950 SUB Select 
20960 COM /Ido/ I dp 
20970 8EEP 
20980 PRINTER IS 1 



20990 


PRINT 


USIN6 


" 4X , " " 


Select 00 1 ion 


2 ’000 


PRINT 


USING 


" 6X , " * 


0 


Taking date or re-proc e s s 1 ng previous date 


21010 


PRINT 


USING 


"6X ,"" 


1 


Plotting beta on Log-Log *’" 


21020 


PRINT 


USING 


*6X , " " 


2 


Plotting beta on Linear*" 


21030 


PRINT 


USIN6 


" 6X , " " 


3 


Me*e cross-plot coefft file""" 


21040 


PRINT 


USING 


’ 6X ," " 


4 


Re-circulete weter"" 


21050 


PRINT 


USING 


"6X ." " 


5 


Purge*"" 


21060 


PRINT 


USING 


"6X ." " 


6 


T-Droo correction" * 


21070 


PRINT 


USING 


"6X , " " 


7 


Pr 1 nt Uo File"" 


21080 


PRINT 


USING 


" 6X ,"" 


8 


Modify X-Y file"" 


21090 


PRINT 


USING 


" 6X , " " 


9 


Move" " 


21100 


PRINT 


USIN6 


"6X 


10 


Conb/F 1 xuD 


21110 


INPUT 


I do 









21120 IF Ido-0 THEN CALL Mein 
21130 IF Ido-1 THEN CALL Pio* 
21140 IF Ido-2 THEN CALL P 1 in 
21150 IF Ido-3 THEN CALL Coef 
21160 IF Ido-4 THEN CALL Main 
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Z'' - ’? I F lap *5 *mEN CAL_ "urg 
21163 IF lap-6 T riEN CALl. Tacn 

2H9® IF I dp-7 THEN CALL Uoprt 
21200 IF Ido-8 THEN CALL Xywod 
21210 IF Idp-9 THEN CALL Move 
21220 IF Ido-10 THEN CALL CoRb 
21230 SUBEND 
21240 SUB Xymod 
21250 PRINTER IS 1 
21260 BEEP 

21270 INPUT 'ENTER FILE NAME " .File* 

21280 ASSI6N •Fiiel TO File* 

21290 BEEP 

21300 INPUT "ENTER NUMBER OF X-Y PAIRS’ t Np 
21310 BEEP 

21320 INPUT ’ENTER NEW FILE NAME’,File2* 

21330 CREATE BOAT File2*,S 
21340 ASSI6N ®File2 TO File2S 
21350 BEEP 

21360 INPUT "ENTER NUMBER OF X-Y PAIRS TO BE DELETED’ .Ndel 
21370 IF Ndel-0 THEN 21410 
21380 FOR 1-1 TO Ndel 
21390 8EEP 

21400 INPUT "ENTER DATA SET NUMBER TO BE DELETED' ,Nd< I ) 

21410 NEXT I 

21420 FOR J-1 TO Np 

21430 ENTER »Filel iX ,Y 

21440 FOR 1-1 TO Ndel 

21450 IF Nd( I >-J THEN 21490 

21460 NEXT I 

21470 OUTPUT 0F l le2 i X ,Y 

21480 PRINT J ,X , Y 

21490 NEXT J 

21500 PRINTER IS 701 , 

21510 ASSIGN CF i lei TO ♦ 

21520 ASSI6N 0File2 TO • 

21530 SUBEND 
21540 SUB Move 
215501 FILE NAME: MOUE 
215601 

2 1570 01 M Bop (66 ) , A < 6 6 ) ,8( 66 ) , C( 66 ) ,0(66 ) ,E< 66 ) t F( 66 ) ,6(66 ) ,H(66 ) ,J(66 ) ,K(66 ) ,L( 
66 ) ,M( 66 ) 

21580 DIM Told* (66 >C 14] ,N(66) t Ur(66) , Ir( 66 ) 

21590 BEEP 

21600 INPUT -OLD FILE TO MOUE* ,D2_f lie* 

21610 ASSI6N *File2 TO 02_file* 

21620 ENTER *File2 iNrun ,0«te* ,Ldtcl ,Ldtc2 ,It t 

21630 FOR 1-1 TO Nrun 

21640 ENTER *F i le2 t8op( I ) ,Told*( I ) 

21650 ENTER *File2 |A( I ) ,B( I ) ,C( I > ,D( I ) ,E( I > ,F( I ) ,6( I ) ,H< I ) ,J( I ) ,K< I ) # L< I ) ,M( I >,N 
(I ) 

21660 ENTER •File2iVr( I > ,Ir( I > 

21670 NEXT I 

21680 ASSI6N *File2 TO • 

21690 BEEP 

21700 INPUT 'SHIFT DISK AND HIT CONTINUE" ,0k 
21710 BEEP 

21720 INPUT -INPUT BOAT SIZE' .Size 
21730 CREATE BOAT D2_f i le* ,S 1 ze 
21740 ASSI6N •File] TO D2_file* 

21750 OUTPUT *F 1 1 e 1 i Nrun ,Det e* ,Ld tc 1 ,Ldt c2 , 1 1 1 

21760 FOR 1-1 TO Nrun 

21770 OUTPUT •Fiiel ;Bop( I ) ,Told*( I > 

21780 OUTPUT fFtlel jA( I > ,d( I ) ,C( I ) ,D( I ) .E( I > ,F( I ) ,6( I > ,H( I ) .J( I ) ,K( I ) ,L( I ) ,M( I ) , 
N( I ) 

21790 OUTPUT #Filel iVp( I > ,Ir( I ) 
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2 1 800 nE*" : 

21810 ASSIGN S^iiel TO • 

21820“ RENAME * TEST * TO 02_filaS 

21330 BEEP 2000,. 2 

21840 BEEP 4000, .2 

21350 3EEP 4000, .2 

21860 PRINT *OAT A FILE MOUEO* 

21870 SU3END 
21880 SU8 Conb 
2 1 390 1 FILE NAME: COM8 
21900! 

21910 DIM Enf( 12) 

21920 SEEP 

21930 INPUT *OLD FILE TO FIXUP* ,D2_f i le* 

21940 ASSIGN ®File2 TO 02_fileS 
21950 OI_f ile*-*TEST" 

21960 CREATE 80AT Ol_file*,20 

21970 ASSIGN ©F i la 1 TO 01_filaS 

21380 ENTER Ie2 iNrun .OateS ,Ldtc 1 ,Ldtc2 , It t 

21990 Nrunw-20 

22000 IF K-0 THEN OUTPUT 0F i le 1 i Nrunw .Data* ,Ldtc l ,Ldtc2 . 1 1 1 
22010 p OR 1-1 TO Nrun 

22020 ENTER #Fi le2;8op .ToldS ,E«f ( • ) ,Ur , Ir 
22030 OUTPUT 3F i la 1 1 8op .ToidS .E«f ( * > ,Vr . I r 
22040 NEXT I 

22050 ASSIGN ©Fila2 TO • 

220601 RENAME ‘TEST * TO D2_fileS 

22070 8EEP 2000 ,.2 

22080 SEEP 4000,. 2 

22090 8EEP 4000. .2 

22100 8EEP 

22110 INPUT * WANT TO AOO ANOTHER FILE ( 1 -Y ,0-N >*> " .Oka 
22120 IF Oka-1 THEN 
22130 K-1 
22140 3EEP 

22150 INPUT *6IUE NEU FILE NAME* ,Nf i laf 
22160 ASSIGN ©F i I e2 TO NfilaS 
22170 GOTO 21980 
22130 ENO IF 

22190 ASSI6N ©Cile2 TO • 

22200 SU8EN0 



Oa t a : 1 9 Jan 1 993 

NOTE: Program nawa : DRP6B 
Disk nuwbar - 00 
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